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INTRODUCTION 


An experiment in Electromyography (EMG) utilizing surface electrode tech- 
niques was developed for the Apollo-Soyuz Test Project. The objective of the 
study was to obtain quantitative measure of changes in muscle characteristics as a 
result of exposure to weightlessness (see Appendix A). This report describes the 
computer program, EMGAN, which was written to provide first order data reduc- 
tion for the experiment. 

EMG signals are produced by the membrane depolarization of muscle fibers 
.during a muscle contraction* Surface electrodes detect a spatially summated 
signal from a large number of muscle fibers commonly called an interference 
pattern. An interference pattern is usually so complex that analysis through signal 
morphology is extremely difficult if not impossible. However, because of the ease 
of application and noninvative character of surface electrodes and the often de- 
sirable feature of studying large groups of muscle fibers, methods have been sought 
to evaluate the surface interference pattern. 

In the time domain a valuable technique to characterize an interference pattern 
has been the root-mean-square (RMS) value of the signal. A similar technique has 
been to integrate the interference pattern over a given time length. A linear rela- 
tionship is usually found between the integrated EMG or the RMS value and the force 
exerted by the muscle up to some maximum. 

More recently with the increased access to the digital computers and the avail- 
ability of an algorithm to economically calculate discrete Fourier transforms, it has 



become common to process EMG interference patterns in the frequency domain. 
Muscle fatigue and certain myopathic conditions can be recognized through changes 
in muscle frequency spectra. 

DATA PREPARATION 

The protocol for the AS TP Myography required that two types of data be record- 
ed, digitized, and processed: the force being exerted by a muscle group and the 
surface EMG signal. The data was originally recorded on a 4-channel Tandberg 
analog recorder with one-channel of IRIG 0 time code, one-channel of force data, 
and two-channels of simultaneous EMG signal (two different muscles). 

A digital tape was prepared from the analog data using the General A/D program 
(see Appendix H) available on the CF-16A mini-computer located in the Sigma 3 
Computer Room. Analog force data was low pass filtered at 10 Hz before being 
digitized at 20 samples/second. Analog EMG data was low pass filtered at 400 
Hz. Consistent with this frequency limit a digitizing rate of 1000 samples/second 
was chosen as the minimum necessary for power spectral density analysis 12). 

The digitized data was organized into half second records with 535 words per 
record. The first twenty-five words were header information with ten words of 
digitized force data interspaced into 500 words of EMG data (see data format in 
Appendix B). Only one-channel of EMG data was digitized with its associated 
force data at any one time. Although not required by EMGAN, a data tape was 
usually constructed with the first file containing force calibration and subsequent 
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files containing EMG data from each muscle with the appropriate EMG calibration 
signals. 

PROGRAM CHARACTERISTICS 

EMGAN exists in two versions. One version was written to execute on the 
Biomedical Sciences Division Xerox Sigma 3 Computer using the Computek 400/20 
CRT terminal for graphical output. A second version was written to execute on 
the Cardiovascular Laboratory's Data General Nova 820 mini-computer. Both 
programs are essentially identical to the user. However, there are a number of 
programming differences made necessary by differences in Xerox and Data 
General's Fortran IV, supporting software packages (graphical and tape I/O), 
and peripheral hardware. Whenever these differences are sufficiently important 
in describing the program they will be noted. Complete program listings of both 
versions with their associated load maps are given in the Appendices F and G. 

An EMG system tape exists for the Sigma 3 Computer which contains an 
appropriately sized operating system. EMGAN exists in absolute binary as a file 
in the User Processor area of the Rad. A source listing of EMGAN exists both on 
magnetic tape stored in the Sigma 3 Computer Room and on punched cards in the 
Cardiovascular Laboratory. Both source listings and absolute binary files of the 
Nova version of EMGAN are stored on disk pack and on magnetic tape in the 
Cardiovascular Laboratory. The absolute binary files are named EMGAN. SV and 
EMGAN. OL (both are needed for program execution). The source listing is lo- 
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cated in separate files for the main program and each subroutine with each file 
named for the subroutine it contains. How to execute a run with either computer 
version is explained in a later section. The computer peripherals required to 
fully execute CMGAN are a card reader, line printer, two 9-track tape drives, 
and the Computek terminal (Sigma) or the Tektronix 4014 terminal (Nova). 

EMGAN's execution is totally controlled by the data card information receiv- 
ed. Force and EMG data slices are processed with separate data cards. The 
maximum length of a force data slice is 120 seconds with previous force calibra- 
tion required. After the raw data is converted into pounds it is averaged over a 
user selected interval for output to the line printer. The maximum length of an 
EMG data slice is four seconds with previous EMG calibration required. After 
scaling into microvolts, an integrated value is found for the EMG data slice. The 
discrete Fourier transform of the data is taken. The raw power spectral density 
(PSD) of the data up to 400 Hz is calculated and smoothed over a frequency band- 
width selected by the user. The raw PSD with associated header information is 
outputed to magnetic tape for use in any subsequent processing program. The 
smoothed PSD is outputed to the line printer along with accessory information 
such as a normalized PSD, percent contribution of each bandwidth to the whole 
spectrum, the mean and standard deviation of the spectrum, and time domain 
integration and mean square value of the data slice. 

EMGAN run time is variable and depends on the number and kind of data 
cards, the length of the EMG data slices, whether or not the 60 Hz filter is used. 
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data location on the input tape (processing data in sequential order on the tape 
is much faster than skipping back and forth on the input tape), the number of 
data plots, and whether execution is on the Sigma or Nova. For the ASTP 
E MG experiments a typical run consisting of calibration cards, 16 force data 
slice cards, and 28 4-second EMG data slice cards (no filtering or plotting) re- 
quired approximately 25 minutes on the Sigma and 35 minutes on the Nova. 

EMGAN DESCRIPTION 

In describing the various subroutines, differences in the Nova version from the 
Sigma version will be set off with a ***. A table of formulas used in the program 
is listed in Appendix D. 

Main Program; 

The Main Program's basic task is to load the needed overlays in the proper 
sequence. A flowchart is given in Appendix C. The program is organized into a 
large do loop with one pass through the loop for each data card to be processed. 

A check is made to determine if calibration data has been received before the 
respective data type, force or EMG, can be processed. If the user should request 
the processing of data without calibration, the program outputs an error message 
and halts. 

*** The magnetic tape units are initialized and the output tape file is open- 
ed. Upon completion of the program the tape units are released. 


Subroutine CINPUT 


Subroutine CINPUT inputs all data card information. Data card format is 
given in Appendix B. As each card is read it is outputed to the line printer to 
aid the user in verification. If the number of data cards to be processed is read 
to be greater than 200, an appropriate error message is given and the program 
halts. All data card information is stored in a scratch file, ICARD, in the user 
data area of the RAD. As each data card is to be processed ICARD is read back 
into main memory. The program variables are then set to the corresponding data 
card values. 

*** The scratch file, ICARD, is created on the disk pack by the Nova pro- 
gram if it does not already exist. For the Sigma version the user should create 
ICARD (1 record, 6000 bytes, random format) if it does not already exist before 
attempting to execute EMGAN. 

Subroutine TINPUT 

Subroutine TINPUT acquires the digitized data from the input tape. Tape 
files are numbered from 00 to 99. The routine keeps track of the current file 
being accessed. If a data card requests another file, TINPUT skips the tape 
either forward or backward to the beginning of the requested file. Subroutine 
FIND is called to position the tape at the start of the data slice within the file. 
TINPUT reads either force or EMG data according to data card request. After 
the data has been read, the tape is backed up to the beginning of the data slice. 
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If an EOT is encountered during a search for a requested file, an error message is 
given and the program halts. If an EOF or EOT is encountered during the input of 
a data slice, a error message is given and the program halts. 

*** Nova Fortran calls for magnetic tape I/O are considerably different from 
the QlNOUT package used on the Sigma. Because Nova error checking is more 
comprehensive, error messages given by EMGAN resulting from tape l/O problems 
contain much more information. 

Subroutine FIND 

Subroutine FIND searches the currently accessed tape file for the start time of 
the data slice. If an EOF is encountered before the start time is found, the tape 
is rewound to the beginning of the file. A second search is then made through the 
entire file. If the start time is still not found, an error message is given and the 
program halts. 

Subroutine QUT60 

Subroutine OUT60 is a digital notch filter which can be called to remove 60 
Hz and its odd harmonics from EMG data. If at all possible filtering should be 
avoided, since no distinction can be made between 60 Hz present in the true data 
and that which may have been picked up from extraneous sources. Since a large 
number of sines and cosines are generated, considerable time is added to the pro- 
cessing of data slices. 
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Subroutine DCAL 


Subroutine DCAL processes force and EMG calibration data to obtain the 
scale factors needed to change the data from raw counts to engineering units. 
After the scale factors have been calculated, the calibration is converted to 
engineering units for plotting. 

A force calibration data slice consists of two DC levels, the first level corre- 
sponding to zero pounds of force, and a second level corresponding to a known 
amount of force (98 pounds in the ASTP experiment). The data slice is searched 
for a jump in the data to indicate the beginning of the second level. When the 
jump is found, the average number of counts three seconds before the jump is 
subtracted from the average number of counts three seconds after the jump. The 
difference is divided by the number of pounds represented by the jump to obtain 
the scale factor. If the jump is not detected, an error message is outputed and 
the program halts . 

The EMG scale factor is calculated using the fact that for a zero mean, 
stationary signal the standard deviation of the signal is equal to its RMS value. 
(For the ASTP experiment a nominal 20 Hz, 350 microvolt RMS signal was used 
for EMG calibration). The calibration data first has any DC mean subtracted 
from it and the standard deviation of the data in counts is found. The scale factor 
is computed by dividing the standard deviation by the RMS magnitude of the 
cal ibration signal . 
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Subroutine GRAPH 


Subroutine GRAPH plots on the Computek terminal, if requested by the data 
card, force and EMG calibration data, force data, and EMG data both in the 
time and frequency domain. Scaling for force and EMG data is variable accord- 
ing to the maximum amplitude of the data. Force data is plotted without smooth- 
ing (such as is done for the printed output). When time domain EMG data is 
plotted, only 1000 evenly spaced points are plotted no matter what the length of 
the data slice. The PSD plotted is smoothed and normalized to the maximum value 
of the spectrum. The return key must be struck for the program to continue when- 
ever the terminal's bell is rung. 

***As of the writing of this report the graphics package for the Tektronix 
4014 terminal is not available. No plotting capability presently exists when 
executing EMG AN on the Nova. However, if plots are requested when executing 
on the Nova, no run error occurs. 

Subroutine MODLINE, WORDS, and BELL 

These routines refer to the Computek terminal and respectively, plot data 
arrays, output alphanumeric characters, and ring the Computek bell. 

Subroutine EMG 

Subroutine EMG ubtracts any DC offset from the unsealed EMG data. The 
data is then scaled into engineering units and the largest absolute value of the 
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data array is found for use in plotting of the data. Simpson's rule is used to find 
the integrated value of the data in units of microvolt-seconds. The mean square 
value of data is calculated. 

Subroutine FORCE 

Subroutine FORCE scales the force data into engineering units. Compensation 
is made for force signals recorded at a higher gain than the force calibration. 

The first second of force data is assumed to be the zero force level • The user 
should insure that this is so for correct results. The maximum force in the data 
slice is found for use in plotting the data. Force data is then averaged over 
intervals requested by the data card for output to the line printer. 

Subroutine FFTPSD 

Subroutine FFTPSD calculates the raw power spectral density of an EMG data 
slice up to 400 Hz. It is suggested that the user read (2) for information on data 
analysis using PSD techniques and (1) for an excellent explanation of discrete 
Fourier transforms and the fast Fourier transform algorithm. A cosine taper is 
applied to the first and last tenths of the data slice to reduce leakage. The 
Fourier transform of the data is then taken and the raw power spectral density is 
computed. The PSD amplitudes are corrected for reduction caused by the cosine 
taper. The raw PSD is outputed to magnetic tape for use in subsequent processing 
programs. The PSD is smoothed over the frequency bandwidth requested by the 
data card and the maximum PSD value found. A second PSD array is computed. 
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normalized to the maximum PSD value. The area under the PSD is found and used 
to calculate the percentage each bandwidth contains of the total power and the 
cumulative of percentage of total power with increasing frequency. The expected 
value and variance of the PSD are computed. 

Viewed statistically the raw PSD has very poor reliability. Each calculated 
value is an inconsistent estimate of the true value and has a possible random 
error of 100%. Smoothing the PSD can greatly reduce the random error and pro- 
vide a better estimate over the smoothed bandwidth (2). The normalized standard 
error is calculated to provide information on the amount of random error associated 
with each PSD value. 

Subroutine WINDOW 

Subroutine WINDOW applies a cosine taper to the first and last tenths of the 
data slice to reduce leakage in the discrete Fourier transform. 

Subroutines RFORT and FORT 

Subroutine RFORT and FORT take the Fourier transform of the EMG data. FORT 
contains the actual fast Fourier transform algorithm and assumes a transform of N 
complex data points. RFORT compensates for the actual case of a transform of 
2N real data points. The FORT computes an initial sine/cosine array and refers to 
this array when computing the transform. Repetitive computation of identical sines 
and cosines is thereby avoided, substantially increasing the speed of the transform 
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at the cost of an extra array of size N/4. Both of these subroutines were received 
from Mr, Jack McBryde of Lockheed. The only modification was to not divide 
each data value by the number of data points (an alternate definition of the dis- 
crete Fourier transform). EMGAN can transform 1, 2, or 4 seconds of EMG data 
which corresponds to 1024, 2048, and 4096 data points, respectively. Computa- 
tion time for the Fourier transform algorithm are as follows: 


of Points 

Sigma Time 

Nova Time 

1024 

6 seconds 

10 seconds 

2048 

13 seconds 

22 seconds 

4096 

29 seconds 

47 seconds 


*** Subroutines RFORT and FORT are named RFFT and FFT to avoid confusion 
with Nova Command Line Interpreter calls. 

Subroutine PSDSAVE 

Subroutine PSDSAVE outputs to magnetic tape the raw PSD of each EMG data 
slice. The PSD values are preceded by a 25-word header record. The format of 
the header recc'd is located in Appendix B. The PSD of a one-second slice of 
data has 512 real values, a two-second slice has 1024 real values, and a four- 
second slice has 2048 real values. A PSD record contains 2048 integer words; 
therefore, a one-second PSD fills 1/2 of one record, a two-second PSD fills one 
full record, and a four-second slice fills two records. Two EOF’s are written on 
the tape after the raw PSD is output and the tape then backs up to just before the 
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EOF's. At the end of an EMGAN run the output tape consists of a file containing 
all of raw PSD's from the run terminated by two EOF's. 

Subroutine PRINT 

Subroutine PRINT outputs to the line printer all of the processed force and 
EMG data. A page of header and calibration information is printed whenever 
the print switch on a calibration data card is set. The length of a force or PSD 
output is checked and if sufficiently long a two column list of the processed data 
is printed instead of a one column list. The frequency corresponding to each 
PSD value is the center of the bandwidth. 

*** Subroutine PRINT is named POUT to avoid confusion with Nova Command 
Line Interpreter calls. 

SAMPLE EMGAN RUN 

To execute an EMGAN run on the Sigma 3 first load the EMG system tape. 
Place the data cards in the card reader and power up the Computek terminal . 

Mount the input data tape on unit 0 and the PSD output on unit 1 . Assign the 
user input device to the teletype and type in "IEMGAN" . 

To execute an EMGAN run on the Nova load the EMG disk pack and bring 
up the system. Place the data cards in the card reader. Mount the input data 
tape on unit 0 and the output tape on unit 1. Type in "EMGAN". A message 
will be outputed to the terminal, ' LOAD $CDR, STRIKE ANY KEY", upon 
entering any keyboard character, the data cards will be read and the program 
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executed. 


An example of an EMGAN input card deck with the printed and plotted out- 
put is in Appendix E. 

PROGRAM IMPROVEMENTS 

There are several additions that would be desirable in future work with 
EMGAN. The major addition would be a subroutine to check the EMG data 
slice prior to transform for stationarity and normal distribution. Suggestions for 
the necessary statistical techniques are found in (2). 

An alternative FFT algorithm might be used which calculates sines and cosines 
as needed rather than establishing a table. The core storage saved would allow 
an 8196 real data point transform to be performed on the Sigma and, with proper 
sizing of operating system, on the Nova. However a price would be paid in 
increased processing time considerably beyond twice the transform time of 4096 
real data point array. 

Of a more minor nature subroutine GRAPH could be modified to output the 
muscle name with the EMG and PSD plots to allow easier later identification. 
Although force and EMG data are now processed separately, it could be useful to 
obtain the average force value over the period of the EMG data slice. The force 
value could be output with the EMG‘s PSD fo more convenient association be- 
tween the two . 
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APPENDIX A 


AS TP MYOGRAPHY DESCRIPTION 



1 .0 General Objective 


The general objective of the experiment is to continue the study efforts begun 
in the Skylab program to identify and describe antigravity muscle dysfunction 
characteristics and consequences resulting from spaceflight. 


2 .0 Specific Objective 

The purpose of the experiment is to assess changes that occur following a 
period of disuse. Some studies suggest that the first few days of exposure to 0-g 
may be significant to the ultimate muscle deconditioning resulting from longer 
missions. The relationship between muscle capability, in terms of strength or 
tension, fatigability, and muscle electrical activity will be investigated, as well 
as the differential effects of spaceflight disuse on "fast" and "slow" muscles. 

3 .0 Significance 

Muscle function and condition may well be a critical determinant of man's 
capacity to endure the effects of long duration weightlessness as well as the 
readaptation to the earth's gravity, or to the gravity of other planets. Data 
collected in the experiment will aid in quantifying the muscle deconditioning 
which results from weightlessness. 

In addition, an important spinoff will be the extension of an already con- 
siderable ground-based body of knowledge about the characteristics and conse- 
quences of muscle disuse. Heretofore, research in this area has relied on the 
"contrived" methods of surgical section and limb or torso fixation to produce the 
disuse effects. Spaceflight uniquely provides a “pure" form of neuromuscular 
system disuse . 


4.0 Method 

A . Concept 

The measurement of muscle electrical activity is known as electromyography 
(EMG). The state of muscle function can be described by EMG measurements com- 
bined with knowledge of the force being exerted by the muscle. 

Electromyographic studies have shown that skeletal muscle undergo chat ges 
in capability and composition when subjected to periods of disuse, i.e., periods 
of time when the contractile mechanisms are not subjected to the stresses and 
forces normally encountered. Also changes in biochemical constituents, as seen in 
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Skylab, such as calcium, potassium, etc., and alterations in enzymic constiuents 
such as ATP, and acetylcholinesterase have been shown, by ground studies to 
affect normal muscle function. 

The EMG results of Skylab 3 skeletal muscle assessment provide ample 
evidence that normal muscle function is significantly altered by periods of weight- 
lessness (disuse) of 56 days or more. 

To investigate the effects of a shorter period of weightlessness on muscle 
function, a standardized test protocol and a muscle stress device will be used. 

The muscle stress apparatus will provide for preplanned isometric muscle forces 
from the calf muscles and the arm muscles. The standardized test protocol will 
include measures of muscle strength, muscle endurance, and muscle fatigability. 
Pre- and postflight measures only will be taken. 

B. Procedure 

The test procedure is identified in Figure 1. Two muscles each from the leg 
and arm will be instrumented and the muscle action potentials recorded. These 
muscles are the brachial bicep, brach ioradial is, gastrocnemius and soleus. The 
procedure will require about 11 minutes to complete. 

C. Data 

The EMG data will be recorded on magnetic tape for time and frequency 
domain analysis at the JSC Cardiovascular Laboratory. 
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PROCEDURE 


A. Apply surface electrodes 

B. Seat subject in muscle stress device prepared for muscle stress test 

C. Three short (1 to 2 seconds) efforts to determine the maximum voluntary con- 
traction (MV C) 

D. Remainder of calf procedure 

1 . 10 seconds at 10% MV C - 20 seconds rest 

2. 10 seconds at 20% MVC - 20 seconds rest 

3. 10 seconds at 30% MVC - 20 seconds rest 

4. 60 seconds at 50% MVC 

E. Reset muscle stress device for arm stress test 

F. Three short (1 to 2 seconds) efforts to determine the maximum voluntary 
contraction 

G. Remainder of arm procedure 

1 . 10 seconds at 10% MVC - 20 seconds rest 

2. 10 seconds at 20% MVC - 20 seconds rest 

3. 10 seconds at 30% MVC - 20 seconds rest 

4. 60 seconds at 50% MVC 

H. Remove surface electrodes 
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APPENDIX B 


INPUT/OUTPUT MAGNETIC TAPE AND DATA CARD FORMATS 



EMG DATA INPUT TAPE RECORD 




Hours 


3 


Minutes 


Seconds Mili- | Subject Exp. 

Seconds No. Month 




Digital 

Month 



Digital 

Year 


13 14 


Flight Run 
Ref. Day No. 


15 


Analog 

Tape 

No. 



19 


Not Used 



ill 

Not Used 



Not Used Not Used 


EMG Dcfa 
Data Samp l Scnl p| e 


231 

EMG 

Data 

Semple 


29 I 

”EMG 

Data 

Sample 



Data 

Sample 



129j 

130 

Fo-’ce 

EMG 

Dcra 

Dcta 

Scmole 

Sample 



Force 

Data 

Sample 



132 ] 
EMG 


Data 

Sa.mole 

! 





ml 

U- . 

Force 

EMG 


Data 

Data 


Sample 

Sample 




EMG 

Data 

Semple 













































HEADER RECORD FOR PSD OUTPUT 


Wotd * 


1 - Subject Number 

2 - Experiment Month 

3 - Experiment Day 

4 - Experiment Year 

5 - Digitizing Month 

6 - Digitizing Day 

7 - Digitizing Year 

8 - Flight Reference Day 

9 - Run Number 

10 - Analog Tape Number 

11 - Sample Rate for EMG Signal (Samp/sec) 

12 - Sample Rate for Force Signal (Samp/sec) 

13 - RMS Amplitude of Sine Wave Cal (in microvolts) 

14 - Not Used 

15 - Amount of Force Cal 

16 - Start Hour 

17 - Start Minute 

18 - Start Second 

19 - Length of Data Slice (sec) 

20 - Muscle ID Number 

21 - ^ of Following Records with PSD Values 

22 - 25 - Not Used 
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DATA CARD FORMAT 


first Card In Data Card Deck 

cc: 1-5 Number of data slices (max = 200) 

Data Slice Cards 

cc: 1-3 File data slice located in (00 to 99) 


cc: 4-6 Data Type 

1 = Force Cal 

2 = EMG Cal 

3 = EMG Data 

4 = Force Data 


cc: 7-9 HR start time of data slice 

cc: 10-12 Min start time of data slice 


cc: 13-15 
cc: 16-18 


cc: 19-21 


cc: 22-24 


cc: 25-27 


Sec start time of data slice 

Filter switch 
-1 = No filter 

0 = Filter 60 Hz 

1 = Filter 60, 180 and 300 Hz 

Length of data slice in seconds 

Force data or cal - any integer number up to 120 seconds 
EMG Cal - 1, 2, 3 or 4 seconds 
EMG Data - 1, 2, or 4 seconds 

Plot switch 

0 = Plot data 

1 = Suppress plotting 

Print switch 

0 = Print data 

1 = Suppress printing 
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cc: 28-30 Averaging interval (force and EMG data cards only) 

Force Data = over any number of seconds up to the length 
of data slice 

EMG Data = (Frequency smoothing only) any integer number 
up to 400 

cc: 31-33 Force signal gain with respect to force calibration (force 

cards only 

cc: 34-36 Muscle Type (cal data cards only) 

1 = Brachial Bicep 

2 = Brachioradialis 

3 = Gastrocnemius 

4 = Soleus 
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appendix c 


main program flow chart 




































EliGAH Flowchart 
Figure IE 

C-' 






APPENDIX D 


TABLE OF FORMULAS 



TABLE OF FORMULAS 


1 • Force Calibration Variable (counts/ pound) 

FORSLP = (H/20 - L/20)/FLBCAL 

where H = sum of samples for one second of force calibration level 

L = sum of samples for one second of zero force level 
FLBCAL= amount of force calibration in pounds 

2. EMG Calibration Variable (counts/microvolt) 

EMGAL =J y^X^/(NPTS-l)/MAGTUD 

where = ^ um va ^ ues squared 

NPTS = number of EMG cal signal samples 
MAGTUD = RMS voltage of EMG cal signal 

3. Mean Square (Variance) 

EMGVAR /(NPTS-1) 

2 

where ]>^X. = sum data values squared 
NPTS = number of data points 

4. Simpson's Rule (without error term) 

VOLTSEC = H/3 (2SUMO+4 SUME - Data (1) + Data (n)) 
where H = spacing between the data values 
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SUMO = sum of the odd numbered data values excluding first data value 
SUME = sum of the even numbered data values excluding last data value 
Data (1)= first data value 
Data (n)= last data value 

5. Scaling of EMG Data (microvolts) 

SEDATA(i) = (USEDATA(i) - OFFSET)/EMGCAL 
th 

where USEDATA(i) = i unsealed data value 

OFFSET = DC value of EMG data array 
EMGVAL = EMG calibration value 


6. Scaling of Force Data (pounds) 

SFDATA(i) =(USFDATA(i) - BASE) * GAIN/FORSLP 
th 

where USFDATA(i) = i unsealed force data value 
BASE = zero force level 

GAIN = ratio of force calibration gain to force data gain 
FORSLP = force calibration variable 


7. Discrete Fourier Transform 

N-l 


nk 


X(n) = y x(k) (EXP (-j2 VN)) n = 0, 1, ... .N-1 


K=0 


th 


where x (k) =k value of untransformed data 
n = number of data points 
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8. Power Spectral Density (magnitude) 

PSD (n) = 8H/NPTS [Re(X(n)) 2 + im (X(n)) 2 ) 

where H = spacing between data samples 

NPTS = number of data samples 
th 

Re(X(n)) = real part of n Fourier transform value 

th 

lm(X(n)) = imaginary part of n Fourier transform value 


9. Expected (Mean) Value of PSD 


EXPVAL =V* F (PSD(n) - PERCNT(n)/100) 


where PSD (n) = value of n^ 1 power spectral density bandwidth 

th 

PERCNT (n) = percent of total power contributed by n bandwidth 

th 

F = center frequency of n bandwidth 


10. Standard Deviation of PSD 


STDEV =JX/ F ' EXPVAL)2 * ( pERCNT ( n )/ 10 °) 
where EXPVAL = expected value of PSD 


th 


PERCNT (n) - percent of total power contributed by n bandwidth 

th 

F = center frequency of n bandwidth 
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APPENDIX E 


EXAMPLE OF EMG AN INPUT CARD DECK WITH PLOTTED 


AND PRINTED OUTPUT 




* * ★ ★ DATA CAHJS * * * * 

NO. OF DATA SLICES= 11 
DATA SLICES 


00 

1 

17 

41 

0 -1 

30 

0 

1 

0 

1 

0 

0 

01 

i 

17 

14 

22 -1 

4 

0 

0 

0 

0 

3 

0 

01 

4 

17 

lb 

le -1 

20 

0 

u 

1 

4 

0 

0 

Cl 

J 

17 

1 6 

22 -i 

4 

0 

0 

10 

0 

0 

0 

01 

a 

17 

lb 

55 -1 

20 

0 

0 

1 

2 

0 

0 

01 

3 

17 

17 

1 - i 

4 

0 

0 

10 

0 

0 

0 

01 

a 

17 

17 

35 -1 

70 

0 

0 

1 

1 

0 

0 

01 

3 

17 

17 

41 -1 

4 

0 

0 

10 

0 

0 

0 

01 

3 

17 

17 

57 -1 

4 

0 

o 

10 

0 

0 

0 

01 

3 

17 

18 

14 -1 

4 

0 

0 

10 

0 

0 

0 

01 

3 

17 

16 

31 -1 

4 

0 

0 

10 

0 

0 

0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
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CARDIOVASCULAR LARUKATUR> 
EMG DATA PROCESSING PWuoRAM 


*** header INFORMATION*** 


SUBJECT NO.! 

44 

EXPERIMENT DATE: 

6/28/75 

RUN NO.: 

3 

analog tape NO.: 

9 

MUSCLE? 

GASTROCNEMIUS 



EMG SAMPLE RATE 

(SAMP/SEC): 1000 

FORCE SAMPLE RATE 

(SAMp/SEC) : 


FLIGHT REFERENCE DAY: F-15 

DIGITIZING DATE: 6/29/75 


FORCE CAL 

CAL WEIGHT (LBS): 98.0 

EMG CAL 


** CALIBRATION ''ATA ** 


COUNTS/LB: 39. AVG. BASELINE COUNT: -3662. CAL GAIN: l. 


RMS AMPLITUDE (MICROVOLTS): 350 


COUNTS/MICROVOLT: 6.44 



*** FORCE DATA 


*** 


Ti«t- 17: Is: lb 


DATA LEMGTh (SECS)- £0 AVERAGE FCRCt INTERVAL (SECS)- 1. 

4 

interval average force (lbsj 



2 

3 

a 

5 

6 
7 
6 
9 

10 
1 1 
12 
13 
la 

15 

16 
17 
16 

19 

20 


0.0 
0.0 
0.0 
0.1 
27 .7 
a 1 . 1 
<*1.9 

41.5 
41 .4 
aO.9 
a2.5 

42. 2 

al.3 
a 1 . to 

42.3 

41.6 

21.7 
0 . 0 

- 0.0 

o.c 




*•* POkER SPECTRAL DENSITY Oft b*G JAT A *** 


START TIME- 17:16:22 
BANDWIDTH (HZ)- 9.766 
MEAN (HZ)- 199.2 


DATA LENGTH (SECS)- 4 

4 

NORMALIZED STANDARD ERROR- U. 15 * 
STANDARD Ot V 1 AT IL'N (hZ ) - 60 . ft 


IMtGWAltU EMG lMCROVOLT*SEC)- 
lNrtl>wAltD PSD (V1LR0V0LT**2)- 
EMG VARIANCE (VlL*UVULT**2)- 


FRED 


PSD 

PSD 

X OF 

CUM X 

(HZ) 


(MMV W/HZ) 

NO»y 

TOTAL 

total 

4.8b 


1.39ft 

0.0306 

0.1 7 

0.17 

14. b5 


6 . 7b6 

0. 10A? 

0 .84 

1.01 

29. ai 


15.152 

0.3319 

1 .69 

2.89 

39.ib 


23. 153 

0.507? 

2.87 

5.75 

43.95 


16.556 

0.40*5 

2.30 

8.05 

5 3.71 


30.225 

0 . 66? 1 

3.74 

11.79 

63.4b 


34.020 

0.7052 

4.21 

16.00 

73.20 


27.169 

0.5951 

3.36 

19.36 

b.i.01 


95.652 

1.0000 

5.65 

25. U1 

92.77 


32.585 

0.7138 

4.03 

29.05 

102.59 

9 

39.069 

0 . 860 o 

4.89 

33.93 

112.30 


92.109 

0.9220 

5.21 

39.14 

122.07 


90.910 

0.69fc2 

5.06 

44.21 

131.89 


95.400 

0.99u5 

5.62 

49.83 

101 .60 


39.921 

0.8635 

4 . 68 

54.71 

151.37 


43.597 

0.9550 

5.40 

© 

o 

o 

161.13 


36.901 

o.eios 

4.5ft 

64.68 

170.90 


39.035 

O.7055 

4.21 

e>8.89 

1 60 • ob 


20.588 

0 . o5 1 0 

2.55 

71.44 

190.03 


21 .723 

0 . o 7 56 

2.69 

74.13 

200.20 


23.019 

0.5CO1 

2.85 

76.98 

209.96 


27.384 

0.599ft 

3.39 

80.37 

219.73 


20.216 

0 • 04 28 

2.50 

82.87 

229.49 


15.688 

0.30 36 

1.90 

84.81 

239.26 


17.412 

0.381O 

2.16 

86.97 

299.02 


13.022 

0.2900 

1 .6b 

88.63 

258.79 


9.829 

0.2153 

1.22 

89.84 

266.55 


10.397 

0.2277 

1.29 

91.13 J 

276.32 


9.821 

0.2151 

0.1317 

1.22 

92.35 

266.09 


6.011 

0.7a 

93.09 

297.85 


8.662 

0.1697 

• 1.07 

94.1b 

307.62 


7.959 

0.1743 

0.99 

95.15 

317.33 


6.416 

0.1U05 

0.79 

95.94 

327.15 


5.909 

0.1294 

0.73 

9b. b7 

336.91 


7.114 

0.155ft 

0.88 

97.55 

346.68 


5.590 

0.1224 

0.69 

98.25 

356.05 


3.670 

0.0803 

0 • 4 r . 

98.70 

366.21 


4.024 

O.Oftftl 

0.50 

99.20 

375.93 


2.570 

0.0564 

0.32 

99.52 

365.79 


2.234 

0.0009 

0.28 

99.79 

395.51 


1.670 

0.0366 

0.21 

100.00 


0.2676E 3 

0.7890E 4 

0.6019E U 





17:14:22 17.14:24 



1CC9 21*3 J CJ3 ***3 SSS 

rir>i ia HiLUseccscs 

exp OAT* 4/23/75 RWf« 5 








33 # 44 


ftUtf !> 



E*C DATA 


1 £T3 1CC3 1 

133 

TG~* 2Z & 3 


UC3 

<5533 


Tins 

in iiiLureccNOs 




3V>J 44 

EXP 

OATS 6/23/7S 


JtUXJ 3 
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APPENDIX F 


SIGMA SOURCE LISTING AND LOAD MAP 



DJOB 

DASS U0*6 
DEDITOR 
-LIST 
-GEN 

cxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxtxxxxxxxxxxxxxxxxxxxxxxxxtxxxxxxxxxxxxxxtxxx 

c 

c 

C PROGRAM EMGAN 

C 

C AUTHOR UILLIAM N HURSTA 

C 

C PURPOSE TO PROUIDE FIRST ORDER REDUCTION OF ENG AND FORCE DATA 

C 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxzxxxxxxxxxxxxxxxxxxxx 

c 

c 

DIMENSION IVARSU(10) 

INTEGER FILNOU , DTVPE . FILNUM -FILTSU.PLOTSU 
INTEGER FCALSU ,ECALSU,UHOA,PRNTSU 
COMMON DATA ( 8590 ) , I HEAD 09 ) . CAL( 2. 4 ) 

EQUIVALENCE (DTVPE- IUARSUC 1 ) ) , (FILNUM . IUARSU( S ) ) 

EQUIVALENCE (FILTSU . IVARSLK 3 ) ) . (NSECS. IVARSU(4)) 

EQUIVALENCE (PLOTSU. IVARSU( 5 ) ) . (PRNTSU . IVARSUt 6 ) ) 

EQUIVALENCE ( ISPAN. IVARSU(7 ) ) . ( IFGAIN, IVARSUC3 ) ) 

FILNOU-1 

FCALSU-9 

ECALSU-9 

UHOA-0 

NSLICE-9 

C 

c GO GET NUMBER OF DATA SLICES FOR PRESENT RUN 

C 

CALL SEGLD (1) 

CALL CINPUT ( IVARSU . NSLICE - I . START ) 

C 

C MAIN PROGRAM DO LOOP 



OOO OOO OOO OOO OOOO o o o 


r 

DO 200 I-l.NSLICE 

GO GET DATA CARD INFO FOR PRESENT DATA SLICE 


CALL SEGLD (1) 

CALL CINPUT (IUARSU.NSLICE, I, START) 

SPAN 'FLOAT ( ISPAN ) 

IF ( DTYPE EO - 4 ) FDGAIN-FLOAT( IFGAIN ) 

CHECK TO SEE IF CALIBRATION DATA ACQUIRED BEFORE PROCESSING FIRST DATA 
SLICE 

IF ( DTYPE EQ 1 ) FCALSU-1 
IF ( DTYPE EQ 2) ECALSU-1 
IF ( DTYPE EQ 3 AND ECALSU NE.l) UHOA-1 
IF(DTYPE EQ 4 AND FCALSU NE 1) UHOA-1 
IF(UHOA EQ. 1 ) GO TO 90 

GO GET THE DATA FROfl TAPE 

CALL SEGLD ( 2 ) 

CALL T INPUT ( DTYPE , FILNOU , FILNUN , START. NSECS ) 

DECIDE UHERE TO GO FOR EACH DATA TYPE 

GO TO (20. 10. 30. 70), DTYPE 

FILTER 60 HZ FOR ENG CAL'’ 

10 IF ( FILTSU EQ -1) GO TO 20 

CALL OUT60 (DTYPE. NSECS. FILTSU) 

CALCULATE CALIBRATION UARIABLES OF ENG OR FORCE DATA 

20 CALL SEGLD(3 ) 

CALL DCAL ( DTYPE . NSECS . IFGAIN ) 

C 



ooo ooo ooo ooo ooo ooo 


PLOT EfIG OR FORCE CAL DATA? 


25 IF ( PLOTSU EQ 1) GO TO 26 
CALL SEGLD ( 4 ) 

CALL GRAPH (DTYPE. START. NSECS. SPAN, STDERR) 

PRINT OUT HEADER? 

26 IF ( PRNTSU EQ 1 ) GO TO 200 
GO TO 90 

FILTER 60 K2 FOR EflG DATA SLICE? 

30 IF (FILTSU EQ -1) GO TO 40 

CALL OUT60 (DTYPE-NSECS, FILTSU) 

SCALE EfIG DATA. CALCULATE INTEGRATED VALUE AND MEAN SQUARE VALUE 

40 CALL SEGLD (5) 

CALL ENG ( NSECS , VOLTSEC ) 

PLOT EfIG DATA ON COMPUTE*? 

IF (PLOTSU EQ 1) GO TO 50 
CALL SEGLD ( 4 ) 

CALL GRAPH (DTYPE. START, NSECS. SPAN. STDERR) 

TAKE THE FOURIER TRANSFORM OF THE DATA AND FIND PSD UP TO 400 HZ 
50 CALL SEGLD (6) 

CALL FFTPSD (NSECS . SPAN . PSDSUM . STDERR ) 

DTYPE-5 

PLOT THE PSD ON THE COMPUTEK? 

IF ( PLOTSU EQ 1 ) GO TO 60 
CALL SEGLD( 4 ) 

CALL GPAPH (DTYPE. START, NSECS. SPAN. STDERR) 



o o o o o o o o o ooooooooooooo 


PRINT OUT RESULTS ON LINE PRINTER 7 

60 IF ( PRNTSU EQ 1 ) GO TO 200 
GO TO 90 

SCALE FORCE DATA FOR PLOTTING AND AUERAGE IT FOR PRINTING. 

70 CALL SEGLD (5) 

CALL FORCE ( DTYPE . NSECS . SPAN , FDGAIN ) 

PLOT FORCE DATA ON COMPUTE*? 

IF ( PLOTSU . EG . 1 ) GO TO 30 
CALL SEGLD ( 4 ) 

CALL GRAPH (DTYPE, START, NSECS. SPAN, STDERR ) 

PRINT OUT FORCE DATA ON LINE PRINTER? 

80 IF ( PRNTSU EG 1 ) GO TO 200 
90 CALL SEGLD(7 ) 

CALL PRINT ( DTYPE , START , NSECS , SPAN , UOLTSEC . PSDSUM . STDERR , UHOA ) 

200 CONTINUE 
END 

SUBROUTINE CINPUT ( IUARSU , NSLICE . I , START ) 

SUBROUTINE CINPUT AQUIRES THE CARD DATA AND STORES IT ON A RAD FILE 

ARGUMENTS IUARSU- ARRAY HOLDING THE UARIABLSS AND SWITCHES FOR EACH 

DATA SLICE (SEE EOUIUALENCE STATEMENTS IN RAIN 
PROGRAM) 

NSLICE- NUMBER OF DATA SLICES IN PRESENT RUN 
I - IDENTIFIES PRESENT DATA SLICE BEING PROCESSED 
START - BEGIN TIME OF PRESENT DATA SLICE 



o o o o o o 


DIMENSION IUARSU CIO) , ICARD (200 15 ) 

DIMENSION IRDIOf 10 ) . IRDF ( 6 ) 

COMMON DATA (8500 ) . IHEAD < 20 ) . CAL( 2 . 4 ) 

EQUIVALENCE ( ICARD - DATAi 1 ) ) 

DATA IRDIO ( 1 ) /'RDV 

DATA IRDF /Z4000 . ' RD ' . 1 . 6000 . 0 . 0/ 

CALL DEFINE (IRDF. IRDF) 

IF ( NSLICE EQ 0) GO TO 30 
GO TO 60 

INPUT NUMBER OF DATA SLICES AND READ DATA SLICE INFO FOR PRESENT RUN 

30 READ (105.20) NSLICE 
UR ITE (108.35) NSLICE 
20 FORMAT ( 15 ) 

35 FORMAT ( 1H1 , ' titt DATA CARDS *«*',///,' NO OF DATA SLICES- '. IS. / 
0//.5X, 'DATA SLICES',//) 

IF ( NSLICE GT 200) GO TO 100 
DO 50 J-l. NSLICE 

READ (105.40) (ICARD(J.K).K-l. 15) 

40 FORMAT ( 1513 ) 

UR ITE ( 108. 40) ( ICARD (J.K).K- 1,15) 

50 CONTINUE 

CALL QOPEN( IRDIO, ICARD. 3000. 1.0) 

CALL QURITE( IRDIO. 3000. INDIC) 

CALL QCLOSE (IRDIO, 2) 

GO TO 200 

INITIALIZE VARIABLES AND SUITCHES FOR EACH DATA SLICE 

60 CALL QOPEN(IRDIO. ICARD, 3000, 0,0) 

CALL QREAD(IRDIO. 3000. INDIC) 

CALL QCLOSE (IRDIO. 2) 

IUARSU ( 1 ) • ICARD (1.2) 

IUARSU(2)*ICARD(I, 1 ) + l 

STAPT *3600 tFLOAT ( ICARD ( 1 . 3 ) ) + 60 tFLOAT (ICARD (1 . 4 ) )+FLOAT ( ICARD( I . 

*5 ) ) 

IUARSU (3 )* ICARD ( 1.6) 



ooooooooooooooo 


IUARSU( 4 )■ ICARD ( I , 7 ) 

IUARSU (5 ) * I CARD ( I - 8 ) 

IUARSU< 6 )»ICARr(I,9) 

IUARSU(7)-ICARD(I,10) 

IUARSU ( 8 ) ■ ICARD ( 1 . 11 ) 

IHEAD( 16 )• ICARD ( 1.3) 

IHEAD ( 17 ) • ICARD ( I - 4 ) 

I HEAD ( 18 ) »ICARD ( I - 5 ) 

IHEAD ( 19 ) ■ ICARDC 1,7) 

IF ( IUARSU( 1 )-2) 90-90,200 
90 IHEAD(20)-ICARD(I,12) 

GO TO 200 

100 URITE ( 108 . 1 10 ) NSLICE 

110 FORMAT ( /// , ' ERROR READ HUMBER OF DATA SLICES TO BE'. 16-' CANN 
SOT HAUE MORE THAN 200 ' ) 

STOP 

200 RETURN 
END 

SUBROUTINE TINPUT ( DTYPE , FILNOU - FILNUM, START , NSECS ) 

SUBROUTINE TINPUT READS FORCE AND EMG DATA IN FROM TAPE 

CALLING ARGUMENTS DTYPE - DATA TYPE 

FILNOU- FILE ON TAPE CURRENTLY BEING ACCESSED 

FILNUM- FILE ON TAPE WHERE DESIRED DATA IS LOCATED 

START - TIME IN TOTALED SECONDS OF BEGINNING OF 
DESIRED DATA SLICE 

NSECS - NUMBER OF SECONDS OF DATA TO BE READ 


INTEGER FILMOUE- FILNUM, DTYPE. FILNOU 
DIMENSION INTP(10), IARRAY(535) 

COMMON DATA(8500 ) , IHEAD (20 1 . CAL (2, 4 ) 
EQUIVALENCE ( I ARRAY , DATA ( 8230 ) ) 



o o o ooooooooo 


intp( 1 )-e 

CALL QOPEN (INTP IARRAY , 535 , 2 , 0 > 

POSITION TAPE AT THE BEGINNING OF THE DESIRED DATA FILE 

IF (FILflOUE LT 0) FILflOUE-FILflOUE-1 
CALL OFSKIP (INTP, FILflOUE, INDIC) 

GO TO (300.5). INDIC/4 

IF (FILMOUE LT 0) CALL OFSKIP ( INTP , 1 , INDIC ) 

5 IF ( FILflOUE LT 0) FILf10UE-FILf!0U£ + l 
FILflOUE-FILNUfl-FILNOU 
FILNOU-FILMOU+FILflOUE 

POSITION TAPE AT THE BEGINNING OF THE DATA SLICE. 

CALL FIND (FILNUfl, START, INTP) 

DECIDE UHERE TO GO FOR EACH DATA TYPE 
GO TO (10,50,50, 10). i)TYPE 
ACQUIRE FORCE CAL OR FORCE DATA 
10 JJ-1 

DO 30 J* 1 , 23NSECS 

CALL QREAD ( INTP , 535 , INDIC ) 

IF( INDIC NE 0) GO TO 200 

DO 20 K-27,535.51 

DATA( JJ ) -FLOAT ( I ARRAY (K )/3 ) 

JJ-JJ+1 
20 CONTINUE 
30 CONTINUE 

IF ( DTYPE EG 1) GO TO 1000 
GO TO 150 

ACQUIRE EfIG CAL OR ENG DATA 

50 JJ-1 



ooooo ooo 


DO 30 J*1 ,2*NSECS+1 

CALL QPEAD ( INTP , 535 . INDIC ) 

IF ( INDIC NE 0) GO TO 200 
DO 70 K»1 ■ 10 
L*2S+K+50*(K-1 ) 

DO 60 KK-L.L+50 

IF(KK EQ L+l ) GO TO 60 

DATA ( JJ )■ FLOAT ( I ARRAY (KK )/8 ) 

JJ-JJ+1 
60 CONTINUE 
70 CONTINUE 
80 CONTINUE 

IF ( DTYPE EQ 2) GO TO 1000 
150 IBACK*-(2*NSECS+2) 

CALL QSKIPCINTP. IBACK, INDIC) 

GO TO 1000 

ERROR MESSAGES 

200 FILNUM-FILNUM-1 

URITEC 108 , 210 ) INDIC . DTYPE . FILNUM 

210 FORMAT ( /// , ' ERROR INDIC- '.14. 'DURING INPUT OF DATA TYPE'. 14 .' L 
0OCATED IN FILE' . 14) 

STOP 

300 FILNUn-FILNUM-1 

URITE( 108 . 310 ) FILNUM 

310 FORMAT!///,' FATAL ERROR EOT ENCOUNTERED DURING SEARCH FOR FILE' 

0.13) 

STOP 

1000 CONTINUE 

CALL QCLOSE (INTP.0) 

RETURN 

SUBROUTINE FIND ( FILNUM , START , INTP ) 

SUBROUTINE FIND POSITIONS THE TAPE AT THE BEGINNING OF A DATA SLICE 
CALLING ARGUMENTS FILNUM- FILE ON TAPE UHERE DATA SLICE IS LOCATED 



oooo ooo ooo 


START - TIME IN TOTALED SECONDS OF THE BEGINNING OF 
DESIRED DATA SLICE 


C 
0 
r 

C INTP - CONTROL BLOCK FOR OINOUT TAPE ROUTINES 

C 

DIMENSION INTP( 10 ) , IARRAYC535 ) 

INTEGER FILNUM 
REAL NOU 

COMMON DATA (8500). IHEAD ( 20 ) , CAL(2 . 4 ) 

EQUIUALENCE (IARRAY,DATA(8230) ) 

IPASS-0 

READ EACH RECORD AND COMPARE ITS TIME LABEL UITH THE START TIME 

10 CALL OREAD ( INTP . 535 , INDIC ) 

IF ( INDIC NE 0) GO TO 20 

NOU* 3600 *ABS(FL0AT(IARRAY(2) ) )+60 XABS (FLOAT ( I ARRAY ( 3 ) ) )+ABS(FLOA 
#T( I ARRAY ( 4 ) ) ) 

IF(NOU.NE START) GO TO 10 

ACQUIRE THE HEADER INFORMATION FOR THE PRESENT DATA SLICE 

DO 15 J*l, 15 
IHEAD( J)*IARRAY( J+5) 

15 CONTINUE 

CALL QSKIP ( INTP . -1 , INDIC ) 

GO TO 100 

DID NOT FIND THE START TIME BEFORE ENCOUNTERING AN EOF IF ONLY ONE PASS 
HAS BEEN MADE THRU THE DATA. BACK UP AND TRY AGAIN OTHERUISE . STOPD 

20 IF( IPASS EQ 1 ) GO TO 40 
CALL QFSKIP(INTP.-2, INDIC) 

IF ( INDIC EQ 8) GO TO 35 
CALL QFSKIP ( INTP , 1 , INDIC ) 

35 IPASS* 1 
GO TO 10 

40 FILNUM*FILNUM-1 



ooooooooooo 


URITE( 108,50 ) FILNUM 

50 FORMATC///, '.CANNOT FIND START TIP1E IN FILE'. 14) 

STOP 

100 CONTINUE 
RETURN 
END 

SUBROUTINE OUT60 ( DTYPE , MSECS , FIITSU ) 

SUBROUTINE OUTS0 IS A DIGIT NOTCH FILTER DEPENDING ON THE DATA CARD 
REQUEST, ONLY 60 H2 IS REMOVED OR ALL HARMONICS OF 60 H2 UP TO 360 H2 

CALLING ARGUMENTS DTYPE - TYPE OF DATA 

MSECS - NUMBER OF SECS OF DATA TO BE FILTERED 

FILTSU- SUITCH UHICH INDICATES UHETHER ONLY 60 H2 
IS TO BE REMOUED OR ALL HARMONICS OF 60 H2 


INTEGER DTYPE. FILTSU 

COMMON DATA(S500). IHEAD(20).CAL(2,4 ) 

EQUIUALENCE ( IRATE . IHEAD ( 1 1 ) ) 

LOOPS* 1 

TUOP 1*6 28318530717 
IF (FILTSU EQ 1 ) L00PS-5 
NPTS*1024*NSECS 
DO 30 J* 1 . LOOPS. 2 
A-0 0 
8*0 0 

U* (60*J )*TUOPI/FLOAT( IRATE ) 

DO 10 K • 1 . NPTS 

A*A+DATA(K)*COS(K*U) 

B*B+DATA(K)tSIN(K*U) 

10 CONTINUE 

A* AJ2 /FLOAT ( NPTS ) 

B*B*2 /FLOAT ( NPTS ) 

DO 20 K ■ 1 . NPTS 

DATA(K )*DATA(K )-A*COS(K*U )-B*SIN(K*U) 
20 CONTINUE 



ooo ooo ooooooooooo 


30 CONTINUE 
100 RETURN 
END 

SUBROUTINE DCAL( DTYPE . NSECS , IFGAIN ) 

SUBROUTINE DCAL CALCULATES THE SCALE FACTORS FOR THE FORCE AND ENG DATA 
THE CALIBRATION CONSTANT FOR INTEGRATED ENG AREA IS ALSO FOUND ENG 
THE CAL DATA IS THEN CONVERTED INTO UNITS OF POUNDS OR NICROUOLTS FOP 
PLOTTING ON THE COMPUTE* TERMINAL 

CALLING ARGUMENTS’ DTYPE - TYPE OF DATA 

NSECS - LENGTH OF CAL DATA 


DIMENSION HIGH( 10). LOU (10) -POUNDS (2600) 

INTEGER DTYPE. FL3CAL 
REAL LOUSUO 

COMMON DATA( 3500 ) . IHEAD ( 20 ) . CAL ( 2.4) 

EQUIVALENCE (FL3CAL , IHEAD ( 15 ) ) , ( IRATE , IHEAD ( 11)) 

EQUIVALENCE (MAGTUD . IHEAD ( 13 ) ), ( ICPS. IHEAD ( 14 ) ) 

EQUIVALENCE (POUNDS. DATA( 1 ) ) 

IF ( DTYPE EQ 2 ) GO TO 100 

FIND FORCE CALIBRATION VARIABLES 

CALU.l )*FL0AT(FL2CAL ) 

CAL (1,4) -FLOAT (IFGAIN) 

LOOK FOR JUMP IN DATA TO INDICATE FORCE CAL 

BEGIN-DATA( 80 ) 

DO 10 J-81.29SNSECS-80 
IF(DATA(J).GT BEGIN+500 ) GO TO 30 
10 CONTINUE 

URITE (108,20) 

20 FORMAT(///. 'COULD NOT FIND THE FORCE CAL ...CHECK DATA SLICE TIMES 

•' ) 



OOOO OOO O OO O 


STOP 


FOUND THE JUP1P, NOU COMPUTE THE AUERAGE BEFORE AND AFTER THE JUMP 
CAL UEIGHT IN COUNTS • AUEPAGE AFTER - AUERAGE BEFOPE 
FORCE SLOPE ■ COUNTS/POUND 

ZERO FORCE * AUERAGE COUNTS JUST BEFORE CAL JUNP 

30 HISUN-0 0 
LOUSUM-0 0 
DO 40 K-1.29 

HISUi1-HISUf1+DATA( J+K+60 ) 

LOUSUN-LOUSUN+DATA( J-K-60) 

40 CONTINUE 

FORSLP- (HISUN/20 -LOUSUfl/20 )/FLOAT(FLBCAL ) 

CAL ( 1 . 2 ) -FORSLP 

ZEROFOR-LOUSUfl/20 

CAL(1.3)-ZER0F0R 

SCALE FORCE CALIBRATION DATA 

DO 50 J* 1 - 29SNSECS 

POUNDS ( J )• (DATA< J )-ZEROFOR )/FORSL? 

50 CONTINUE 
GO TO 200 

FIND ENG CALIBRATION UARIABLES 

FIRST REHOUE ANY DC OFFSET IN THE CAL DATA 

100 SUN-0 0 

INDEX- 1000 tNSECS 
DO 110 J-l. INDEX 
SUf1-SUn+DATA( J ) 

110 CONTINUE 

OFFSET -SUN /FLO AT ( INDEX ) 

DO 129 J-l. INDEX 

DATA( J )-DATA( J ) -OFF SET 

DATA( 4000+ J ) -ABS( DATA( J )-OFF SET ) 

120 CONTINUE 
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C CALCULATE STANDARD DEVIATION AND OBTAIN SCALE FACTOR FOR ZERO MEAN 

STATIONARY SIGNAL STANDARD DEUIATION. RMS UALUE 

XSQUAR.0 0 
DO 140 .1-1 , INDEX 
XSQUAR*XSQUAR+DATA( J )XX2. 

140 CONTINUE 

EilGCAL ‘SORT ( XSQUAR/ ( INDEX- 1 ) )/FLOAT ( MAGTUD ) 

CAL(2. 1 )* FLOAT ( NAGTUD ) 

CAL (2,2 ) "EMGCAL 

SCALE CAL DATA FOR PLOTTING 

DO 180 J-l, INDEX 
DATA( J )-DATA( J )/EMGCAL 
180 CONTINUE 
200 RETURN 
END 

SUBROUTINE GRAPH (DTYPE , START , NSECS , SPAN , STDERR ) 

SUBROUTINE GRAPH PLOTS ON THE COMPUTE* TERMINAL THE DATA SLICE BEING 
PROCESSED ALONG UITH UARIOUS HEADER INFORMATION 

CALLING ARGUMENTS DTYPE - TYPE OF DATA 

START - BEGIN TIME OF PRESENT DATA SLICE 

NSECS - LENGTH OF DATA SLICE IN SECONDS 

SPAN - AVERAGING INTERVAL FOR PSD 

STDERR- NORMALIZED STANDARD ERROR FOR PSD ESTIMATE 

DIMENSION Z( 15 ) , DATE( 3 ). POUNDS (2600 ) , PSDN(410 ) 

DIMENSION FREQ (410 ) ,CALDATA(2000 ) 

INTEGER DTYPE 

COMMON DATA ( 8500 ). IHEAD ( 20 ) , CAL( 2 . 4 ) 
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EQUIVALENCE ( POUNDS , DATA ( 1 ) ) . ''PSDN, DATA(8000 ) ) 
EQUIVALENCE ( FREQ , DATA ( 5*00 )). (FLBCAL . CAL ( 1 . 1 ) ) 
EQUIVALENCE (CALDATA , DATA( 1 ) ). ( EMGNAX , DATA ( 8500 ) ) 
EQUIVALENCE (FORNAX . DATA< 8500 ) ) 

SUB-FLOAT(IHEAD(l ) ) 

RUN -FLOAT ( I HE ADC 9 ) ) 

DATE( 1 ) - FLOAT ( IHEAD ( 2 ) ) 

DATE ( 2 ) -FLOAT ( I HEAD ( 3 ) ) 

DATE ( 3 ) - FLOAT( IHEAD( 4 ) ) 

FLBCAL-FLOATC IHEAD ( 15 ) ) 

SHOUR-FLOATC IFIX($T ART/3600 )) 

SNIN -FLOAT ( IFIXC ( START-SHOUR*3600 . )/60 )) 
SSEC-START-SHOURS3600 -SMIN*60 
FINISH- ST ART+FLO AT (NSECS ) 

FHOUR- FLOAT (IFIX( FINISH/3600 )) 

FNIN- FLOAT ( IFIXC ( FINISH-FHOUR*3600 )/60 )) 
FSEC-FINISH-FHOUR*3600 -FMIN*S0 

INITIALIZE PLOT AND OUTPUT HEADER INFORMATION 

CALL INITAL (Z) 

CALL MODSETCZ. 4, 10 ) 

CALL «ODSET(Z,5. 10 ) 

CALL OBJECTCZ. 75 .1008 ,100. .770 ) 

CALL UORDSCZ.660 .780 .16.' • ■ •) 

CALL UORDS(Z. 100 ,30 , 4 » ' SUB* ' ) 

CALL UORDSCZ, 440 ,30, 16. 'EXP DATE / / ') 

CALL UORDS ( Z , 840, 30, 4 , ' RUN? ' ) 

CALL NUMBER (Z. 660 ,780 ,SHOUR,2.0) 

CALL NUNBERCZ. 696 ,780 ,Si1IN,2.0) 

CALL NUf1BER(Z. 732 ,780 .SSEC.2.0) 

CALL NUMBER ( Z , 780 ,780 .FHOUR. 2.0) 

CALL NUMBER ( Z . 816 ,780 .FHIN,2,0) 

CALL NUf1BER(Z,85£ ,780 .FSEC.2.0) 

CALL NUP1BER( Z. 160 ,30 .SUB, 3.0) 

CALL NUMBER* Z . 548 .30 , DATE ( 1 ) . 2 . 0 ) 

CALL NUMBER(Z.584 ,30 , DATE ( 2 ) . 2 . 0 ) 

CALL NUMBER(Z . 620 ,30 , DATE ( 3 ) . 2 . 0 ' 
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CALL NUnBER C Z . 830 .30 .RUN, 2.0) 

C 

C DEC T DE UHERE TO GO FOR EACH DATA TYPE 

C 

GO TO (10.40.40.10,70),DTYPE 

LABEL AXES FOR FORCE DATA 

10 CALL MOBSET (2,2, 1 ) 

NPTS-202NSECS 
PTS-FLOF.T(NPTS) 

CALL UORDS (2.13 . 380 . 4 , ' LBS ') 

CALL U0RDS(2.25 ,355 ,2, 'OF') 

CALL UORDS(2,0 ,330 ,6, 'FORCE ') 

CALL UORDS<2.450 ,55 ,16, 'TINE IN SECONDS ') 

IF ( DTYPE EG 4) GO TO 30 

FINISH UP PLOT FOR FORCE CALIBRATION DATA 

CALL UORDS(2. 100 .780 . 16. 'FORCE CALIBRATION') 

CALL U0RDS(7 . 350 ,0 ,26. 'AMOUNT OF FORCE CAL (LBS) •') 
CALL NUMBER <2,674.,0 .FLBCAL.3,0) 

CALL SUBJEC ( 2 . 0 . FLOATC NSECS ) . 0 . . 200 > 

CALL GRID(2. 2 .25 ) 

CALL LABEL(2. 1.2 ,2,0) 

CALL LABEL ( 2 , 2 , 50 . . 3 . 0 ) 

CALL SUBJEC (2.0 .PTS.0 .290 ) 

CALL M0DLINE(2.NPTS, 1 . POUNDS) 

CALL BELL(2) 

CALL PAGE ( 2 ) 

GO TO 90 

FINISH UP PLOT FOR FORCE DATA 

30 CALL U0RDS(2. 100 .780 .10, 'FORCE DATA') 

FMAX*200 

IF (FORMAX LT 80 ) FMAX«100 
CALL SUBJEC (Z.0 . FLOAT ( NSECS ). 0 , FMAX ) 
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CALL GRID (Z 2 .10 ) 

TICK*2 +FL0AT ( NSECS/31 ) 

CALL LABEL (Z • 1 . TICK . 3 , 0 ' 

CALL LABEL (Z.2 20 .3,0) 

CALL SUBJEC (2.0 .PTS 0 , FMAX ) 

CALL MODLINE(Z.NPTS. 1, POUNDS) 

CALL BELL ( 2 ) 

CALL PAGE ( 2 ) 

GO TO 90 

FINISH UP PLOT FOR EMG DATA 

40 CALL UORDS ( 2 . 430 ,55 ,20. 'TIME IN MILLISECONDS') 
CALL MODSET (2.2.1 . ) 

CALL UORDS ( 2 , 0 . 440 . , 6 , ' MICRO ') 

CALL UORDS (2.0 .410 ,6,'UOLTS ') 

NPTS-1000 

SECMIL -FLOAT ( 10O0XNSECS ) 

IF ( DTYPE EQ 3) GO TO 50 

EMG CAL DATA 

CALL U0RDS(2, 100 .780 ,8, 'EMG CAL ') 

CALL SUBJEC (Z.0 , SECMIL . -1000 . , 1000 ) 

CALL GRID (Z.250 . 250 . ) 

CALL LABEL ( 2 , 1,500 .4.0) 

CALL LABEL (2,2,500 .50) 

CALL MODLINE ( 2 . NPTS .MSECS . CALDATA ) 

GO TO 60 

EMG DATA 

50 CALL UORDS ( 2 . 100 .780 .8. 'EMG DATA') 
ULIM*(FLOAT(IFIX( EMGMAX/500 ))+l )*500 
BLIM--ULIM 

CALL SUBJEC (Z.0 . SECMIL , BLIM . ULIM ) 

HACK *ULIM/4 

CALL LAPEL ( 2 1.500 .4,0) 



CALL GRID ( Z . 2S0 .HACK) 

uNun-uLin/2 

CALL LABEL (Z . 2 , UNUH . 5 . 0 ) 

CALL MODLINE(Z. NPTS. MSECS. DATA) 

60 CALL BELL ( Z ) 

CALL PAGE ( Z ) 

GO TO 90 
C 

C FINISH UP PLOT FOR PSD DATA 

C 

70 CALL UORDS(Z. 100 . 780 . 22 . ' POUER SPECTRAL DENSITY') 
CALL U0RDS(Z,490 . 55 18 , 'FREQUENCY IN HERTZ') 

CALL UORDS (Z.0 .425 .6. 'MICRO ') 

CALL UORDS ( Z , 0 . . 405 . 6 . 'UOLTS ') 

CALL UORDS (Z.60 .410 ,2. '2 ') 

CALL UORDS (Z.0 0.385 ,6.'X1000 ') 

CALL UORDS ( Z . 100 . 0 ,16. 'BANDUIDTH (HZ) -') 

CALL UORDS (Z.440 ,0 , 26 . 'NORMALIZED STANDARD ERROR-') 
CALL NUMBER ( Z , 232 .0 .SPAN, 2, 2) 

CALL NUMBER (Z.824 ,0 .STDERR.1,3) 

CALL SUBJEC (Z.0 .400 .0 .1 1) 

CALL GRID (Z.25 , 05) 

CALL LABEL (Z.1.25 .3,0) 

CALL LABEL (Z.2. .10001,1.1) 

C 

C CREATE X-AXIS FOR PSD UALUES 

C 

NPTS*IFIX(400 /SPAN) 

DO SO J-l.NPTS 
FREG( J)-SPANSFLOAT( J) 

80 CONTINUE 

CALL LINES(Z.NPTS.FREQ.PSDN) 

CALL BELL(Z) 

CALL PAGE ( Z ) 

90 RETURN 
END 

SUBROUTINE MODLINE ( Z , NPTS . MODX , Y ) 

C 
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SUBROUTINE flODLINE PLOTS A Y-ARRAY OF DATA POINTS UITHOUT REQUIRING A 
CORRESPONDING X-ARRAY INSTEAD OF AN X-ARRAY, THE ARRAY POSITION OF A 
Y-ELEMENT IS USED FOR THE X UALUE 

CALLING ARGUMENTS 2 - THE PLOT ARRAY 

NPTS - THE NUMBER OF POINTS TO BE PLOTTED 

MODX - X UALUE INCREMENT (X-l.EUERY Y POINT; X-2. 
EUERY OTHER Y POINT; ETC ) 

Y - THE Y-ARRAY 

SPECIAL CONSIDERATION : NPTSSMODX SHOULD NOT BE GREATER THAN THE LENGTH 

OF THE Y-ARRAY 

DIMENSION 2(15). Y(l) 

CALL POINTS (2.1, MODX , Y ) 

J2 • IURDXY (2(2), 1) 

DO 70 1-2, NPTS 
IY* IJMODX 
XX-FLOAT (I*MODX) 

GO TO ( 40 , 50 ) , J2 
40 YY-UORDXY ( Y , IY ) 

GO TO 60 
50 YY-Y(IY) 

60 CALL UCTOR( SCALE(2 , XX , 0 ) , SCALE (2, YY , 1 ) , 1 ) 

70 CONTINUE 
RETURN 
END 

SUBROUTINE UORDS (2,X,Y,NL,L) 

DIMENSION 2 ( 15 ) . L( 50 ) , M( 50 ) 

NU-NL/2 
DO 10 I-l.NU 
10 M( I )«L( I ) 

CALL PASTOR (0.M.NU) 

CALL LEGEND (2 - X, Y , M , NL ) 

RETURN 
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END 

SUBPOUTINE BELL (2) 

CALL NODE (2.0) 

CALL OUT (7,1 ) 

CALL SCRIBE 
DO 10 1-1,300 . 

DO 10 J-1,60 
10 K-I+J 

CALL OUT (2.1) 

CALL SCRIBE 
30 CALL IN (I) 

IF ( I NE 13) GO TO 20 

RETURN 

END 

SUBROUTINE EfIG ( NSECS , UOLTSEC ) 

SUBROUTINE EfIG SCALES ENG DATA FOR PLOTTING ON THE COMPUTER AND FINDS 
THE INTEGRATED EMG VALUE FOR THE DATA SLICE . 

CALLING ARGUMENTS NSECS - LENGTH OF DATA SLICE 

VOLTSEC- INTEGRATED EMG VALUE IN MICROVOLTtSEC 

INTEGER DTVPE 

COMMON DATA(8500), IHEAD(20).CAL(2,4) 

EQUIVALENCE ( IRATE , IHEAD ( 1 1 ) ) , ( EMGCAL , CAL (2.2)) 

EQUIVALENCE ( EMGMAX . DATA( 8500 ) ) . ( EMGVAR . DATA( 3499 ) ) 

SUM-0 0 

CALCULATE OFFSET. SUBTRACT FROM DATA AND APPLY SCALE FACTOR 

DO 10 J* 1 . 1024TNSECS 
SUM»SUM+DATA( J ) 

10 CONTINUE 

OFFSET»SUM/( 1024 *FLOAT( NSECS ) ) 

CAL(2.4 J-OFFSET 
EMGMAX-0 0 

DO 20 J-1,1024*NSECS 



c 

CALI RFORT ( DATA . NEXP , S . -1 . IFERR ) 

C 

C CALCULATE RAU POUER SPECTRAL DENSITY 

C 

H-l /FLOAT ( IRATE ) 

K * 1 

DO 10 J-3.N+1.2 

PSDOC )"8 XH/PTSt ( DATA( J )**2+DATA( J+l )»2 ) 
K-K+l 

10 CONTINUE 
C 

C CORRECT PSD AMPLITUDES FOR UINDOU REDUCTION 

C 

DO 15 J-l.N/2 
PSD(J)-1 1429SPSD(J) 

15 CONTINUE 
C 

C OUTPUT RAU PSD TO TAPE 

C 

CALL PSDSAUE (NSECS) 

C 

C CALCULATE SMOOTHED PSD 

C 

F0- 1 . /(PTS/FLOAT ( IRATE ) ) 

IOTA- IFIX( SPAN/F0 ) 

SPAN "FLOAT ( IOTA )SF0 

PSDMAX-0 0 

K-l 

DO 30 J-l .N/2. IOTA 
TEMP-0 0 
DO 20 JJ-1 , IOTA 
TEf1P-TEMP+PSD( J+JJ-1 ) 

20 CONTINUE 

PSD CK )• TEMP/FLOAT ( IOTA ) 

IF ( FLOAT (J/ IOTA )TSPAN GT -400 +SPAN ) GO TO 25 
IF (PSD(K ) GT PSDMAX ) PSDMAX-PSDdO 
25 K-K + l 
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30 CONTINUE 

CALCULATE THE NORMALIZED PSD 

Lin- IFIX( 400 /SPAN ) + l 
DO -40 J-l.LIM 
PSDN< J)-PSD( J)/PSDHAX 
•4e CONTINUE 

INTEGRATE THE SMOOTHED PSD 

PSDSUN-0 0 

DO 50 J» 1 - LIM 

PSDSUM-P:>DSU?1+SPANSPSD(J) 

50 CONTINUE 

CALCULATE X EACH BANDWIDTH IS OF TOTAL AND FIND THE CUNLATIUE X OF TOTAL 
DO 60 J-l , LIM 

PERCNT ( J ) -SPAN*PSD( J )/PSDSUMX100 
IF < J NE 1 ) GO TO 55 
CUi’PCT ( 1 ) » PERCNT ( 1 ) 

GO TO 66 

55 CUnPCT(J)»CUMPCT(J-l )+PERCNT(J) 

60 CONTINUE 

CALCULATE THE EXPECTED UALUE AND UARIANCE OF THE PSD 

F-SPAN/2 
EXPUAL-0 0 
DO 70 J-l, LIM 

EXPO AL- EXPO AL+( PERCNT ( J )/100 )2F 
F-F+SPAN 
70 CONTINUE 
F-SPAN/2 
OAROAL-0 0 
DO 80 J-l, LIM 

OAROAL*OAROAL+( PERCNT ( J )/100 )T(F-EXPUAL )T22 
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F-F+SPAN 
S0 CONTINUE 

S T DEU-SQRT(UARUAL> 

FIND STANDARDIZED NORMAL ERROR 

STDERR-SQRTCl /FLOAT ( IOTA ) ) 

RETURN 

END 

SUBROUTINE UINDOU (NSECS) 

SUBROUTINE UINDOU APPLIES A COSINE TAPER TO THE FIRST AND LAST TENTHS OF 
THE DATA TO REDUCE LEAKAGE 

CALLING ARGUMENTS ; NSECS - LENGTH OF THE DATA 

COMMON DATA ( 3500 ) , I HEAD ( 20 ) . CALCS. 4 ) 

PI*3 1415927 

NPTS- 1024SNSECS 

IEDGE-NPTS/10 

TTOTAL-1 024ZFLOAT (NSECS ) 

K-l 

DO 10 J-l.IEDGE 
Tl- 0013FLOATCJ) 

Cl- S*<1 -COSCPISTl/C 1ST70TAL ) ) ) 

DATAC J)-CIJDATA(J) 

DATACNPTS-J+l )-ClSDATA(NPTS-J+l ) 

10 CONTINUE 
RETURN 
END 

SUBROUTINE RFORTC A . M . S , IFS . IFERR ) 

ONE-DIMENSIONAL REAL FINITE FOURIER TRANSFORM 


FOURIER TRANSFORM SUBROUTINE FOR REAL DATA 
PROGRAMMED IN SYSTEM/360, BASIC PROGRAMMING SUPPORT, 
FORTRAN IU. (SEE FORM C2S-6504 ) 

THIS DECK IS SET UP FOR IBSYS ON THE IBM 7094 




THIS PROGRAM USES THE SUBROUTINE FORT TO COMPUTE COMPLEX 
FOURIER TRANSFORMS OF REAL DATA PK FORT S D A NO 3465 IS 
AVAILABLE THROUGH SHARE 

THE FOURIER SERIES IS 

X(J)« SUM OUER IC-0 TO N. OF COC )*EXP(2tPI*I SJ*K/N ) 

J-0.1.2 N-l 

WHERE I-SCRT(-i) AHD WHERE C(X) IS COMPLEX 

SINCE X(J) IS REAL. C(SC)« CCNJG(C(N-K) ) THEREFORE ONLY 

C(!O.K- 0 . 1 .,.. N/3 ARE COMPUTED AND/OR USED 

ARGUHENTS- 

A IS INITIALLY THE INPUT ARRAY, X, WHEN COMPUTING A FOURIER 
TRANSFORM AND C UHEH COMPUTING A FOURIER SERIES A IS REPLACES BY 
THE OUTPUT ARRAY. C IN 7HE FORMER CASE. X IN THE LATTER 
THE X VECTOR CONTAINS THE REAL DATA X(0 ) , X< 1 ) . , X( N- 1 ) 

THE C VECTOR CONTAINS THE COMPLEX FOURIER AMPLITUDES 

C(0),C(1) C ( ft/2 ) . THE COMPLEX VECTOR C IS STORED ACCORDING 

TO THE NORMAL FORTRAN IV CONVENTION FOR STORING COMPLEX NUMBERS 
IE. REAL PARTS IN ALTERNATE CELLS STARTING WITH THE FIRST. 
IMAGINARY PARTS IN ALTERNATE CELLS STARTING UITH THE SECOND 

TO ADHERE TO FORTRAN RULES. X(0>. X(l) ARE REFERRED TO AS 

X(i).X(2) RES? . IN THE PROGRAMS ALSO. CO).C(l), . ARE 

REFERRED TO AS C(1),C(2> RES? . IF C IS DESIGNATED AS 

COMPLEX IN A TYPE STATEMENT 

M GIVES N-2S3M 

THE ARGUMENTS S. IFS. AND IFERR ARE THE SAME AS IN THE 
SUBROUTINE FORT AND THE USER IS REFERRED TO THE COMMENT CARDS 
IN FORT THEIR EXPLANATION 

DIMENSION STATEMENTS- THE DIMENSIONS OF ARRAYS A AND S SHOULD 
BE N+2 AND N/4 RESP FOR THE LARGEST N TO BE USED FOR 
EXAMPLE. IF THE LARGEST M IS 13. THEN. N»8192 AND ONE SHOULD 


ORIGINAL PAGE IS 
OF POOR QUALITY; 
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HAUE THE DIMENSION STATEMENT- 
DIMENSION A( 8194 ) , S(2048) 

IF ONE WISHES TO SPECIFY A TO BE COMPLEX BY A TYPE STATEMENT. 
ONE SHOULD GIUE IT A DIMENSION OF N/2 +1 , FOR THE LARGEST N 
DIMENSION A(8500 ) . S( 1024 ) 

10 IFERRS • 9 
N-22SN 
NU2 -N/2 
NU4R1 - N/4 - 1 
P1P11 • M - 1 

IF < ABS( IF3 )-l ) 40,40.20 
20 IF (MP-M) 30,50.50 
30 IFERRS • 1 
40 NP • N 
HP - M 

CALL FCST( A. M , S , 0 . IFERR1 ) 

IFERRS • IFERRS + IFERR1 
50 (CD • NP / N 
KT • KD 
NPU4 - NP / 4 
IF (IPS) S0.80.90 
60 CALL FCRT( A. MM1 , S. -2 , IFERR2 ) 

IFERRS - IFERRS + IFERR2 
DO 70 K-1.NU4M1 
J-NU2-K 

A1R» A(22X+1) + A(2SJ+1) 

All «A(22X+2 )-A(23J+2 ) 

A2R"A(22K+2)+A(2$J+2) 

A2I »A(25J+1 )-A(23K+l ) 

KK7 » NPU4-XT 
AUR*A2R2S(KKT ) + A2I5SOCT) 

AUI • A2ISSOCKT )- A2R*S(KT> 

A(22K + 1 )• ( A1R+AUR )/4 
A(22K+2)*(A1I+AUI )/4 
A( 23J+1 ) • ( A1R-AUR )/4 
A(2tJ + 2)"(AUI-AlI )/4 
70 KT-KT+KD 
T * A ( 1 ) 
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30 


A(1 )*(T+A(2) J/2 
A < N+ 1 ) * (T-A*:2 ) )/2 
A ( 2 ) *0 
A ' N + 2 ) » 0 

A( NU2 + 1 ) ■ 5*A(NU2+1) 

A f NU2 + 2 ) ■ - S * A ( NU2+2 ) 

IFERP • IFEPRS 
RETURN 
90 DO 100 K-1.NU4M1 
J-NU2-K 

A1R-A(2r01 ) + A(2*J + 1 ) 

AlI»A(2*K+2)-A(2*J+2) 

AUR*A ( 2TK + 1 )-A(2*J+l ) 

AUI *A ( 2*K+2 )+A(2*J*2 ) 

KKT • NPU4 - KT 
A2R-AUR*S(KKT) - AUI*S(KT) 

A2I*AUR*S(KT) + AUItS i KKT ) 

A ( 2*K+ 1 ) ■ AIR - A2I 
A ( 2*K+2 ) • All + A2R 
A( 2TJ + 1 ) ■ AIR + A2I 
A(2*J+2) ■ A2R - All 
100 KT • KT + KD 
T • A ( 1 ) 

A ( 1 ) » T + A ( N+ 1 ) 

A(2 ) ■ T - A ( N+ 1 ) 

A ( NU2+ 1 ) • 2 tA(NU2*l ) 

A(NU2+2) *-2 *A(NU2+2) 

CALL c’ORT ( A , MM1 . S . 2 , I FERRE ) 

IFERRS • IFERRS+IFERR2 

GO TO 80 

END 

SUBROUTINE FORT ( A . M . S , IFS . IFERR ) 

FORT. ONE-DINENSIONAL FINITE COnPLEX FOURIER TRANSFORM 


FOURIER TRANSFORM SUBROUTINE- PROGRAMMED IN SYSTEM/360. 
BASIC PROGRAMMING SUPPORT. FORTRAN IU FORM C28-6504 
THIS DECK SE T UP FOP IBSYS ON IBM 7094 



POES EITHER FOURIER SYNTHESIS I E .CONFUTES COMPLEX FOURIER SEPI r 
GI'JEN A UECTOR OF N COMPLEX FOURIEP AMPLITUDES , OR . GIUEN ft UECTOP 
OF COMPLEX DATA X DOES FOURIER ANALYSIS. COMPUTING AMPLITUDES 
A IS A COMPLEX UECTOR OF LENGTH N-2**M COMPLEX NOS. OR 2*N PEAL 
C NUMBERS A IS TO BE SET BY USER 

C M IS AN INTEGER 0 LT M LE 13. SET BY USER 

C S IS A UECTOR S(J)- S IN ( 2*P I * J/NP ). J-1.2. . .NP/4-1. 

C COMPUTED BY PROGRAM 

C IFS IS A PARAMETER TO BE SET BY USER AS FOLLOUS- 

C IFS-0 TO SET NP*22~M AND SET UP SINE TABLE 

C IFS* 1 TO SET N-NP*2*tM, SET UP SIN TABLE. AND DO FOURIER 

C SYNTHESIS. REPLACING THE UECTOR A BY 
C 

C XCJ)* SUM OUER K-0.N-1 OF A(K ) *EXP(2*PI*I/N )**( J*K ) . 

C J*0 . N-l , UHERE I *SQRT ( -1 ) 

C 

C THE X'S ARE STORED UITH RE X(J) IN CELL 2SJ+1 

C AND IM X(J) IN CELL 2XJ+2 FOR J-0 , 1 , 2 . . , N-l 

C THE A ' S ARE STORED IN THE SAME MANNER 

C 

C IFS--1 TO SET N*NP*2SJM. SET UP SIN TABLE. AND DO FOURIER 

C ANALYSIS. TAKING THE INPUT UECTOR A AS X AND 

C REPLACING IT BY THE A SATISFYING THE ABOUE FOURIER SERIES 

C IFS>0 SET UP SIN TABLE AND RETURN 

C IFS-+2 TO DO FOURIER SYNTHESIS ONLY, UITH A PRE-COMPUTED S 

C IFS*-2 TO DO FOURIER ANALYSIS ONLY, UITH A PRE-COMPUTED S 

C 

C IFERR IS SET BY PROGRAM TO- 

C *0 IF NO ERROR DETECTED 

C *1 IF M IS OUT OF RANGE , OR. UHEN IFS-+2.-2. THE 

C PRE-COr°UTED S TABLE IS NOT LARGE ENOUGH 

C *-l UHEN IFS *+1,-1, MEANS ONE IS RECOMPUTING S TABLE 

C UNNECESSARILY 

C 

C NOTE- AS STATED ABOUE. THE MAXIMUM UALUE OF M FOR THIS PROGRAM 

C ON THE IBM 7094 IS 13 FOR 360 MACHINES HAUING GREATER STORAGE 

C CAPACITY. ONE MAY INCREASE THIS LIMIT BY Rrot.acrv^ -• 



0 STATEMENT 3 BELOU BY L0G2 N. UHEPE n 15 THE MAX NO OF 
•: COMPLEX NUMBERS ONE CAN STORE IN HIGH-SPEED COPE ONE PIUS'" 

C ALSO ADD MORE DO STATEMENTS TO THE BINARY SORT ROUTINE 

FOLLOUING STATEMENT 24 AND CHANGE THE EQUIUALENCE STATEMENTS 
C FOR THE K ' 5 

r 

DIMENSION A(8500).S(1024),K(15) 

IF (M) 20.20,10 
10 IF (N-1S) 40,40.20 
20 IFERR* 1 
30 RETURN 
40 IFERR-0 
N • 2tJM 

IF (IABS(IFS)-l ) 440,440.50 

C UE ARE DOING TRANSFORM ONLY SEE IF PRE-COMPUTED 
C S TABLE IS SUFFICIENTLY LARGE 
50 IF (N-NP) 70.70.60 
60 IFERR* 1 
GO TO 4**0 

C SCRAMBLE A. BY SANDE'S METHOD 

70 K( 1 )*2*N 
DO 80 L*2. M 
80 K(L)*K(L-1 )/2 
DO 90 L*M , 14 
90 < ( L + l ) *2 

C THE FOLLOUING 15 STATEMENTS ARE TO COMPENSATE FOR A WEAKNESS IN 

C THE FORTRAN U COMPILER 

K 1 *K ( 15 ) 

K2 * K ( 1 4 ) 

K3 - K C 1 3 ) 

K4 *K ( 12 ) 

K5 *K( 1 1 ) 

K6 *K ( 10 ) 

K7 *K( 9 ) 

K8 *K ( 3 ) 

KO *K( 7 ) 

K10*K(6> 

K 1 1 * K ( 5 ) 


K12*K( 4 ) 

K13*K(3 ) 

K14*K(2 ) 

K1 5 *K ( 1 ) 

N2 *K ( 1 ) 

C NOTE EQUIVALENCE OF KL AND K ( i 4-L ) 
C BINARY SORT- 



IJ 

*2 

100 

JI 

•2 

110 

J 2 

•JI 

120 

J3 

■ J2 

130 

J 4 

• J3 

140 

J5 

• J4 

150 

J6 

• J5 

160 

J7 

• J6 

170 

J8 

• J7 

180 

J9 

• J8 

190 

J10 

• J9 

200 

Jll 

•J10 

210 

J12 

•Jll 

220 

J13 

• J12 

230 

J 1 4 

• J13 


240 JI-J14 

2S0 IF (IJ-JI) 260.270,270 
260 T • A ( I J - 1 ) 

A( IJ-1 )*A( JI-1 ) 

A ( J I - 1 ) ■ T 
T-A(IJ) 

A( I J )«A( JI ) 

A( JI )-T 
270 IJ-IJ+2 
JI * JI+K14 

IF (JI LE K15 ) GO TO 250 
J 1 4 • J 1 4+K 13 

IF ( J 1 4 LE K14) GO TO 240 
J13*J13+<12 

IF ( J13 LE K13) GO TO 230 
J12* J12+K1 1 




IF ( J12 LE K12) GO TO 220 
J11-J11+K10 

IF (Jll LE.K11) GO TO 210 
J 10* J10+K9 


IF 

(J10.LE.Kli 

3.) 

GO 

TO 200 

J9 

* J9+K8 





IF 

C J9 LE 

K9 ) 

GO 

TO 

190 

J8' 

■ J8+K7 





IF 

(JS LE 

K8 ) 

GO 

TO 

180 

J7’ 

• J7+K6 





IF 

(J7 LE 

K7 ) 

GO 

TO 

170 

J6* 

■J6+KS 





IF 

(J6.LE 

K6 ) 

GO 

TO 

160 

JS- 

> J5+K4 





IF 

(JS LE 

KS ) 

GO 

TO 

150 

J4> 

■J4+K3 





IF 

( J4 LE 

K4) 

GO 

TO 

140 

J3- 

■ J3+K2 





IF 

(J3.LE 

K3 ) 

GO 

TO 

130 

J2- 

■ J2+K1 





IF 

( J2. LE 

K2) 

GO 

TO 

120 

Jl- 

1 Jl+2 





IF 

(J1 LE 

K1 ) 

GO 

TO 

110 


IF ( IFS ) 2S0 , 20 . 300 

C DOING FOURIER ANALYSIS, SO DIO BY N AND CONJUGATE 
2S0 FN • N 

DO 290 I-l.N 
A(25I-1 )*A(22I-1 ) 

290 A(2SI )*-A(25I ) 

C SPECIAL CASE- L*1 
300 DO 310 I-l.N, 2 
T • A ( 25 1 - 1 ) 

A(25I-1 ) -T + A ( 251 + 1 ) 

A(25l + 1 )-T-A(25I + l ) 

T * A ( 25 1 ) 

A ( 25 1 ) • T + A(25I+2) 

310 A^ 251+2 )■ T - A(25I+2) 

IF (H-1 ) 20,30.320 



C SET FOP 1*2 
320 IEXP1-2 
C LEXPl-a«(L-l) 

LEXP *8 

C LEXP ■ 2** < L + 1 ) 

MPl * 2«nT 

C NPL • HP* 2**-L 
330 DO 390 L-2.J1 
C SPECIAL CASE- J*0 
DO 340 1-2. N2. LEXP 
II - 1 + LEXP1 
12*11+ LEXP1 
13 * I2+LEXP1 
T* A ( I - 1 ) 

A( 1-1 ) * T +A( 12-1 ) 

AC 12-1 ) ■ T-A( 12-1 ) 

T * A ( I ) 

A( I ) * T+A ( 12 ) 

A( 12 ) * T-A ( 12 ) 

T* -A( 13 ) 

TI • A( 13-1 ) 

AC 13-1 ) • AC 1 1-1 ) - T 
A( 13 ) » A ( 1 1 ) - TI 

AC 1 1-1 ) • AC 1 1 - 1 ) +T 
340 A( II ) • A( II ) +TI 

IF (L-2 ) 380.380.3S0 
350 KLAST *N2-LEXP 
JJ-MPL 

DO 370 J*4, LEXF1 . 2 
NPJJ-MT-JJ 
UR*S(NPJJ ) 

UI*S(JJ) 

ILAST * J+<LAST 

DO 3S0 I-J. ILAST. LEXP 

11*1 +LEXP 1 

12*1 1+LEXP1 

13* I2+LEXP1 

T«A(I2-1)*UR-A(I2)*UI 
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IhRRAV( J1*0 
10 CONTINUE 
DO 20 J* 1 • 20 
IARRAV( J '• IHEAD' J ^ 

20 CONTINUE 

NPECS*NSECS/2 ' 

IRNRECS EQ 0) NRECS-1 

I ARRAY ( 2 1 )-NRECS 

CALL QURITE( IOUTP. 2S, INDIC) 

OUTPUT F0LL0UING DATA RECORDS CONTAINING THE RAU PSD UALUES 
K *0 

DO 50 J ■ 1 , NRECS 
DO 40 JJ»1 , 1024 
DATA ( 4099 + JJ )»DATA(K+JJ ) 

40 CONTINUE 

CALL OURITE ( IOUTP . 2043 . INDIC ) 

K-K+1024 
50 CONTINUE 

URITE TUO EOFS AND BACK UP TO JUST BEFORE THEN 

CALL QUEOF ( IOUTP , 1 , INDIC ) 

CALL QCLOSE ( IOUTP, 0) 

RETURN 

END 

SUBROUTINE PRINT ( DTYPE . START . NSECS , SPAN, UOLTSEC , PSDSUN , S7DERR . UHO 
* A ) 

SUBROUTINE PRINT OUTPUTS TO THE LINE PRINTER HEADER AND CALIBRATION DATA. 
AUERAGED FORCE DATA, AND SMOOTHED PSD SPEC-RUM 

CALLING ARGUMENTS DTYPE - TYPE OF DATA 

START - BEGIN TIME OF PRESENT DATA SLICE 

NSECS - LENGTH OF DATA SLICE 


r 

C 

r 

r 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SPAN - FOP FORCE DATA NUMBER OF SECS FOPCE DATA I 
AUERAGED OUER, FOP PSD DATA FREQUENCY 
INTEPUAL FOR SMOOTHING ESTIMh t E? 

OOLTSEC- INTEGRATED EMG OALUE (TIME DOMAIN, 

PSDSUM- INTEGRATED EMG OALUE (FPEOUENCV DOMAIN) 

STDERR- NORMALIZED STANDARD EPROP OF ESTIMATE FOP 
PSD UALUES 

UHOA - FLAG SET UHEN DATA PROCESSING ATTEMPTED 
UITHOUT CALIBRATION 


INTEGER DTYPE , UHOA 

DIMENSION I ED ATE ( 3 ) . IDDATE ( 3 ) , MUSCLE( ?. 4 ) , IFOPRC 2 ) 

DIMENSION AOGL3S(130).PSD(410),ITIME(3),PSDN(410),PERCNT(410) 
DIMENSION CUMPCT (410 ) 

COMMON DATA ( 8500 ) . IHEAD ( 20 ) , CAL ( 2 . 4 ) 

EQUIOALENCE ( ISU3N0. IHEAD ( 1 ) ) . ( IEDATE. IHEADO ) ) 

(IDDATE, IHEAD<5)), (IFLITE. IHEAD(8) ) 

( IRUN, IHEADO) ). (MUS, IHEAD ( 29 ) ) 

( I AT APE , IHEAD( 10 ) ) , ( ISAMPE, IHEAD( 11 )) 

(ISAMPH. IHEAD ( 12) ), ( AOGL3S, DATA ( 5909 ) ) 

(PSD, DATA ( 1 ) ), (PSDN, DATA (3900 ) ) 

( PERCNT,DATA( 7000 )), (CUMPCT, DATA( 7500) ) 

( EXPOAL. DATA ( 6500 )). (STDEU.DATA( 6501 ) ) 

(EMGOAR, DATA (8499 ) ) 

/'BR'.'AC'.'HI'. 'AL'. ' B' , 'IC'. 'EP', ' B',' 

' G' , 'AS' , 'TR' , 'OC' . 'NE' . 'MI' , 'US' , ' S 

'/ 


EQUIOALENCE 
EQUIOALENCE 
EQUIOALENCE 
EQUIOALENCE 
EQUIOALENCE 
EQUIOALENCE 
EQUIOALENCE 
EQUIOALENCE 
DATA MUSCLE 
, 'AL' . 'IS' 
e . 'S ' ' , 


, 'RA' , 'DI 
'OL' . 'EU' 


DATA IFORR /' F ' , 'R+'/ 

IF ( UHOA EQ 1) GO TO 150 
ITIME ( 1 )*IFIX(ST ART/3600 ) 

ITIME(2)*IFIX((START-FLOAT(ITIME(1))*3600 )/60. ) 
ITIME f 3)*IFIX(START-FLOAT(ITIME(l))t3600 -FLOAT ( ITIME ( 2 ) )t6C ) 
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BRANCH TO PROPER PRINT OUT SECTION FOP THE DATA TYPE 
IF ( DTYPE-4 ) 5.50, 100 
PRINT OUT HEADER 


S IP-1 

IF (IFLITE GE 0) IP-2 
URITE( 188 . 10 ) 

10 FORMAT( 1H1 ,//, 4SX. 'CARDIOVASCULAR LABORATORY' . // . 47X. "ENG DATA PRO 

•CESSING PROGRAM ',/,47X. ' - ------- ------ ) 

*URITE( 108.20 )ISU8N0. IEDATE, IFORR( IP ). IFLITE . IRUN . IATAPE . IDDATE . (HU 
eSCLE < J , MUS ) . J- 1 . 7 ) , I SAMPE . ISARPM 

20 FORMAT( //// , 47X , '*** HEADER INFORMATION*** ’,///. ' SUBJECT NO 7 
ex, 13. 18X. 'EXPERIMENT DATE ' . 6X, 12. '/' , 12, '/ ' , 12 . 10X. 'FLIGHT REFERE 
•NCE DAY '.A3. 13.//. ' RUN NO ■'. 11X. 12, 19X. 'ANALOG TAPE NO 7X , I 
#3. 14X, 'DIGITIZING DATE ' . 6X , 12, ' / ' . 12 . ' / ' . 12 . // . ' MUSCLE ', 9X.7A2 
e.//.' EMG SAMPLE RATE ( SAMP/SEC ) '. 15 , 8X . 'FORCE SAMPLE RATE (SAMP/ 
*SEC)‘ '.15) 

URITE(10S.30)(CAL(1. J). J-1,4) 

30 FORMAT (////, 49X , '** CALIBRATION DATA **',///.' FORCE CAL'./.' 

% CAL UEIGHT <L3S)='.F7 1 . 16X, 'COUNTS/IB ', 2X , FS 0,5X.'A 

eUG BASELINE COUNT '.2X.F6 O.SX.'CAL GAIN '.2X.F3 0) 

URITE( 108.40) (CAL(2, J). J-1.2) 

40 FORMAT ( /// . ' EMG CAL',/.' './/,' RMS AMPLITUDE (MICROVOLTS) 

0 ' , F6 0 , 7X , ' COUNTS/MICROVOLT • ' ,F5 2.//) 

GO TO 200 

PRINT OUT FORCE DATA 

50 WRITE ( 108 . 60 ) ITIME.NSECS.SPAN 

60 FORMAT (1H1.51X. '*** FORCE DATA ***' ///5X , 'START TIME-'. 13.' '.12.' 

0 ' , 12, 15X, 'DATA LENGTH (SECS )-'. 14 , 15X . 'AVERAGE FORCE INTERVAL (SE 
«CS)-' , F3 0,//) 

LOOPS -NSECS/IFIX( SPAN) 

IF< LOOPS GT 49) GO TO 66 
URITE( 103 . 62 ) 
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62 FORMAT ( 47X, ' INTERVAL ' , 5X . 'AUER AGE FORCE (LBS ) '/, 47X. ' '.c ; . 

» ' ',//) 

DO 65 J-l. LOOPS 

UR1TE( 108,64 ) J , AUGLBS ( J ) 

64 F0RMAT(49X. 13, 15X,F5 1 ) 

65 CONTINUE 
GO TO 200 

PRINT OUT PSD DATA 

66 IF( (LOOPS/2 )*2 NE LOOPS) LOOPS-LOOPS-1 
URITE ( 108 ,63) 

68 FORMAT ( 15X , ' INTERUAL ' . 5X , ' AUER AGE FORCE ( LBS ) ' , 20X , ' INTERUAL ' , 5X , ' 

GAUERAGE FORCE (LBS)' ,/, 15X, ' ' ,5X. ' ' ,20x 

9 ' ' , 5X , ' ',//) 

DO 75 J-l, LOOPS/2 
JJ-J+L00PS/2 

URITE ( 103 . 73 ) J , AUGLBS ( J ) . J J , AUGLBS ( JJ ) 

73 FORMAT (17X.I3.15X.F5 1 , 29X , 13 . 1SX , FS 1 ) 

75 CONTINUE 
GO TO 200 

100 URITE( 108. 110) ITIME.NSECS.UOLTSZC.SPAN, STDERR , PSDSUM . EXPUAL , STDEU 
0 , EMGUAR 

110 FORMAT ( 1H1 . 39X . 'ttt POWER SPECTRAL DENSITY OF EMG DATA *SS'/// 4X, ' 
GSTART TIME-' . 13, ' ■ ' , 12. ' • ' . 12, 11X, 'DATA LENGTH (SECS )- ' . 14 . 14X . ' IN 
0TEGRATED EMG (MICR0U0LT3SEC >-' .E12 4//4X, 'BANDWIDTH (HZ)-',F6 3.10 
GX. 'NORMALIZED STANDARD ERROR-', FS 3 - 5X, ' INTEGRATED PSD (MICR0U0L7S 
Q*2)-' ,E13 4.//.4X. 'MEAN (HZ )-', FS 1 . 15X. 'STANDARD DEUIATION 

G( HZ )- ' . F6 l.SX.'EMG UARIANCE ( MICR0U0LTS3:2 )- ' , E15 4,/) 

LOOPS- IF IX (400 /SPAN ) + l 
FREQ-SPAN/2 

IF ( LOOPS GT 46) GO TO 130 
URITE( 108. 115) 

115 FORMAT ( 30X , ' FREQ' . 13X, 'PSD' . 11X. 'PSD' , 10X, '% OF',10X.'CUM V/.30X 
O' (HZ)' ,SX. ' (MMU**2/HZ)' ,8X. 'NORM' ,9X, 'TOTAL' ,9X. 'TOTAL' ./, 31X. ' 

0 ' ,8X, ' ' ,8X. ' ' ,9X, ' ' ,9X, ' ' ,/) 

DO 125 J-l, LOOPS 

URITE' 108.120) FREQ ,PSD(J).PSDN(J). PERCNT ( J ) . CUMPCT ( J ) 
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120 FORMAT (30X.FS 2-7X.G10 5 SX.F6 4.7X.F6 2.8X.F6 2) 

FREQ=FREQ+SPAN 

125 CONTINUE 
GO TO 200 

130 IF (LOOPS/2 )*2 NE LOOPS) LOOPS‘LOOPS-1 
PPEQ2 ‘FLOAT ( LOOPS/2 )*SPAN+FREQ 
UPITE ( 108. 135) 

135 FORMAT ( ' FREQ' . 8X.'PSD'. 9X . ' PSD ' . 8X , 'X 0F',6X,'CUM X' , 22X, 'FREQ 
§'. 9X , ' PSD ' . 8X. 'PSD' , 8X, 'X 0F',6X,'CUM X ' /cX , ' ( HZ ) ' , 4X , ' ( MMU3X2/ 
OHS ) ' . 5X , 'NORM' ,7X, 'TOTAL' ,SX, ' TOTAL ',22X, '(HZ)' ,5X, ' ( MMU3X2/HZ ) ' - 5 

OX, 'NORM' ,6X, ' TOTAL '.6X. 'TOTAL' ' ' , 4X. ' ',SX, ' — 

§-' ,?X, ' ' ,5X. ' ' ,22X, ' ' ,5X, ' ' ,SX, ' ' ,6X, 

A' ' , 6X , ' ',/) 

DO 145 J-l, LOOPS/2 
JJ.J+L00PS/2 

URITE(108. 140) FREQ, PSD ( J ) , PSDN( J ) , PERCNT( J ) , CUMPCT( J ) , FRE02 , PSD( J 
0 J ) , PSDN ( J J ) , PERCNT ( J J ) , CUMPCTC JJ ) 

140 FORMAT ( F7 2 , 3X , G18 5 , 5X , F6 4 , SX , FS 2 , 5X , FS 2 - 20X , F6 2 - 4X , G1 0 5,5X, 
0-6 4.5X.F6 2,5X,F6 2) 

FREQ»FREQ+SPAN 

FREQ2*FREQ2+SPAN 

145 CONTINUE 
GO TO 200 

ERROR MESSAGE 

150 URITE< 108, 160) 

160 FORMAT ( ' ERROR DATA ANALYSIS ATTEMPTED WITHOUT CALIBRAT\-NN MsST 
0 CHECK DATA CARD ORDER ' ) 

STOP 

200 RETURN 
END 
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(,***** *****«****************,.********** *•>.********************************** 

C 

C 

c program: t -1 GAN 

c 

AUTHOR* WILLIAM f.% HURST A 

PuRPOSt! ru PROVIDE FIRST ORDER REDUCTION OF F.MG AND FORCE DATA. 


EXTERNAL 0 V 1 , 0 V 2 , 0 V 3 , OV 4 , 0 VS , 0 V6 , 0 V 7 
INTEGER DTYPE,FILNUM,FILTS*,PLOrSv< 
lNTtGER FCAUSn , EC ALSft / ftHOA , PKNTSft , OVL 
DIMENSION ITFILECa) 

COMMON DATUMS (4300) , 1 HEAD (20) , CAL (2, 4) , IVARSaC 10) 

COM^ON/IGO/ t!U3CLE(7#4) , I FORK ( 2 ) 

EQUIVALENCE ( 0 T YPE , I V ARSV\ Cl ) ) ) (F ILNUM , I V ARSft ( 2 ) ) 

EWUI VALENCE (FILTSft, IVARSft(3) ) # (NSECS , I V ARSft (« ) ) 

EQUIVALENCE (PLOTS*, I VARS*. (5) ) , (PRNTSft , I V AHS.v (6) ) 

EQUIVALENCE (I SPAN, IVARS.NC7)) , ( I F G A I N , I V A R S ft ( 8 ) ) 

MT0 = 2 
UVL = 0 
FCAL3N=0 
ECALSrt=0 
ft H 0 A = 0 
NSL1CE=0 

INITIALIZE MAG TAPE UNITS AND OPEN PSD OUTPUT TAPE FILE 

call rwir ("MT0 M ,0,IER) 

CALL X iM X T ("MU", 9, IE R) 

CALL MTUPD (Mro# ,, MTls0 H #0#IER) 

OPEN OVERLAY FILE 

CALL OVOPN (OVL, M EMGAN.OL M , IER) 

* 

HEAD IN OATA CAR03 

CALL OVLOD ( 0 VL . OV 1 / 0 , IER) 

CALL CINPUT CNSLICE/I* START) 

DO 200 I = 1 » NSL ICE 

GO GET OATA CARO INFO FOR PRESENT OATA SLICE 

CALL OVLOO (OVL,OV1,0,IER) 

CALL CINPUT (NSLICE/ I > START ) 

SPAN=FLOAT (ISPAN) 

IF(oTTPE.tO.a) FDGAINsrLOAT(IFGAIN) 

CHECK TO SEE IF CALIBRATION OATA ACQUIRED BEFORE PROCESSING FIRST 0 A 1 
SLICE 

IF (UTYPE.EQ. 1 ) FCALS/,’ = 1 
1FU/TVPE.EQ.2) ECALSi\ = 1 
irar»p;..EG.s.Aj;.ECALSi > . .ne.d v-hoa=i 

IT (UTTPE.EO.a.AiVO.FCALSW.NE.l) WM0A=1 
IF (NMClA.EO, 1 ) GO TO 90 

GO GET THE DATA FROM TAPE 


4 - 


CALL 0 >/LOU (UVL,HV2,0, IEP) 

CALL TInPUT(DT YPE, ITF ILE , F ILNUM# S T AR T , NStCS, 1 ) 

C 

C DEC IOC «iHERE TO GO FOR EACH OATA TYPE. 

C 

GO TO (20, 10, 30, 70) ,OTYPE 
C 

C FILTER 60 HZ FOR EMG CAL? 

C 

10 IF (FILT3rt.EQ.-I ) GO TO 20 

CALL OUT 60 (0 f YPE , NSECS , F ILTjW ) 

C 

C CALCULATE CALIBRATION VARIABLES OF EMG OR FORCE OATA. 

C 

20 CALL UVLOO (0 VL # QV3, 0 # IER ) 

CALL DC aL (0TYPE,N3ECS, IFGAIN) 

C 

C PLOT EMG OR FORCE CAL OATA? 

C 

IF (PLOTSK.EQ. 1) GO TO 25 
CALL OVLUO (OVL,UV4,0,IER) 

CALL GRAPH ( OT YPE , S I AH T , NSECS , SPAN, STOERR ) 

C 

C PRINT OUI HEADER i 

C 

25 IF(PRNTSrt.E').l) GO TO 200 
GO TO 90 
C 

C FILTER 60 HZ FOR EMG DATA SLICE? 

C 

30 IF(FlLTSW.EO.-l) GO TO AO 

CALL CTUT60 (OTYP'c NSECS, FILTSR) 

c 

C SCALE EMG DATA, CALCULATE INTEGRATED VALUE ANO MEAN SQUARE VALUE 

C 

00 CALL OVLOO (OVL , 075 , 0 , IER ) 

• CALL EMG (NSECS, VOLTSEC) 

C 

C PLOT EMG DATA ? 

C 

IF ( PLOT 31V .EQ. 1 ) GO TO 50 
CALL OVLOD CO VL , 0V«, 0, IER ) 

CALL GRAPH CUT YPE , S TAR T , NSECS , SPAN , S TDERR ) 

C 

C TAKE THE FOURIER TRANSFORM UF THE DATA AND FIND PSD Ur TO 400 HZ. 

C 

50 CALL OVLOD (U'/L , 0V6 , 0 , I ER ) 

CALL FFIPSD (NSECS, SPAN, PSOSUM, STOERR) 

D T Y P E = 5 
C 

C PLOT THE PSO ? 

C 

IF(PLOTSn.EQ.I) GO TO 60 
CALL OVLOD (OVL# 0V4 » 0 , IER ) 

CALL GRAPH (DT YPE, ST ART # NSECS # SPAN , STOERR ) 

C 

C PRINT OUT RESULTS ON LINE PRINTER? 

C 

60 IFCPRNTSW.EQ.l) GO TO 200 
GO TO 90 
C 

C SCALE FORCE DATA FOR PLOTTING AND AVERAGE IT FOR PHINTING. 

C 



I . i U- I * * • 

70 CALI. OVUM (UVLrOVS^O, IFR) 

CALL FORCE ( D W PE , NSECS , SPAN , FOG A.' ) 

C 

C PLOT FORCE DATA V 

c 

If- (PLUT&h .EO . I ) GO I 0 PO 
CALL 0 V LUO tUVL/UV'l, 0/ ILH) 

CALL GRAPH (0rYPE,SfAHr,NSLC3,GPAN,$TDERR) 

C 

C PRINT OUT FORCE DATA ON LINE PRINTER? 

C 

ttO IF (PRNTSu.EQ. 1 ) GO TO 200 
PO CALL 0 V LOO (OVL/GV7* 0# IER) 

CALL POUT COT YPE , START , NSEC S, SPAN, VOL TSEC , PSD3UM, STDEWR , *HO A ) 
200 COMIINUE 

CALL RLSK ("MTO'MER) 

CAl.L HLSL C M MT 1 " , IER) 

CALL CLOSE (OVL,IER) 

END 


ooo ooo nrsonnnnr. n 


OVERLAY UVl 

SuDROUUNr. CINPUT (NSLICE, I ,3TAR| ) 

SUBROUTINE. CINPUT A0UIRE3 THE CAkD OATA AND STORES IT 00 A RAO FILE 
ARGUMENTS! ■; SLICE- NUMBER O' OATA SLICES 

I - lOtNT IKIES PRESENT DATA SLICE BEING PROCESSED 

START - BEGIN TIME OF PRESENT DATA SLICE 
DIMENSION 1CARCH200, IS) 

Common datumsUsOO) , IHEAD(20) ,CaL(2,«),ivars«(10) 

COMMON/ 1 GO/ MUSCLE (7,4) , IF0RRC2) 

EGO I VALENCE ( ICAkO, DATUMS ( 1 ) ) 

I i<il EGER COR,OSK 
■ LP I = 1 2 
CDR = 9 
USK=1 

IF(NSLICE) 120,20,60 

INPUT NUMBER OF DATA SLICES AND READ DATA SLICE INFO FOR PRESENT RUN 

20 RE AU (COR > 30 ) NSLICE 
30 FORMAT CIS) 

wRI Tt (LPT > 35) NSLICE 

3b FORMAT l 1 H 1 , ' **** OATA LARDS ****',///,' NO. OF DATA SLICES- ' > IS, / 
3//,SX, 'DATA SLICES' ,//) 

IF (NSLICE. GT. 200) GO TO 100 
IFCNSLICt.LT. 1) GO TO 120 
DO SO .J = l, NSLICE 
kEAU(CDR, 40) (iCARU(J,K),K=l,l5) 

40 FORMAT ( IX, A2, 1413) 

NR I TE (LPT, 40) (ICARD(J,K),K = 1, IS) 

50 LUNTINUE 

CALL FOPEN (DSK, "ICARD ",6000) 

WRITE BINARY (DSK) ICARD 
CALL CLOSE (DSK, IER) 

GO TO 200 

INITIALIZE VARIABLES AND SWITCHES FOR EACH DATA SLICE 

60 CALL FOPEN ( DSK , " I C ARU " , 600 0 ) 

READ BINARY (DSK) ICARD 
CALL CLOSt (DSK, IER) 

IVARSrt(l )= ICARD (I, 2) 

IVARSrt(2)=ICARD(I,l) 

START=3600.»FLO«T(ICARO(i,3))+60.*rLOAT(ICARD(I,4))+rLOAT(ICARD(I, 
oiS) ) 

1VARSN(3)=ICAR0(I,B) 

IVARS.v(4) = ICAR0(I,7) 

1VARSn(S)=ICARO(I,8) 

I V A K: S >1 ( 6 ) = I C A R 0 ( 1 , 9 ) 

IVARSNC7)=ICAH0(1, 10) 

IvARSv. (I3) = ICAR0(I, 11) 
lH£AO(lo)=ICARO(I,3) 

IHEAD(17)=ICARD(I,4) 

InEADC 1 B)sICaR0(I,S) 

IHEADC 19)=ICARd(I -7) 

IF ( I VARSw ( 1 ; -2) 90,90,200 
90 IHEAO(20) = 1CARD( I , 12) 

GO TO 200 


C 


C E N R U R MESSAGES 

C 

100 VJUTECLHT, 110) N3LICE 

llu FORMA f C /// / 9 FATAL ERROR • • • READ NUMBER OF DATA SLICES TO BEM6/'. 

h CaM NOT HAVE FORE I HAH 200.') 

3 F UR 

120 *-.RIlfc(LPl# lio) N SLICE 

130 FORMAT (/// # ' FATAL ERROR • . .READ U OF DATA CANOS 10 BE NEGATIVE. 

»* OF SLICES= Mo) 

STOP 

200 IF Cl .EO.NSLICE) CALL Dr ILlfl ( " IC ARO'MfcR ) 

RETURN 

END 


ORIGINAL P ACr ro 

0F **» <25®* 
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OVERLAY UVP 

SUBROUTINE TINPUf (UTYPE, 1 fFIlfc,F ft. JM,START,N5ECS, I) 

SUBROUTINE TINPUT Mb ADS F URCE A n\J EMG DATA IN FROM TAPE. 

CALLING ARG'J^NISS D1YPE - OnfM TYPE 

ITFILE- FILE CURRENTLY BEING ACCESSED ON TAPE 

FILNUM- FILE ON TAPE fcHERE DESIRED DATA IS LOLA. 

START - TIME IN TOTALED SECONDS OF BEGINNING OF 
DESIRED DATA SLICE 

NSECS - NUMBER OF SECONDS OF DATA TO BE HEAD . 

I - DATA SLICE NUMBER 

INTEGER FILNUM,OTYPE 

DIMENSION I ARRAY (536) , ITFlLE(4) , I TIME (3) 

COMMON DArUMS(4S00) # IHEADC20) #CAL(2#4) #IVARSft(iG) 

LU‘iMON/IGD/ MUSCLEC 7, cl) ,IFQRR(2) 

EQUIVALENCE ( I ARRAY# DATUMS (4230) ) / ( I T IME# I ARRAY (2) ) 

LP T = 1 2 
MTI = i 

I RE aD= i 0 30K 
IGOOOsiiiOSK 
IBREC=40000K 

CHECK TO SEE IF DESIRED FILE IS CURRENTLY BEING ACCESSED. 

IF NUT, EXCEPT UN FIRST DATA SLICE CLUSE PRESENT FILE BEFORE 
OPEN TNG NE.Y ONE. 

IF(I.EQ.l) GO TO 1 

IF* (F ILNUM .EQ .ITFILEC3)) GO TO 5 

CALL CLOSE (MTI,I6R) 

1 ITFILE(l)s»rtT M 
ITFILE(2) = "0S " 

1 TF IlE ( 3 ) sF ILNUM 
ITFILt (4 J ="<0><U> M 

POSITION TAPE AT THE BEGINNING OF ThE DESIRED DATA FILE. 

CALL MTOPD (MTI#ITFILE#0#XER) 

IF(IER.Nt.l) GO TO 300 

POSITION TAPE AT THE BEGINNING OF THE DATA SLICE. 
b CALL FIND (FILNUM# START) 

DECIDE rJHERE TO GO FOR EACH DATA TYPE. 

GO TO (10,50,50, 10), DTYPE 
ACQUIRE FORCE CAL UR FORCE DATA 


10 JJ = l 

NRtCS=2*N3ECS 
DO 30 J=1,NREC3 

CALL MTDIO (MTI,IREAD,IARRAY#ISTAT#IER#NrtRDS) 

IF (IER.NE. 1 .0R.ISTAT.NE.IG00D.0R.NWRDS.NE.535) GO TO 200 

DO 20 K=27 ,555,51 

DATUMS ( JJ) =FLOAT ( 1 ARRAY (K) )/8. 
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JJ=JJ+1 
20 CONTINUE 
30 CONTINUE 

It- (DTYPE.EG. 1 ) GO T ) 1000 
GO TO I5u 

AC NO I HE EMG CAL OK EMG OATA 
50 JJ=1 

NREC3=2*NSEC3+1 
00 BO J=1,NREC3 

CALL MTOIO (f!TI,IR£AD,IARHAY,ISTAT,ItR,i'H'iROS) 
IFCltK.NE.t.OH.ISTAr.NE.IGOUO.OR.NwRDS.NE.535) GO TO 200 
00 70 K=t , 10 
L-25+K+50* (K-l ) 

Lb0=L+50 

00 60 KK=L » L50 

IF (AK.EQ.Ltl ) GO TO 60 

DATUMS (JJ) =FL0AT(IAHKAY(KK))/8, 

JJ=JJ+1 
60 CONTINUE 
70 CONTINUE 
60 CONTINUE 

IF CUTYPE.EQ.2) GO TO 1000 
150 IUACR=i8RtC+2*NSECS+2 

CALL MTOIO (MTI,IBACK,IARRAY,ISTAT, IER) 

GU TO 1000 

ERROR MESSAGES 

200 WRITE (.LPT ,210) OT YPE , F ILNUM , I ST AT , IER , Il.vKOS , I TIME 
210 FORMAT!///, " FATAL ERROR .. .DURING INRUi' OF DATA TYPE", 12, "IN FILE 
it" , A2 , / , 1 OX , " T APE STATUS (8) = ",018,/, 1 OX, "FORTRAN ERROR = " ,17,/, 1 OX , " 
o)« Ur WORDS RE AO = " , 1 6 , / , 1 OX , "U AT A T IME = " , OX , 3 I i ) 

STOP 

300 wRITECLPT ,310) IER.FILNUM 

310 FORMAT ( /// , " FATAL ERROR ... FORTRAN ERROR CODE = ",I3," DURING SEARCH 
o, FOR FILE " » 42) , 

STOP 

1000 CONTINUE 
RETURN 
END 
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SUBROUTINE FIND (FlLNUNt* $ TART) 


SUfiHUUT INE KIND POSITIONS THE TAPE AT THE BEGINNING OF A PAT A SUCK. 

CALLING ARGUMENTS! FILNUM- FIlI: UN I'APt v.MERE DATA SLICE IS LOLA I 

START - TIME IN TOTALED SECONDS OF Thk bEGIivNlf. 
DESIRED DATA SLICE. 


DIMENSION I ARRAY (536) 

INTEGER FILNUM 
NtAL NOrt 

COMMON DATUMS (4500) , I Hfc AD ( 2 0 ) , C AL ( 2 ,4 ) , I V ARSft ( 1 0 ) 

COMMON/ 1 GO/ MU3CLE(7,4) # IFURR(2) 

EQUIVALENCE (IAKnAY,DATUMS(4230) ) 

IPASSsO 
LPT =12 
M T 1=3 

IEOF=-7327iK 
IFFILE=30000K 
Itfr ILE=40000K 
IkE AD= 1 030K 

READ EACH RECORD AND COMPARE ITS TIME LABEL WITH THE START TIME. 

10 Call MTDIO (m1I,IRLAD,IAkRAY,ISTAT,IER) 

IF CISTAT.EO. I EOF ) GO TO 20 

NOv’.s3600.*At33(FLOA[ (IARRAYC2) ) ) +60 . * ABS C FLOA T ( I ARR A Y ( 3 ) ) ) ♦ ABS (FLOA 
a)T ( I ARRAY ( 4 ) ) ) 

IF (NOrt.NE. START) GO TO 10 

ACQUIRE THE HEADER INFORMATION FOR THE PRESENT DATA SLICE. 

DO 15 J = l, 15 
IHEaO(vI) = IARRAY ( J ■* 5 ) 

15 CONTINUE 

I6REC= IBP ILEtl 

CALL MTDIO (MTI , IOREC , I ARRAY# ISTAT , IEM) 

GO TO 100 

DID NUT FIND THE START TIME BEFORE ENCOUNTERING AN EOF. IF ONLY ONE P 
HAS BEEN MADE THRU THE DATA, BACK UP AND TRY AGAIN. OTHERWISE# STOP) 

20 IFCIPASS.EO.l) GO TO 40 
DO 30 J = 1 , 2 

40 CALL MTDIO ( MT I , I DF HE , I ARR A Y , 1ST A T , I ER ) 

CALL MTDIO (MTI, IFF ILE, l ARRAY , ISTAT, IER) 

IP ASS= 1 
GO TO 10 

40 ARI TE(LPT,50) FILNUM 

50 FORMAT (///, ' CANNOT FIND START TIME IN FILE # ,A2) 

STOP 

100 CONTINUE 
RETURN 
END 


i 
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SUbi<OU r iNt UUJSO (DTYPE,NGLCa,i r ILT j.-.J 

SUBkOU TINE OUT"0 IS A 01GII4L NOTCH FILTER. DEPENDING ON THE DATA CAR 
RtulJEG I / ONLY GO Hi IS REMOVED OH ALL HARMONICS OF SO HZ OH TO 360 

CALLING AliOcT'ivY f S S * DTYl’i. •• TYPE Ul DATA 

HatCS - NUMbtK OF SLCS OF CAIA lu ut FILTERED 

FILTSM- Sy.ITCH WHICH INDICATES WHETHER ONLY SO 
IS TO HE REMOvEO OH ALL HARMONICS OF 60 


INTEGER D T YPE » FILTSM 

COMMON DA TOMS (0300 ) t IHE AD (20) , C AC (2, 4 ) , I V AH8W ( 1 0 ) 
COMMON/ IGO/ HOSCLE ( 7 » 0 ) , IF0rtH(2) 
ton VALENCE ( IRATE, IHE AO ( ID) 

L00PS=1 

T.',0PI = 6. 28310530717 
IF CFILTSY..EQ.1) L00FS = 5 
t.PTS = 102A*NSECS 
DO 30 J = 1 » LOOPS » 2 
A = 0.0 
t'*0 » 0 

ivs(60*J)*r«0PI /FLOAT (IRATE) 

00 10 K=1,NPT3 
AsAtOATUMS(K) *CUS(K*W) 
c = 6 + 0A f UMS(K) *3IN(K*W) 

10 CONTINUE 

A=A*2. /FLOAT (NPTS) 

BsB*2 . /float (npts) 

00 20 K=1,NPTS 

DATUMS (K) = 0ATUM8(K)-A*C0S(R»H)-6'*3IN(K*I1) 

20 CONTINUE 
30 CONTINUE 
RETURN 
END 


ORIGINAL PAG,: it 
DP POOR QUal TV 
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OVERLAY UV3 

SUbKOUT IKE 1 >Cal(D TYPE* MSECS* IF GAIN) 

SUbROuT lNE OCAL CALCULATES THE SCALE FACTORS r OR THE FORCE AND EMG DA i 

TMc CALIttRAriO J CONS IAN I FOR .INTEGRA f HO AREA IS ALSO FOUND • 

Tub CAL D M a\ u) rni-.l-iLOKVt.MLD INTO UNITS UF FuUNDS OH M IC HU VOLTS F G;< 
PLOTTING ON THE COMPUTER T E H M INAL. 

CALLING ARGUMENTS: DTYPE - TYPE OF DATA 

MSECS - LENGTH OF CAL DATA 

IFGAIM- GAIN APPLIED TO FORCE CAL SIGNAL 


DIMENSION MIGH( 10) * LOW ( 10) * POUNDS (2600) 

INTEGER DTYPE* FISCAL 
RtAL LUfiSUM 

COMMON lMTU*SU500) * IHEAD<20) *CAL(2*4) * IVARSW( 10) 

COMMON/ I GO/ MUSCLE (7*4) * IFORR(2) 

EQUIVALENCE (FLBCAL * I HE AD ( 1 S) ) * ( IRATE* THE AD ( 1 1 ) ) 

EQUIVALENCE (MAGTUD* IMEAD ( 13) ) * (POUNDS * DATUMS C 1 ) ) 

LPT=12 

IF (DTYPE. EQ. 2) GO TO 100 

FIND FORCE CALIBRATION VARIABLES 

C A L ( 1 * l)=FL0Ar (FLBCAL) 

C A L ( 1 * 4 ) = F L U A T (IF GAIN) 

LOO^ FUR JUMP IN DATA TO INDICATE FORCE CAL. 

0 E G I N = D A T U MS ( 8 0 ) 

NPTS=2u*NSEC3-80 
DO 10 Js8i*NPT3 

IF (DATUMS (J) «GT •GEGIN + 500 • ) GO TO 30 
10 CONTINUE 

WRITE (LPT/20) 

20 FORMAT ( /// * 'COULD NOT FIND THE FORCE CAL. ChECK DATA SLICE TIMES 

a)') 

STOP 

FOUND THE JUMP, NOW COMPUTE THE AVERAGE BEFORE AND AFTER THE JUMP. 
CAL WEIGHT IN COUNTS = AVERAGE AFTER - AVERAGE BEFORE 
FORCE SLOPE = COUNTS/POUND 

ZERO FORCE = AVERAGE COUNTS JUST BEFORE CAL JUMP 
30 MISUM-O.O 

LOw-SUMsO.O 

DO 40 K=l/20 

HISUM = HISUM + 0ATIJMS( J+K + 60) 

LOwSUM=LONSUMtOATUMS( J-K-60) 

40 CONTINUE 

FOkSLHs (H i SUM/20. -LOW5UM/20.) /FLOAT (FLBCAL) 

CAL(1*2)=FURSLP 
ZtROFORsLOv’/ SUM/20, 

CAL ( 1 * 3 ) sZERQFOR 

SCALE FORCE CALIBRATION DATA 

NP TS~20*NbEC3 
UO BO J=l*NPTS 

PGU L?;3 ( J ) = ( DATUMS ( J ) -ZEROFQR ) /FORSLP 
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, ■_ I i i u. _ i - j t * - - r - i . u— r ~ w - » 

5 0 CONTINUf: 

(iU TU dOO 

t 

F I NO tMG CALIBRATION VARIABLES 

FIRST REMOVE A >NY DC OFFSET IN THE CAL DATA 

100 SUH-O.O 

INDtX= 1000 *NSEC 3 
DO no J=l, INDEX 
buMsSUMtDATJMS(J) 

110 CONTINUE 

UFF 3 ET=SUM/FL 0 AT ( INDEX ) 

DO 120 J = l, INDEX 
DATUMS U)=DAruMS(J) -OFFSET 
120 CONTINUE 

CALCULATE STANDARD DEVIATION AND OBTAIN SCALE FACTOR. FOR ZERO MEAN, 

stationary s i dial standard deviatiun=rms value 

XSGUAR = 0 .0 

do iro j = i, index 

XSuOAR = XSQUaR+DA TUMS ( J ) **2 
1 R 0 CONTINUE 

EMGCAL=SQn T(XSOUAR/( INDEX- 1 n /FLOAT (MAGTUD) 

CAH 2 , I)=f LOATMAGTUD) 

CAL ( 2 , 2 ) =EMGCAL 

SCALE CAL DATA FOR PLOTTING. 

DO ISO J = 1 > INDEX 
DA f Ui-'S ( J ) sOA TUMS ( J ) /EMGCAL 
ISO CONTINUE 
200 RETURN 
END 
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OVERLAY 0V4 

SUBROUTINE GRAPH (0 ! YPE , 31 AH T , IlStCS, SP All , 31 DERR ) 

SUBROUTINE GRAPH PLOTS ON f HL COMPUTER TERM INAL THE OATA SLICE UEING 

PROCESSED ALONG n IT h VAHIOUS HtAOER t NFOft *• A I ION • 


CALLING ARGUMENTS! 01 YPt 


TYPE OP OATA 


START - BEGIN TIME OP PRESENT DATA SLICE 

NSEC3 - LENGTH OF OATA SLICE IN SECONDS 

SPAN - AVERAGING INTERVAL FOR PSO 

STDERR- NORMALIZED STANDARD ERROR FOR PSD ESTIMA 

DIMENSION DATE (3) , POUNDS (2600) ,PSDN(410) 

OIVENSION FREQ (4 10) , CAL OAT A (2000) 

INTEGER OTYPE 

COMMON DATU>!S(4500), IHEAD(20),CAL(2,4), IVARSa(IO) 

COMMON/ IGO/ MUSCLE (7,4) , IF0RRC2) 

EQUIVALENCE (ROUNDS, DATUMSC 1 ) ) , ( PSON , DATUMS ( 4000 ) ) 

EQUIVALENCE (FREQ, DA TUMS (iSOO) ) , ( P LOCAL , C AL ( 1 , 1 ) ) 

EQUIVALENCE (CAL DATA, DATUM 3(1) ) , (EMGMAX , DA TUMS ( 4499 ) ) 

EQUIVALENCE (FORMAX , DATUMS ( 4500 ) ) 

RETURN 

END 




OVERLAY UV5 

SUBROUTINE EMG (NotCS , VOL TSEC) 

SUbRUU TINE EMG SCALES EMG DATA FOR PLOTTING ON T Mb COMPUTE* AND FINDS 

Ink INTEGRA fuu t»’U VAI.UK FuW l';t DATA LI. • 0 K • 

CALLING ARGUMENT S t Notes - LENGTH UK DATA SLICE 

VuL TSEC- INTEGRATED EMG VALUE IN M I CPU vOL T * SEC • 

COMMON DATUMS (4500) ,IHEAD(20) /CAL (2, 4) ^ IVARSv. ( 10) 

COMMON/ I GO/ MUSCLE ( 7 , 4) , IK JRR(2) 

EQUIVALENCE (IRATE, 1 HE ADC 1 1 )) , (EMGCAL/ CAL(2, 2) ) 

EQUIVALENCE (E v l UMAX, DATUM $(4499) ) , (EMGVAR , DA TUMS ( 4500 ) ) 

CALCULATE OFFSET, SUbTRACT FROM DATA AND APPLY SCALE FACTOR 


5 U M s 0 • 0 

nPTS=I 024*NSEC3 
DU 10 J=l,NPTS 
SUM=SUM+DATUMS( J) 

CONTINUE 

OFF SET -SUM/FLO A T (NPTS) 

CAL (2,4) =OFFSE T 
EMGMAXsO.O 
DO 20 J=l, NPTS 

DATUMS (J)= (DA TUM3 (J ) -OFFSET ) /EMGCAL 
IF (AiiS(DATUM3(J) ) • G T . EMGMAX ) EMGM A X = D A TUMS ( J ) 

CONTINUE 

FIND .INTEGRATED EMG VALUE 

NPTS=NSECS*i024-i 
Hslt/FLUAT (IRATE) 

SUMOrO.O 

SUME = 0 • 0 

NPTS2=NPTS-2 

DO 40 I = 1 , NPTS2 , 2 

SUME=SUME + A6S (DATUMS (1 + I) ) 

SUMOsSUMO+ASSCDATUMSCI) ) 

CUN 1 INUE 

VOL T StC = (2 ,0*SUMO+4.0*SUME- ASS (DATUMS ( l )) +A8S (DATUMS (NPTS) ) )*H/3.0 

CALCULATE THE VARIANCE OF THE DATA 

EMGVARsO.O 

DO 50 J = 1 , M P T S 

EMGVARs EMG VAR FDA TUMS ( J) **2 

CONTINUE 

EMGVaR=EMGVAR/FLUAT (NPTS- i) 

RETURN 

END 
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bUbKUUII.it FUKCE (OfYPE, MSECS, i P A > ; , H U G « 1 1-1 ) 

SUUHOUT INE FORCE SCALES ALL FORCE DAfA FUR PLOTTING ON THE COMPUTER A 
AVERAGES FORCE OAfA 10VE* INTERVALS DEUR.'iMEO BY DATA CARO HEOUtST 1 
FOR PR I N I EL OUTPUT. 

CALLING ARGUMENTS; D1TPE - TYPt OF DATA 

NSECS - LENGTH OF DATA SLICE 

SPAR - AVERAGING INTERVAL FOR PRINTED OUTPUT 

FDGAUm- AMOUNT OF GAIN WHICH WAS APPLIED TO FORi 
DATA 

DIMENSION POUNDS ( 2600 ) , AVGLBS(liO) 

INTEGER OTYPE 

COMMON DATUMS (4S 00) , I HE AO ( 20 ) , C AL ( 2 , 4) , I V ARSh ( 1 0 ) 

COMMON/ 1 GO/ MUSCLE ( 7 , 4 ) , IFORR(2) 

EQUIVALENCE ( POUMUS » DAT UK’S ( 1 ) ) , ( A VGLBS , DATUMS ( 4000 ) ) 

EQUIVALENCE (FCGAI(i,CAL( 1 ,4 ) ) , UERUFOR,CAL ( 1 , 3) ) 

EQUIVALENCE (FORSLP, CALC 1 ,2) ) , (FORMAX , DATUMS (450 0 ) ) 

GAIN = FCGAI,J/FDGAIN 

FIND ZERO BASELINE BEFORE BEGINNING OF FORCE DATA 

SUM=0.0 
DO iO J = l , 2 0 
SOKsSOM+POONOS ( J ) 

30 CONTINUE 

BASEsSUM/20. 

SCALE FORCE OATA INTO POUNOS 

FORMAX=0.0 
NPT3=2Q*NSEC3 
DO 40 Jsl/NPTS 

POUNDS ( J)=(POUNDS(J)>BASE)*GAIN/FORSLP 
IF ( POUNDS ( J ). G T . FORMAX ) FORMAX-POUNOS ( J) 

40 CONTINUE 

AVERAGE FORCE DATA OVER SPAN INTERVAL 

L00PS= MSECS/ IF IX(SPAN) 

NPIS=20*IFIX(SPAN) 

DO 60 J= l, LOOPS 
SUM=0.0 

DO SO K=1,NPTS 

SUM=SUM+P0UN03( (J-l) *NP TS+K ) 

50 CONTINUE 

AVGL8S(J)=SOM/FL0AT(NPTS) 

60 CONTINUE 
RETURN 
END • 


QUALITY 
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OVERLAY (JvG 

SUBROUT I tit F F T P S D (NSECS, SPAN, PSD ■ .STUF.KR) 

SUBROUTINE FF'fPSU T A K 1 5 THE FOURIER TRANSFORM OF EMG DATA /. l.y CALCULAi 
THE RA.T POvFR 3PECTKAI. DENSITY UP TO nno H7. ACCORDING • IhF. OAT a 
wEliUESI ' ’UO THEO P.’iO .1 •> !. Al.i.. ».* A T Ei) •'•Uji’.l' ;■■.»:> li.oiiivi 
PROPER I ItS . IU REFLECT i MOSt SIAIISliCAL PROrLMIES I Hi: ,i, iAL WED 

STAN0AR0 ERRUR IS ALSO calculated. 

CALLING ARGUMENTS: NSECS - LENGTH OF DATA SLICE 

SPAN - BANDWIDTH OF SMOOTHED PSD 

PSOSUM- AREA UNDER POWER SPECTRUM 

STOERR- NORMALIZED STANOARO ERROR OF PSD 

DIMENSION PSD (50 '48) ,PSDN(A10) , PERCNT ( <1 1 0 ) 

DIMENSION CUMPC f < a 1 0 ) ,S( 1024) 

COMMON DATUMS («SOO) , I HE AD ( 20 ) , C AL ( 2 , 'I ) , I V ARSW ( 1 0 ) 

CUM-iON/IUD/ MUSCLE ( 7 » u ) , IF QRR ( 2 ) 

EQUIVALENCE (P30,0AT0MS(1) ) 

EQUIVALENCE (IRATE, IhEAOdl )) , (PSDN , DATUMS ( U000 ) ) 

EQUIVALENCE ( PERCNT # uA T UMS ( SO 00 ) ) , ( CUMPC T , DA TOMS ( 3500 ) ) 

EuU I VALENCE ( ExPVAL, DA TOMS ( 2500 ) ) , (STDEV# DAT UMS (2501 ) ) 

CONTINUE 

APPLY COSINE NlNCOw TO DATA 
CALL WINDOW (NSECS) 

DETERMINE THE POwER OF 2 AND CORRESPONDING NUMBER OF DATA POINTS 

N£XP=I0+NSEC3/2 
IF (NEXP.EG. 14) NEXPs 1 3 
N=2**NEXP 
P1S=FL0AI(N) 

TAKE FFT OF THE DATA 

CALL RFFT (DATUMS, NEXP # S # - 1 , IFERR ) 

CALCULATE RAw POWER SPECTRAL DENSITY 

H=1 ./FLOAT (IRATE) 

K = 1 

N 1 =N+ 1 

00 10 J=3,Nl,2 

PSD(K)=8.*H/P rS*(0ATUMS(J)**2+0ATUMS(J+l)**2) 

K-K + 1 

10 CONTINUE 

CORRECT PSD AMPLITUDES FOR WINDOW REDUCTION 
N2»N/2 

DO 15 J = 1 / n 2 

PSD ( J ) = 1 . 1'429*PS0(J) 

15 CONTINUE 

OUTPUT RAw PSD TO TAPE 

CALL PSDSAVE (NSECS) 
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C CAlCULAlt SMOOTHED PSD 

C 

FOsU/(PrS/FLOAr(IHArt)) 

IUTA = 1F I X ( SPAT /FO) 

SPANhFLUAI dill A ) * F 0 
HSU^AXn»j,0 
K = 1 

N2=Na/IOTA*IOrA 

DO 50 J=l, N2, IOTA 

TtMpsO.O 

DU 20 JJ=l,IOTA 

TEMP=IEMP+P30(JtJJ-l) 

20 CONTINUE 

PSD CK)=ThMP/FLOAT( IOTA) 

IF IFLUATCJ/IOT A) iSPAN.GT.^OO.+SPAN) GO TO 25 
IF (PSD C K ) . G T .PSD MAX) P$OMAX=PSD (K ) 

25 K=Ktl 
30 CONTINUE 
C 

C CALCULATE THE NORMALIZED PSD 

C 

LIM=IF IX(400«/3PAN)+l 
DO 40 J=lrLlM 
PSi)N(J) =PSD( J) /PSDMAX 
40 CONTINUE 
C 

C INTEGRATE THE SMOOTHED PSD 

C 

PSDSUMsO.O 
DO SO J = WLIM 
PSDSUM=PSD3UM+SPAN*PSD( J) 

50 CONTINUE 
C 

C CALCULATE X EACH BANDWIDTH IS OF TOTAL AND FIND THE CUMLATIVE % OF TOT 

C 

DO 60 J=1,HM 

• PEPCNT ( J )=SPAN*PSD( J)/PSDSUM*iOO. 

IF ( J ,NE . I ) GO TO 55 
CUMPCT (l)sPEHCNT(l) 

GO TO 60 

55 CUMPCT ( J)=CUMPCT ( J~ 1 ) tPERCNT(J) 

60 CONTINUE 
C 

C CALCULATE THE EXPECTED VALUE AND VARIANCE OF THE PSD 

C 

FsSPAN/2. 

EXPvAL=0.0 
DO 70 J=1,LIM 

EXPVAL-EXPVAL-KPEPCNT (J)/100. ) *F 
F=F+SPAN 
70 CONTINUE 
F=SPAN/2. 

V ARV AL = 0 .0 
DU SO Jsl/LIM 

VAWVAL=VaRVAL>(PERCNT(J)/100*)*(F-EXPVAL)**2 

F=FfSPAN 
80 CONTINUE 

STDEV=SORT (VARVAL) 

C 

c find the normalized standard error 
c 

STDERR=SQRT( 1 ./FLOAT (IOTA)) 

RETURN 
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SUBROUTINE rtlNoOU (Nbb.CS) 


SUtthOUT I Nil fllNDO* APPLIES A COSli.i. TAPER TO THE FIRST AND l A > T TENTHS 
THE DA I A TO REDUCE Lt AK AGt - 

C CALLING ARGuiEufSj NS ECS - LENGTH 01* fHb DATA 

C 

COMMON OATUMSC-4SOO) , I HE AD ( <2 0 ) , CALC 2 , 4) , l V APS - ( 10) 

LOMs-ioN/iGo/ muscle (7*4) # ifukncz) 

PI=3. 141S927 
NPTS=1024*N3ECS 
IEDGEsNPTS/10 
f TOT AL = 1 .024 AFLOAT (NSECS) 

K*1 

L'O 10 J=1 , i EDGE 
Tl=.001*FLQAr(J) 

C1k.5*(1.-C03(PX*U/(.1*T TOTAL)) ) 

LATUMS(J)-C1^UATUMS(J) 

DATUMSCNPTS-J>i)=Cl*0ATUMS(NPT3-J+l) 

10 CONTINUE 
RETURN 
END 


o o o 


SUbROUl Int KFFT (A/'WSjlFS, IFLHIT) 

C UNt-DIMEivSIUNxL RtAL FINITE FOURIER TITANS! OHM 


FOURIER TRANSFORM SUHHOUT INF. FUR REAL DATA. 

C PROGRAMMED t ■> GYSTEM/ibO, PROGRAMMING SUPPORT# 

C CUMKan IV# (ofc:. FORM C2 j-6G0'O . 

C I hlS UECK IS St T UP FUR IbSYS ON THE I tlM 7094 

C 

C This PROGRAM U3K3 THt SUbROUTINF. FFT TO COMPUTt COMPUtX 

C FOUKItR TRANSFORMS OF Rt At DATA. PR FURT S.D.A. NO. 3465 IS 

C AVAILABLE THROUGH SHARt. 

C 

C THt FOURItR StRIES IS 

C 

C X ( J J - SUM OVER K = 0 TO N « UF C CK ) *F.XP ( 2 *P I * I * J* K/N ) 

C J=0»1,2»...#N-1 

C 

C WHERE I =S(4Rf ( - 1 ) AND WHERE CCK) IS CUMPLEX. 

C SINCE X ( J ) IS REAL# C(K) - CON JG (C ( N-K ) ). I HEREFORE ONE.' 

C C(K)#K= 0, 1 , . . . ,N/2 ARE COMPUTED AND/OR USED. 

C 

C ARGUMENTS- 

C 

C A IS INITIALLY THE INPUT ARRAY# X, WHEN COMPUTING A FOURIER 

C TRANSFORM AND C WHEN COMPUTING A FOURIER SERIES. A IS REPLACES BY 

C THE OUTPUT ARRAY# C IN THE FORMER CASE, X IN THE LATTER. 

C THE X VECTOR CONTAINS THt REAL DATA X C 0 ) , X (1 ) , . . . , X l N- 1 ) 

C THE C VECTOR CONTAINS THE COMPLEX’ FOURIER AMPLITUDES 

C C ( 0 ) , C ( 1 ) # . . . # C (N/2 ) . THt COMPLEX VECTOR (. IS STORED ACCORDING 

C TO Tht NORMAL FORTRAN IV CONVENTION FUR STORING COMPLEX NUMBERS. 

C I.E.# ‘REAL PARTS IN ALTERNATE CELLS STARTING ITh THE FIRST, 

L IMAGINARY PARTS IN ALTERNATE CELLS STARTING WITH THE StCONO. 

C 10 ADHERE TO FORTRAN RULES# X(0)# XC1),..., ARE REFERRED To AS 

C X ( 1) # X (2) , . . . # RESP. IN THE PROGRAMS. ALSO, C ( 0 ) , C ( 1 ) # . . . ARE 

C REFERRED TO AS C ( 1 ) # C (2 ) , . . . , RESP., IF C IS DESIGNATED AS 

C COMPLEX IN A TYPE STATEMENT. 

C 

C M GIVES N=2**M 

C 

C THE ARGUMENTS S# IFS# AND IFERR ARE THE 3A E AS IN THE 

C SUBROUTINE FFT AND THE USER IS REFERRED TO THE COMMENT CAROS 

C IN FFT FOR THEIR EXPLANATION. 

C DIMENSION STATEMENTS- THE DIMENSIONS OF ARRAYS A AND S SHOULD 

C BE N + 2 AND N/‘T# RESP. FOR THE LARGEST N TO EE USED. FOR 

C EXAMPLE# IF THE LARGEST M IS 13# THEN, N=8192 AND ONE SHOULD 

C HAVE THE DIMENSION STATEmENT- 

C DIMENSION A ( 0 1 9 4 ) , SC2048) 

C IF ONE WISHES TO SPECIFY A TO BE COMPLEX BY A TYPE STATEMENT# 

C ONE SHOULD GIVE IT A DIMENSION OF N/2 >1 # FOR THE LARGtST N. 

C 

C 

DIMENSION AC4500) ,S(1024) 

IF ERRS = 0 
N = 2 * * M 
NV 2 = N / 2 
NV4M1 = N/4 - 1 
MM 1 = M - 1 

IF (IABS(IFS)“l) 40,40,20 
20 IF {Mp-MJ 30,30, GO 
30 IFERRS = 1 
40 NP = N 
MP = M 
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CALL. PM ( A#to#S#0# IFP.MKt ) 
IPtRRS n IFtrtMS + IP'ERRl 
SO KO = nH / N 
KT = KD 
NPV'4 n ,*4P / a 
IP UPS) ? >Q t 80 9 90 
60 CmLL PPT(a, ''il,S,~2UPtRR2) 
IFERRS = IFtRRS ♦ IPERR2 
do 7o K3i,NV'Hi 
J=NV2-K 

A 1 K= A(2*K+i) t A(2*J+1) 

A1 I=A(2*K + 2)-A(2*J + 2) 
A2RsA(2*K*2)+A(2* J + 2) 
A21sA(2*Jtl )-A(2*K + l) 
kk r s NPva-Kr 

A*RsA2R*3(KKT) + A2I*S(KT) 
Aft 1 = A2I*3(KKT)- A2R*S(KT ) 
A(2*Kfl)s(Alrt + A»4R)/<4. 

A(2*K + 2) = (AUtAwI)/4. 

A(2*J + l)-(AlR-A,‘lK)/4. 
A(2*J+2)=( A /J I - A 1 I)/a. 

70 KT=KT+KD 
T = A ( 1 ) 

A(l)=(T+A(d))/2. 

A ( N 1 1 ) •- (T-A(2))/2. 

A(2)=0. 

A(N+2J = 0. 

A(NV2tl) = .5*A(NV2tl) 
A(toV2t2) = -,S * A(iMV2t2) 

60 IFERR = IFERR3 
RETURN 

90 DO i OO K = 1 ,NV4M1 
J = N v 2 - K 

AlR=A(2*Kti) + A ( 2 * J + 1 ) 

AlI = A(2*K + 2)-A(2* J + 2 ) 
A«RsA(2*K+l)-A(2*JfJ) 
AwI=A(2*K+2)+A(2*J+2) 

KKT = i\HVa - K T 
A2RsAU**S(KKT) - A,U*S(KT) 
A2I=AftR*3(KT) ♦ AftI*S(KKT) 

A ( 2 * K + 1 ) = AIR - A2I 
A(2*K*2) = All + A2 k 
A ( 2* J+ 1 ) ^ AIR * A2I 
A(2*J+2) = A2R - All 
100 KT = KT t KD 
T ^ A ( 1 ) 

A ( 1 ) = T + ACN+1) 

A(2) = T - A (N + 1 ) 

A ( N V 2+ 1 ) = 2.*A(NV2ti) 
ACNV2+2) s-2.*A(NV2t2) 

CALL FF r(A,^Ml,S,2,IFERR2) 
IP ERRS = I P ERRS* IFERR2 
bU TO 80 
END 
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SUBROUT INt F F r ( A - 4,8, TFS, IFtHR) 

FuW, ONfc-OlMEHSIONAL FINITE RLEX FOURIER TRANSFORM, 


FOURIEH THANSFQRM SUBROUTINE, PPIHRAMMFIO IN SYSTKM/J60, 

IJ«3lC P s' Ul.RA : ti j'i SC'PpTiHI# iiMl'O'J IV, m.’Ri' L.M 0 4 
1HI3 Deck SLI Of’ FUR IbSYG 0,1 iBM 7094. 

DUES EITHER FOURIER SYNTHESIS# I .E . #C0MPU TtS COMPLEX F0U"tKR ScRIES 
GIVEN A VECTOR OF N COPPLEX FOURIER A!PL I TliUtS# OH # GIVEN A vcCTOR 
OF COMPLEX DATA X DOES FOURIER ANALYSIS# COMPUTING AMPLITUDES. 

A IS A CUMPLEX VECTOR OF LENi.TH N = 2 » * M COMPLEX NOS. OH 2*N HEAL 
NUMBERS . A IS TO UK SET BY USER, 

V IS AN INTEGER 0 . L T ,M. LE . I J , SET BY USER. 

S IS A VECTOR S ( J J - SlN(2*PI*J/NP ), J = l , l , . . . . , NP/4- 1 , 

COMPUTED BY PROGRAM. 

IF S IS A PARAMETER TO BE SET bY USER AS FOLLOwS- 
IFS = 0 TJ SET NP=2**M AND SET OP SINE TABLE. 

IF S= 1 TO SET N=NP=2**M, SET UP SIN TABLE# AND 00 FOURIER 
SYNTHESIS, REPLACING TmE VECTOR A bY 

X(J) = SUM OVER K=0#N-1 OF A(K)*EXP(2*PI*I/N)**(J*K) # 

J = 0#N-1, .(HERE I -SORT C- 1 ) 

THE X'S ARE STORED WITH RE X(J) IN CELL 2*J+1 
AND IM X(J) IN CELL 2*0 + 2 FOR J=0 , 1 , 2 # . . . #N- 1 . 

THE A'S ARE STORED IN THE SAME MANNER . 

IF3 = -1 TO SET N=NP=2**M#SET UP SIN TABLE, AND 00 FOURIER 

ANALYSIS, TAKING THE INPUT VECTOR A AS X AND 
REPLACING iT bY THE A SATISFYING THE ABOVE FOURIER SERIES. 

IFS>0 * SET UP SIN TABLE AND RETURN 

IFS=+2 TO UO FOURIER SYNTHESIS ONLY, WITH A PHE-COMPUTEO 3. 

IFSi-2 TO 00 FOURIER ANALYSIS ONLY# WITH A PRE-COMPUTEO S. 

IFEHR IS SET BY PROGRAM TO* 

=0 IF NO ERROR DETECTED. 

=1 IF M IS OuT OF RANGE., OR, WHEN IFS=+2,-2, THE 
PRE-COMPOTEO 3 TABLE IS NOT LARGE ENOUGH. 

=-l WHEN IPS =+1,-1, MEANS ONE IS RECOMPUTING S TABLE 
UNNECESSARILY. 

NOTE- AS STATED ABOVE, THE MAXIMUM VALUE OF V FOR THIS PROGRAM 
UN THE IBM 709m IS IS. FOR 360 MAC H I it S HAVING GREATER STORAGE 
CAPACITY, ONE MAY INCREASE THIS LIMIT BY REPLACING 13 IN 
STATEMENT 3 BELOW BY L0G2 N, WHERE N IS The MAX, NO. OF 
COMPLEX NUMBERS ONE CAN STORE IN HIGH-SPEtD CORE. ONE MUST 
ALSO AOU MORE DO STATEMENTS TO THE BINARY SORT ROUTINE 
FOLLOWING STATEMENT 24 AND CHANGE ThE EQUIVALtNCE STATEMENTS 
FUR THE K'S. 

DIMENSION A (4500) ,S( 1024) ,K(15) 

IF (M) 20,20,10 
10 IF (M-15) 40,40,20 
20 IFERR=1 
30 RETURN 
40 IF ERR = 0 
N = 2**M 

IF (lABS(IFS)-i) 440,440,50 
C WE ARE DOING TRANSFORM ONLY. SEE IF PRE-COMPUTEO 

C S TABLE IS SUFFICIENTLY LARGE 
50 IF (N-NP) 70,70,60 
60 IF ERR= 1 


GO TO i^ao 

C SCR A MOLE A / iJ Y SANOfS T HOD 

70 Ml )s2*N * 

DO 60 L-2/M 
60 K(L)sK(L-t)/2 
DO */ 0 L~M#19 
90 K(Lfl)s2 

C TMt FOLLOWING IS ST A1 EriENTS ARE TO COMPENSATE FUR A VJEAKHEoS IN 

C 1 HE FORTRAN V COMPILER * 

K1 =K(15) 

K2 =K(19) 

K 3 =K ( 1 3 ) 

K9 zK(12) 

KS "K ( 1 1 ) 

Kb ( 1 0 ) 

* 7 -K ( 9 ) 

Kb SKC6) 

K9 =K(7) 

K10=K(b) 

Ft 1 - K ( 5 ) 

K12=K(9) 

K l 3=K ( 3 ) 

K19=K(2) 

K1S=K(1) 

N2 =K(1) 

C NOTE EQUIVALENCE OK KL AND KC19-L) 

C til NARY 30RT- 

IJ = 2 
J1 =2 
110 J2 ~J1 
120 J 3 = J2 
130 J9 sJi 1 
190 JS = J9 
150 Jb = J5 
lbO J7 = J6 
170 J 6 = J 7 
160 J 9 = Jb 
190 J10=J9 
200 J l l = J l 0 
210 J12 = JU 
220 J 1 3 = J 1 2 
230 J19=J13 
290 JI=J19 

250 IF CIJ-JI) 260,270/270 
260 T=A(IJ-1 ) 

A(IJ-1)=A(JI-1) 

ACJI-l )*T 
T-ACIJ) 

A ( X J)-A(JI) 

A(jn = r 
270 I J- 1 J+2 
JI=J1+K19 

IF CJI.LE.K15) GO TO 250 
J19=J19+K13 

IF (J19.LE.K19) GO TO 290 
J13=J13tKl2 

IF (J13.LE.K13) GO TO 230 
Jl2 = J12tKl 1 
IF (J12.LE.K12) G 
JllsJil+KlO 
IF (Jll.LE.Kli) G 
J 1 0 = J 1 0 + K9 

IF (J10.LE.K10) GO TO 200 


0 TO 220 
0 TO 210 


J9 sJ'HKd 

It (J9.Lt.K4J GO TO 190 
Jb=J0+K7 

It (J6.1.E.K3) GO TO t«0 

J 7 " J I tf.it 

It ( J V • L E « K / ) GO TO 170 
Jfr=J6+K5 

It (Jb.Lt.K6) GO TO 160 
J5=J5+K4 

It (J5.LE.K5) GO TO 150 
Jt|rj4 + K3 

it (jM.Lt.K4) GO TO 1«0 
J3=J3+K2 

It (J3.LE.K3) GO TO 130 
J 2 = J 2 + K 1 

It (J2.Lt.K2) GO TO 120 
J 1 =J 1 +2 

It (Ji.Lt.Ki) GO TO 110 
IF (IKS) 280 , 20 , 300 

C DOING fOUKIEP A'i ALYS I S , SO OIV. BY N AND CONJUGATE . 

280 FN = N 

DO 290 1 = 1, N 
A( 2*1-1 )=A(2*I-1 ) 

290 A(2*I)=-A(2*I) 

C SPECIAL CASE- L= 1 

300 00 310 1=1, N, 2 
T = A(2*l-1) 

A(2*I-t) =T t A(2*I+1) 

A(2*I+t)=T-A(2*I+l) 

T = A (2* I ) 

a ( 2 * i y ~ r + a ( 2 * 1 + 2 ) 

310 A(2*I+2)= T - A(2*I+2) 

If (m-1) 20,30,320 
C SET FON L=2 

320 LE XP l =2 
C LtXPl=2**(L-l) 

LEXP=8 

C LEXP=2**(L+1 ) 

NPL = 2**MT 

C NPL = NP* 2**-L 

DO 390 L = 2,M 
C SPECIAL CASE- J=0 

CO 3M0 I=2,N2,LEXP 
11=1 + LEXP1 
12=11+ LEXP1 
13 =I2+LEXP1 
T = A ( I - 1 ) 

A ( I - 1 ) = T + A (12-1) 

A(I2-1) = T-ACI2-1) 

T = A ( I ) 

A ( i ) = T + A(I2) 

A ( 1 2 ) = T - A ( 1 2 ) 

T = - A ( 1 3 ) 

T 1 = A ( 1 3 — 1 ) 

A ( 1 3- 1 ) = A ( 1 1 - 1 ) - T 
A (13 ) = A ( C 1 ) - TI 

A ( 1 1 - 1 ) = A ( 1 1 - 1 ) + T 
340 A ( 1 1 ) = A( 11 ) +TI 

It ( L-2 ) 3«0, 300, 350 
■ 350 KLAS T =N2-LEXP 
J J =\PL 

DO 370 J= 9, LEX PI, 2 
NP J J =N 1 - J J 


uxsb(jj) 

1lA3T=J+kLAST 

DO 360 I=J,II.A3T/LEXP 

il=I<LEXPl 

I2aU+LEXPi 

13=I2+LEXPl 

T=A(I2-1 )*UK-A( [ d ) I 

TIsA(I2-l)*UI+A(I2)*UR 
A(I2-l)8A(I-l)-T 
A ( 1 a ) =A ( x ) - T I 
A ( I - 1 ) sA(I-i)+r 
A ( I ) sA(I)+TI 
T=-A(I3-i)*UI-A(I3)*UR 
T I = A ( X 3 — 1 ) * U R • A ( 13) * U I 
A(13-i)sA(Il-l)-T 
A ( I 3 ) =A(I1 )-TI 

A(Il-l)=A(Il-l)+T 
360 A ( I 1 ) sA(Xi) *TI 

C END UF I LOOP 

3/0 JJ=JJ+NPL 
C END UF J LOOP 

360 LEXP1=2*LK/P1 
LtXP = 2 * L £ X P 
590 NPLsNPL/2 
C END OF L LOOP 

IF (IPS) 410/20/30 

C DOING FOUKIER ANALYSIS. REPLACE A BY CONJUGATE. 

410 DO 420 1 = 1, N 
420 A(2*I) =-A(2*I) 

GO TO ‘30 
C RETURN 

C MAKE TABLE OF S(J)=SIN(2*PI*J/NP) /J=l/2/ . . . .NT-l/NT-NP/a 

440 N P - N 
MPsM 
NT=N/4 
M T = M - 2 

IF (Mi) 510/510/450 
450 THETAS. 7853961634 
C TMETA = PI/2**CLH) FOR L=1 

JSTEP = NT 

C JSTEP = 2* * ( MT-L + 1 ) FOR L = 1 

JD1F = NT/2 

C JDIF = 2**(MT- L) FUR L=1 

S(JDIF) = SIN (THETA) 

IF (MT-2) 510/470/470 
470 DO 500 L-2/MT 

ThETA = THETA/2, 

JSTEP2 r. JSTEP 
JSTEP a JDIF 
JDIF = JDIF/2 
SXJDIF )=SIN(THETA) 

JC 1 =NT -JDIF 
S( JC1)=C0S(THETA) 

JLAST=NT-JST£P2 
IF (JLAST-JSTEP) 500/480/4Q0 
460 DU 490 JaJSTtPz J-ASTz JSTEP 
J L = N T ** J 
JUsJ t JDIF 

490 5(JD)«3(J)*S(JCi)+S(JDIF)*S(JC) 

500 C U iM T I HUE 
510 IP (IPS) 70/30/70 
END 
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SUBROUTINE P309A7b (*IStCS) 

C 

C SUBROUTINE P30SAVE OUTPUTS TO T*i" THt RA’.v FSO (102a REAL NUMhtR$ PER 

C 2040 ftUHO RECORD) PRECEDED ti V A 25 «r«OHO ME ADtR RECORD. 

C 

C CALLIN'- ARGO *i.* 4 r : N3PC3 - IE. -JOT II OP HM: SERIES DA \ <' . [ U* 

C 

DIMENSION IARRAYC2048) 

COMMON UATUMSUSOO) , IHEAi;(20) rCAL(2r 4) , IVAHS/. ( 10) 

COMMON/ 1 GO/ MUSCLE(7#4) , IFUHRC2) 

ECU I VALENCE ( I A R R A Y rOATUMS(21QO) ) 

MTO = 2 

ltOh=oOOOOK 
!KlTti = 50000K 
ibfile=400ook 
c 

C CREATE AN OUTPUT HEADER RECORD 

C 

LO 10 Jsi,25 
1 A H R A Y ( J ) s 0 
10 CONTINUE 

OU 20 J=l/20 
lAKRAY(J)=XHEAO(J) 

20 con r anus 

NRECSsNSECS/2 
IF (NRECS.fcQ.O) n R E C S = 1 
IARRAY(21)=NR6CS 
Iy.HDs.IRI TEt25 

CALL MTDIO (MTO, InRD, IARRAY , ISTAT , IER) 

C 

C OUTPUT FOLLOWING DATA RECORDS CONTAINING TrlE R Art PSD VALUES 

C * ‘ 

K-0 

DO 50 J=1,NREC3 
DO 40 JJnt , 1024 
D A TUNS (20R9*JJ) = U A TUMS ( K t J J ) 

40 CONTINUE 

XinRD = IRX TE + 2048 

CALL MTDXO ( M f 0 , X rtRD , X ARR A Y , I ST A T , I ER ) 

K=K+1024 
50 CONTINUE 

V.RITE TrtO EOFS AND BACK UP TO JUST BEFORE THEM 


DO 60 J=1 ,2 

60 CALL MTDIO (MTO, IEOF, I ARRAY , ISTAT , IER ) 

X BACKS I bF X LE + 1 
DU 70 J= I # 2 

70 CALL MTDIO ( M TO , I B ACK , I ARR A Y , I ST A T / I ER ) 
RE fURN 
tNl) 
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OVERLAY 0V7 

SUBHUul IMt Hour (DrYPE,STAKT,NScCSf SPAN, VULTSEC , PSU5UM , $ TUERR , toMOA ) 

SUBROUTINE HUH OUTPUTS TO THE LINE PRINTER HEADER AND CALIBRATION DA f 
AVERAGED FORCE DATA, AND SMOOTHED PSD SPECTRUM. 

CALLING ARGUMENTS: DTYHE - TYPE UF DA 1 A 

START - BEGIN TIME OF PRESENT DATA SLICE 

NSECS ~ LENGTH OF DATA SLICE 

SPAN - FOR FORCE DATA NUMBER OF SECS FORCE DAI 
AVERAGED OVER; FOR PSD DATA FREQUENCY 
INTERVAL FOR SMOOTHING ESTIMATES 

VULTSEC- INTEGRATED tMG VALUE (TIME DOMAIN) 

PSDSUM- INTEGRATED ENG VALUE (FREQUENCY DOMAIN) 

STDERR- NORMALIZED STANDARD ERROR OF ESTIMATE F< 
PSD VALUES 

H 0 A - FLAG SET a h £ N DATA PROCESSING ATTEMPTED 
to I THOU T CALIBRATION 


INTEGER DTYPErtfHOA 

DIMENSION I EDA TE (3) , I DD A TE ( 3 ) , CUN PC T ( a 1 0 ) 

DIMENSION AvGLBB(iSO) , PSD (4 10) / ITIME(3) ,PSUh(410) ,PERCNT(410) 
COMMON DATUMS (400 0) , IhEAD(20) , C AL ( 2 , a ) , I V ARSW ( 1 0 ) 

COMMON/ I GO/ MUSCLE (7, a) , IF0RR(2) 


EQUIVALENCE 

equivalence 

equivalence 

equivalence 

equivalence 

equivalence 

equivalence 

EQUIVALENCE 
EQUIVALENCE 
DATA MUSCLE 

a)', ' AL # , ' 1 3 # 
a; , # S ' , ' ' , 

DATA .IFURR 


(IS'JSNO, I HE AD ( 1 ) ) , ( ICOATt , I HE AD ( 2 ) ) 

(IDOATE, I HE AD (5)) , (IFLITE, IHEAD(6) ) 

( ..RUN, I HEAD (9) ) , (MJS, (HEAD (20 ) ) 

(i A TAPE, IHfeflD(lO) ) / ( I SAMP E , IHEAD ( 1 1 ) ) 
CiSAMPM, IHEAOU2) ) , (AVGLBS,DAIU»S(4000)) 
(PSD, DATUMS ( 1 ) ) , ( P S D N , 0 A FUMS(4000) ) 

( P E R C N T , DATUM S(3000) ) , (CUMPCT , DATUMS (3500) ) 
(EXP VAL, DATUMS (2500) ), (STOEV# DATUMS (250 1 ) ) 
(EMGVAR f DATUMS ( 4500 ) ) / ( I T I v E , I ' E A D ( 1 b ) ) 

/ # BR ' , ' AC ' , 'HI', ' AL ' > ' b', 'ICS 'EP', ' 6','. 
' ' , ' G ' , 'AS', ' T R ' , ' 0 C ' , ' N t ' , 'MI', 'US', ' S 

' '/ 

' F ' , 'R+ ' / 


LP T = 1 2 

IF(AHOA.EQ.l) GO TO ISO 


'RA 
OL 9 


, 9 D I 
' EU ' 


BRANCH TO THE PROPER PRINT OUT SECTION FOR THE DATA TYPE 


IF(DTYPE-4) 3, SO, 100 
PRINT OUT HEADER INFORMATION 


S IP = 1 

IF (IFLITE.GE.O) I P = 2 
to R I T E ( L P T , 10) 

10 FORMAT (1H1,//,48X, 'CARDIOVASCULAR LABORATORY ',//, «7X , 'EMG DATA PRO 

oJCESSING PROGRAM ' ,/,U?X, ' ') 

wR I TE (LPT, 20) I 3UUN0, 1EDATE, IFORR( IP) , IFL I TE , I RUN , I A T APE , I DD AT E , (MU 
a> S C L E ( J , M u 3 ) , J •= i * / ) / ISAMPE, ISAMPM 

20 FONMAT(////#47X, '*** HEADER INFORMATION***',///, 9 SUBJECT NO.:', 7 
aJX , 13, 1 8X , 'EXPERIMENT D A TE : ' , 6X , 1 2 f ' / M 2 , ' / ' , I 2,1 0 X , ' F L I GH T RErERE 


o r n o o o 


if.Lt DAY: ', A3, 13,//, ' RUN NO. :’, 1 U, 12. 1>X, 'ANALOG TAPE NO.:', / a , i 
33, lax, 'DIGITIZING DATE: ',bX,I2, '/', 12, '/', 12,//, ' MUSCLES', 9 X,/A2 
«,//,' KMf. 3AMPL& HATE C S AkP/SEC ) : ' , 1 5 , 8X , ' F ORCE SAMPLE RATE (SAMP/ 
aise.o : ',ibj 

AH I It (LPT ,3.1) (CAl.(l,.l),J"l,'l) 

30 FO*>! .11 (////, 4 lit •** C AH UP A T (ON DATA *>',///,' FORCE CAL',/,' 

at"””” ■•',//,' CAL N E I G h T (Lb3)t , ,F/.l,loX,'C0UhTS/Lt>:'>2A>F6.0,bX,'A 
«VG. OA3EL1NE COU M T I ' , 2X , F 6 . 0, bX , 'CAL GAIN: ' , 2 X , F 3 . U ) 

WRITEUP T, 40) (CAL(2,J),J=1,2) 

40 FOR MAT (///,' F.MG CAL',/,' - ',//,' RMS AMPLITUDE (MICROVOLTS) 

<0: ' ,F6.0,7X, ' COd. V T S/MICRO VOLT S',F 5. 2,//) 

GO TO 200 

PRINT OUT FORCE DATA 
bO WRITEUP T, 60) I T I ME , MSECS, SPAN 

60 FOKMAT(lMl,blX, '*** FORCE DATA ***'///5X, 'START T I ME- ' , 1 3 , ' : ' , I 2 , ' 
a)! ' , 12, lbX, 'DATA LtNGTH (SEC3) 14 , 1 5X ,' AVERAGE FORCE INTERVAL (SE 
3CS)-',F3.0,//) 

LOOPSsNSECS/IFIX (SPAN) 

IF (LUUPS.GT.H9) GO TO 66 
WK I TE (LP T , 62 ) 

62 FORMAK47X, 'INTERVAL', 5X, 'AVERAGE FORCE U 33 ) ' / , 47 X , ' ',5X, 

a), - - -•-',//) 

DO o5 J=l, LOOPS 

wR I TE (LPT , 64 ) J , AVGLOS ( J ) 

64 FORMAT ( 49X , Is, 1 5X , Fb . 1) 

6b CONTINUE 
GO TO 200 

66 IF ((L00PS/2W.NE. LOOPS) LOOPS = LOOPS- 1 
WRITE (LPT ,68) 

68 FOrMA I ( 1 bX , ’ In T ER V AL ' , 5X , 'AVERAGE FORCE (LOS) ' ,20X, 'INTERVAL' ,bX, ' 

a) A VEKAGt FORCE (L03) ’,/, 15X, ' ',5X,' ',20X 

3, ' , b X , ' ',//) 

L00PS2= LOO PS/2 
00 7b J= 1 , LUOPS2 
' JJ-J+LU0PS2 

WRITEUP F, 73) J , A VGL6S ( J ) , J J, AVGLbS ( J J ) 

73 PORMAT(t7X,I3,lSX,F5.1,29X,I3,15X,F5.1) 

7b CONTINUE 
GO TO 200 

PRINT OUT PSD DATA 

100 WRITE (LPT, 110) ITIME,N SECS, VOLTSEC, SPAN, STDERR,PSDSUM,EXPVAL,STDEV 
3, EMGVAR 

110 FORMAT ( 1 H 1 , 3 9 X , ' * * * POWER SPECTRAL DENSITY OF EMG DATA ***'///4X,' 
3 3 T A R Y TIME-' , 13, * : ',12, ' : ' , 12, 1 IX, 'DATA LENGTH (SECS)-', 14, 14a, 'IN 
JlTtGRATED EMG ( MICROVOLT* SEC)-', E12.4//4X, 'SAND WIDTH ( hZ ) - ' , F 6 . 3 , l 0 
3X , 'NORMALIZED STANDARD ERH0 k-'»F 6.3,5X* 'INTEGRATED PSD (MICROVOLT* 
3*2) - ' ,E13.4,//,4X, 'MEAN ( HZ ) - ' , F6 . 1 , 1 bX , ' ST AND ARD DEVIATION 

3(riZ)~',F6. 1,8:, 'tMG VARIANCE (MI CROVOL T * *2 ) - ' , E 1 5 . 4 , / ) 

LOOPS=IFIX (400./SPAN)+l 
FREU=SPAN/2. 

IF (LOOPS. GT. 46) GO TO 130 
WRITEUPT, lib) 

lib FORMAT (30X , ' FREQ ' , 1 3X, 'PSO ' , 1 i X , 'PSD ' , 1 OX, 'X OF',lOX,'CUM X'/,30X 
3, ' (HZ) ' ,8X, ' (M •tV**2/HZ) ' ,6X, 'NORM' ,9X, 'TOT AL',9X, 'TOTAL',/, 31X,' 

3 ' , 6X , ' ' ,8A, ' ' ,9X, 9X, '-----' ,/) 

DO 125 Jsl , LOOPS 

WRITE (LPT, 120) FREu,PSD( J) ,PSDN(J) ,PF.RCNT(J) ,CUMPC T (J) 

120 FORMAT (30X»F6.2, /X,F10.3,8X»F6.4»7X,F6.2»6X,F6.2) 

FREu=FREU+SPAN 



non 


1 * ■ j I i .1 — - - » mm i i mmm 

l<ii COTitINUE 

to io 200 

ISO It- C (L00P3/2 ) *2 .NE . LOOPS ) LOOPS -LOOPS- 1 
FREU2=FLUAT (L£0PS/2)*SPAN+FH£0 
MUTE (LPT. 135) 

IV.) »-UKiv,uf(' KPEO', OX,', 'SO', VX , 'PSD , ‘ X OF ' , OX , 'CL' ' ’< * , 22X , 'I Ht'! 

<!)', VX, 'PSD' » OX, 'PoD'/Sa, V. OK', OX, 'CUM '/, ' / 2X , ' (HZ ) ' , '4 X , ' (M V >*2/ 
alMZ ) ' , 5X > 'nORM'./X, 'TUTAL'.SX, 'TOTAL ',22X,'(h2)’,5X, '(MMV**2/HZ) ',5 

A»X, 'N0WM',6X, 'TOTAL ',bX, 'TOTAL',/, ','ix, '- ',5X, * 

- ' , 7X , ’ ',5X, 22X, ' ' , 5X , ' ',5X, 6X, 

a)'- -',6X, 

LOUPS2=LUOPS/2 
DO l«a J" 1 , LU0PS2 
J J = J +LUQPS2 

OKU l (LPT ,1*40) FREG,P8D<J) ,PSDfUJ) ,PERCNT( J) , CUSPC T ( J ) , FREQ2, P30 ( J 
a' J > PSUij ( J J ) , PERCH f (J J) , CUN'PCT ( JJ ) 

140 FORK. AT (F7.2,3X,F10.i,5X,Fb.«,5X,F6.2,SX,F6.2,20X,F6.2,4X,F10.3,SX, 
&F6.4»5X,F6.2»5X,F6.2) 

FREO=FREU+3PAN 
FREG2=FNEQ2+SPAN 
145 CONTINUE 
00 TO 200 

ERROR MESSAGE 

150 WHITE (LPT, 160) 

160 FORMAT!' ERROR. , .DATA ANALYSIS ATTEMPTED tilThOUT CALIBRATION FIRST 
in. CHECK DATA CARD ORDER.') 

S TOP 

200 RETURN 
END 
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4.2.0 LBNP ONLINE WITH SIGMA 3 

REQUIREMENT DELETED 

4. 3.0.0 General Data Acquisition System. 

The DBAS (General Data Acquisition System) acquires, digitizes, 
and records test data from an analog playback. The user is allowed up to 
3 different rates and sixteen channels of input. The only restriction 
are that the rates must be integer multiples of each other and that the 
A/D channels be grouped in the order of highest rate, mid-rate, and slow 
rate. The GDAS is loaded from the nine-track tape and will operate 
independently of theSigma 3 Biomedical System. It will, however, require 
the FM recorder and Systron Donner used by the Sigma 3 unless the user can 
supply his own recorder and Systron Donner. The GDAS data storage and 
output, operator interaction, and subroutines are described in the follow- 
ing sections. 


4.3. 1 .0 Data Storage 

There ar: two types of data storage - that which is used by 
the executive for control purposes; and that which is used as buffers 
for digitized data accumulation. The data used by the executive for 
control is stored in locations 80 i 6 to A5 j 6 of scratch pad (locations 
0-255io). The following table defines these data. 


LOCATION 

VARIABLE 

80 , s (128 i„) 

CLCNT 

81 , 6 ( 1 29 1 o ) 

NDRATE 

82 j 6 (130 io ) 

N0L0 

83 ,6 (131, o) 

L0CNT 

84 i6 ( 1 32 i o ) 

N0MED 

85 i6(133 l0 ) 

MEDCNT 

86,6 (134, o) 

N0HS 


DESCRIPTION 

User input number of microseconds to 
elapse between each A/D read of the 
highest speed data. 

User input number of different rates 
(1. 2, or 3) 

Minus total number of channels (com- 
puted from user input) 

Number of reads of the high speed 
channels before the lowest speed chan- 
nel is read. 

User input minus total number of medium 
and high speed channels. 

Number of reads of the high speed chan- 
nels before the medium speed channels 
are read (user input) 

Number of high speed channels to read - 
user input 


H-l 


LOCATION 

VARIABLE 

DESCRIPTION 

87,6 

(135,0 ) 

BUFZER 

Buffer zero starting address minus 1. 
Fixed at 77Cie 

88,6 

(136, 0 ) 

BUFZND 

Buffer zero end address. 

Set to B64ie, but will be changed if 
the user sets the data buffer length 
to less than 1000. 

CXI 

lO 

<T> 

(137, 0 ) 

BUF0NE 

Buffer one starting address minus 1. 
Fixed at B7Dj 6 - 

« 

< 

CO 

(138, ,) 

NDBUF1 

End address for Buffer one. 

Set at F65 16 , but will be changed if 
the user selects a data buffer length 
less than 1000. 

8B, 6 

(139, 0 ) 

BUFUSE 

Buffer in use pointer. 

Initially set at zero but is changed 
when data buffer 1 is being filled 
with digitized data. 

8C , 6 

(140, 0 ) 

MAXBUF 

Maximum buffer length including header 
information. Fixed at 401 1 6 ( 1 025 i 0 ) 

8D, 6 

(141, o) 

TP0UT 

Tape output flag Initially set to 
zero, but is set to • when a data buf- 
fer has been filled ano is ready for 
output to tape. 

00 

m 

a> 

(142, 0 ) 

D0NCLK 

Flag, initially zero, set to one when- 
ever the Systron Donner clock timer 
has been read and are ready to be 
unpacked. 

8F,6 

(143, 0 ) 

BGNCHN 

Address of the first high speed A/D 
convertor channel. Default value is 
zero. 

90,6 

o 

•3- 

ADCRD 

Flag which indicates an A/D convertor 
read is in progress (=1) 

91,6 

(145, 0 ) 

D0NMS 

Timer read from the Systron 

92,6 

(146, 0 ) 

DNSEL 

Donner when data buffers are 

93 16 

(147, 0 ) 

DNMIN 

switched. These times are in a 

94,6 

(148, „) 

DN0UR 

pseudo BCD format, Complemented. 

95,6 

(149, 0 ) 

STHR 

| Systron Donner timer at which 

96,6 

(150, „) 

STM IN 1 

the user wants to begin digitizing 

97,6 

(151, 0 ) 

STSEL 1 

data. Default values are zero. 

98,6 

(152, 0 ) 

STRMS 
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LOCATION 

VARIABLE 

DESCRIPTION 

99 16 
9A, 6 
9B 1 6 
9C,6 

(153. 0) 

(154. 0) 

(155. 0 ) 

( 1 56 . 0 ) 

STPHR 

STPMIN 

STPSEL 

ST0PMS 

Systron Donner timer at which the 
the user wants to stop digitizing 
data. Default values are 1 EOO 1 6 • 

9D U 

(157, 0 ) 

BUFHD1 

Buffer one header address minus 
one. Used for tape output. 

9E,6 

(158, 0 ) 

BUFHDO 

Buffer zero header address minus 
one. Used for tape output. 

90,6 

(159, 0 ) 

MBUFLG 

Negative of Maximum buffer length 
including header. Set by user with 
a default value of - 1 025 i 0 • 

AOi 6 

( 1 60 1 o ) 

MXDBUF 

Maximum data buffer length. 
Set at 1 OOOi o . 

A1 ,6 

(161, o) 

PAR 

Number of parity errors detected 
during tape output of the digitized 
data . 

A2, 6 

(162, 0 ) 

TIMG 

Number of timing errors detected 
during tape output of the digitized 
data . 

A3, e 

(163, „) 

HDBFLG 

Length of the header portion 
of each data buffer. Set at 25 10 . 

A4, 6 

( 1 64 , o ) 

BUFTTY 

Address of the teletype data buffer. 
Address is 100 16 (256 10 ). 

A5,s 

(165, 0 ) 

STATUS 

Temporary storage for the status word 
of the tape write operation. 

A6,6 

( 1 66, o ) 

TAPEND 

Flag signifying end of operation on 
tape. Writes during test. 

Also used for control purposes, are the addresses of subroutines which are 
general in nature or called by more than one other subroutine. These are 
defined below 

70 u 

(112,o ) 

FETCH 

Used to coordinate data input, e.g., 
reads a buffer from teletype and checks 
for data errors. 

71,6 

(113, ») 

CRLF 

Outputs a carriage return-line feed to 
to the teletype. 

72,6 

(114, o) 

0TL 

Outputs a line to the teletype. 
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73„(115„) 

UNPACK 

Unpacks the System Donner Time from 
pseudo BCD to binary. 

74 16 (116„) 

CjDNVR 

Converts ASCIi' Input number to a signed 
binary integer and checks for overflow 
(number too Urge). 

75 u (117 10 ) 

TDUMP 

Starts tape I/O during a test and writes 
E0F Marks after test. 

76„(118,o) 

0DEL 

Converts a binary integer to a string of 
ASCII digits for output. 

77 l6 (119io) 

ITT 

Inputs one ASCII Character from the 
teletype. 

78 1 * ( 1 20 1 » ) 

0TT 

Outputs one ASCII character to the 
teletype. 

79i 6 (121,o) 

TYPIN 

Routine called to fill the teletype 
input buffer. 

7A 16 (122 10 ) 

REST0R 

Restores the executive data table 
(locations 80 16 to A5 16 ) to it's 
initial value, and zero the digitized 
data buffer areas. 

7B i e ( 1 23 i o ) 

ST0PIT 

Routine to stop the data acquisition 
and restore the acquisition routine to 
its initial state. 

The buffers which are used to accumulate the digitized data consist of two 
blocks of core, each 1 025 i 0 words in length. The first 25 10 words of each 
block are leader data, 20 10 of which may be set by the user. The blocks 
are described below. 

LOCATION 

BUFFER WORD 

PURPOSE 

764 le \ 

B65 ! 6 / 

1 

Buffer record counter incremental by 
the executive. 

765-768 ! 6 \ 
B66-B69 i 6 f 

2 - 5 

Storage for Systron Donner time at 
which the first data for this record 
were acquired. 

7 69-77 C 1 6 1 

B6A-B7D j 6 f 

6-25 

Header information stored by the user. 
The 20 10 words are constant during a 


given acquisition. 
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77D-B64 16 * 26 Digitized data. Note that these end 

B7E-F65 16 / 1025 addresses are only valid for the max- 

imum data buffer length, 1000 10 . Also 
the data will be stored in the order 
that it is read, ie. there is no attempt 
to group the data by channels, rates, etc. 


4. 3. 2.0 Operator Interaction with the DGAS 

The GDAS acts in a conversational manner with the user via the 
teletype. In response to questions and for all data (numeric) input, a 
slash (/) is used to terminate the input. For correction of input prior 
to the slash, the last character input can be deleted by an up arrow (t) 
and the entire line deleted by a back («-) . Of the particular input is 
rejected by the GDAS, the back arrow and a carriage return - line feed will 
be printed. The quantity will be requested again. 

The following inputs are requested by the GDAS. 

Number of different data rates, the 
user must input, 1, 2, or 3. 


Number of microseconds to elapse between 
reads of the highest data rate. For 
example, 1000 samples/second would be 
1000 microseconds; 100 samples/second 
would be 10000 microseconds. The lowest 
rate which can be input for a high speed 
rate is 62500 microseconds/read or 16 
samples/second. This does not preclude 
lower rates if more than one rate is used. 


Length of the data buffer to be filled 
before output to tape. This number should 
be chosen with care-especially with multiple 
data rates-so that the data from each rate 
will appear at the same place in each buf- 
fer. For example, with four channels of 
data at 320 samples/second and eight chan- 
nels at 2 samples/second, a good number 
for buffer length would be 648. (4x160+8). 

The buffer would be dumped each H second. 

ADC CHANNEL Address of the first A/D convertor chan- 

nel to be read, i.e. the first high speed 
channel . 


0 UAb. 

No. Rates 


m 

|3) 


MICROSEC PER HI SPD READ. 


DATA BUF LENGTH, MAX I S 1000 
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GAIN 

NO. HI SPD CHANS. 

NO. MED SPED CHANS. 

NO. LO SPED CHANS. 

NO. HI SPD READS/MEDIUM READ. 

NO. HISPD READS/LOW READ. 
INPUT DATA FOR TAPE HEADER? 


WRITE BEGINNING E0F? 


Gain to be applied to each A/D 
conversion. Gains available are 
times 1, 2, 4, or 8. 

Number of high speed channels to read. 

For one rate, this number will be the 
total number of channels. 

Number of medium speed channels to read. 

The output applies only when there are 
two or three rates. When there are two 
rates, this is the number of lower speed 
channels to read. 

Number of low speed channels. This output 
applies only when there are three rates. 

This input is the number of channels for 
the lowest rate. 

Input the number of times to read the 
fastest rate before reading the mid-rate 
(or in the case of two rates, the low 
rate), e.g.. Suppose the high rate is 
100 samples/second and the mid rate is 
25 samples/second. The number to be 
input would be 4 = 100/25. 

In the case of three data rates, how many 
times are the highest speed channels read 
before the lowest speed channels are read? 
Eg. Suppose the rates are 500, 250, and 25 
samples per second, respectively. The 
response to this input would be 500/25 = 20. 

Yes or no response. This question asks if 
the user desires to store constant numeric 
data in words 6 through 25 in the header of 
each buffer. This data is output to tape 
with the digitized data. See 4. 3. 1.0 for 
further data description. If a yes response 
is given, the GDAS asks for the input to 
word 6, then 7, then eight, etc. If the 
user desires to terminate the header input, 
he may enter letter A followed by a slash. 

This allows the user to write a end of file 
on the beginning of his tape. One reason to 
do so is to allow the tape to space past the 
extra 3 inches at the beginning. This 3 
inch gap is especially important during high 
speed acquisition rates. (See 4. 3. 2.1 for 
more detail . ) 
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START TIMES? 


ST0P TIMES? 


Does the user want to input the Systron 
Donner times (hr, min. Sec) at which the 
GDAS is to start digitizing data? Yes 
or No response required! If yes, the 
user must input the times as requested 
If no, the GDAS uses zero as the default 
value. 

Same as start times except these times 
are when to stop acquiring data. The 
default value is lEOOie if a no response 
is entered. 


A > 

S TO START This is an input which prevents the GDAS 

from looking at the start times. It is 
useful if the default is entered or if 
the user must position an analog tape for 
playback. The quantity S/ will start the 
GDAS. 

When a test is completed, either by reaching the stop times or by an x input, 
the GDAS will allow additional tests to be conducted with the same data 
input. This question is 


CONTINUE TEST WITH NEW TIMES? 

If a yes response is entered, then new start/stop times are requested. Then 
the test may be restarted by using the s/input. If a no response is entered, 
the GDAS will ask of the tape is to be rewound between tests. If yes, it 
rewinds the tapes. In either case, the GDAS initializes itself, i.e., the 
user must enter new inputs for rates, etc. 


4. 3. 2.1. Timing Estimates for GDAS 

The GDAS will digitize data at various rates, write various 
length records depending on the user's inputs. For this reason, the user 
should put some "handle" on whether his specific rates and data buffer 
lengths are feasibile. As mentioned previously, the minimum single rate 
speed is 16 samples/second. To determine the maximum rate, several tests 
were conducted. The results of these tests indicate that the tape speed 
will be the limiting factor. The only control the user has over tape 
speed is to vary his buffer length. A longer buffer will take longer to 
fill than to dump in most cases. The limiting case is: 

Time to fill buffer 
with data in seconds 
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0.046 + 0.0001 (buffer 
length in words). 



The above inequality has been tested with 6 channels of data at 1000 
samples/second using various buffer lengths. The minimum buffer length 
which worked was 690. The time required to fill the buffer and to dump 
the buffer in this case is identical. Note that the above inequality and 
test are based on the GDAS not being required to write a record from the 
BOT (beinning of tape) position. The first operation from BOT requires 
.197 seconds to simply start and stop the tape. Hence, if the user is 
dealing with high rates, he must not start the tape from BOT. A beginning 
dunmy record or an end of file can be used to position past the BOT. 
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4. 3. 3. 0.0 

All routines 

.0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

4. 3. 3. 1.0 

4. 3. 3. 1.1 


4. 3. 3. 1.2 


4. 3. 3. 1.3 


REST0P 


CRLF 


0TL 

FETCH 


TYPIN 


CHKIT 


Subroutine Description 

are documented according to the pattern shown below: 

name 

purpose 

calling sequence 

software/software interfaces 

input data 

output data 

storage required 

description 

flow diagrams 

BEGIN - Pre- and Post Test Control Routine 

The purpose of this program is to control the flow 
of input data from the user and to insure that all 
required inputs are within tolerances. BEGIN also 
starts the test by calling the test routine. 

Calling Sequence - 
JST *BEGIN 

Where BEGIN is a core location which contains the 
value 1 lOie . 

Software/Software Interfaces 

BEGIN uses both external and internal subroutines. 

If no indication is given, the subroutine should 
be assumed to be external . 

Subroutine called to restores all the GDAS variables 
that may be changed by the user to their default 
values. It also resets all the 1/0 devices. 

Subroutine called to issue a carriage return - line 
feed to the teletype. 

Subroutine called to output a message to the teletype. 

Subroutine called to output a message and return the 
numeric response in the A register. It skips forward 
one instruction if the numeric input was free from 
character errors. 

Subroutine called to fill the 16 character teletype 
input buffer, BUFR. 

Internal subroutine used to check for the no t hang*? 
condition in variables whose input is optional. The 
routine FETCH returns an error condition for this 
particular case when possibly no input error has been 
made. 
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E INS - 
ZWEI - 
DREI - 
HEADER 


TPCHK 
ADCHK - 
READ0N - 

UNPACK - 

TEST - 

4. 3. 3. 1.4 

BUFR 

N0RATE 

CLCNT 

MAXBUF 

BUFZER 

BUFZND 

BUF0NE 

flDBUFl 


Subroutine called to input data pertinent to a one 
rate data acquisition. 

Subroutine called to input data pertinent to a two 
rate data acquisition. 

Subroutine called to input data pertinent to a three 
rate data acquisition 

Subroutine called to input and store constant data 
into the 25 word tape buffer header. Only twenty 
words of input are allowed: The other five are 

reserved for record number and time. 

Subroutine called to check for the tape ready condi- 
tion, i.e. tape is on-line and not write protected. 

Subroutine called to check that the A/D convertor is 
not in mannual mode. 

Internal subroutine called to read the Systron Donner 
time code generator to determine when to start the 
test. 

Subroutine called to unpack the Systron Donner time 
from pseudo BCD to binary. 

Subroutine called to start the data acquisition and 
remain in control throughout the test. 

Input Data - The following parameters must be available 
for storage/use by BEGIN 

A 16 word tel type input buffer filled by TYPIN. 

Number of A/D rates 

Number of microseconds between highest rate reads 
of the A/D convertor 

Maximum data buffer length (default lOOOio) 

Starting address minus 1 for data buffer zero. (Set 
at 77C 1 6 ) 

End address of buffer zero (default B64 ig) 

Start address minus 1 for Buffer One (set to B7Dic) 

End address of buffer one (default F65 1 g ) 
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HDBFLG 

MBUFLG 

BGNCHN 

STHR 
STM IN 
STSEC 
STRMS 

j 

STPHR ) 

STPMIN / 

STPSEL 
ST0PMS 

In addition to the above values, the pointers to subroutines FETCH, 
CRLF, REST0R, UNPACK, and TYPIN must be available in Scratch pad 
(locations 0 to 255). BUFR is stored at 256 i o to 271 lfi . The binary 
numbers for the above values are returned from subroutine FETCH. 

The inputs which may be alphabetical are put in the teletype buffer, 
BUFR, and then checked for the appropriate alphabetical vaiue by 
BEGIN. 

The hours, minute, second, and millisecond are input by internal 
subroutine READ0N. 


Header length on data buffer (set at 25) 

Minus buffer length including header (defaul t=-1025i 0 ) 
First A/D channel to read (default = 0) 

Systron Donner Strat hour, minute, second, and 
Millisecond (default is 0) 


Systron Donner stop hour, minute, second, and 
millisecond (default is 1E00 1G ) 


4. 3.3. 1 .5 Output Data 


The specific names in the input section may be changed if there is a 
default value. In addition to the above outputs, the following sub- 
routines have calling arguments as indicated. 


FETCH - 
0TL 


READ0N - 
UNPACK - 
CHKIT - 


The address of the message requesting numeric input 
is transfered through the x-register. 

The terminating character of a message is transfered 
through the S-register. The address pointing to the 
first word of the message is stored in the first 
location after the call to 0TL. 

The code for the time to be returned (hours, minutes, 
seconds, milliseconds) is passed through the x register. 

The pseudo BCD time code in passed to this subroutine 
through the A register. 

The first character of the teletype buffer, BUFR, 
must be available. Since CHKIT is an internal sub- 
routine within +255 of BUFR, no special action is 
required. 
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Any single character to be output before filling the 
teletype buffer with input must be passed through 
the A register. 


Storage Required - BEGIN requires 183jg or 387jg 
words of memory. 


Description 

The tasks of BEGIN are to: obtain all the required 
and optional inputs from the user; insure that the tape and A/D 
converter are ready; and to start the acquisition at the specified 
Systron Donner time (If the default values are used, the acquisition 
starts immediately) . BEGIN Accomplish. The first task by calling 
the appropriate support subroutines to get a message to the user 
requesting input and return a binary number for those numeric inputs. 
In the case of a question, a yes or no response is required. If an 
error is detected in either the optional or mandatory input, the 
request for data is typed again. In the specific case of the number 
of microseconds per high speed read, BEGIN cannot determine if the 
error was legitimate. It requests that the user verify that the 
input is okay. The only correct input which can trigger the veri- 
fication is if the number of microseconds is greater than 32767 and 
less than 65535. 


T Y P I N 

4 . 3 - 3 . 1 . 6 . 
4 . 3 . 3 . 1 .7 


The second task is accomplished by calling the A/D 
convertor and tape check routines. The A/D convertor checks only 
for the convertor being in manual mode. If so, the routine outputs 
diagnostic message and continues checking until the convertor is 
placed in automatic mode. The tape check routine checks for tape 
on-line and not write protected. If either/or both conditions are 
not true, a diagnostic is output. The condition(s) must be connected 
before BEGIN will continue. 

After the tape and A/D convertor are checked, BEGIN 
will then start the test sequence by asking about writing a beginning 
end of file. The user can write an end of file to position past the 
first three inches of tape (recommended for high data rates). BEGIN 
will then request an S to begin checking the Systron Donner times 
against the start times. This allows the user a chance to review 
his input data before actually starting the test. 

After the test acquisition is completed, the user 
is asked about continuing the acquisition with new times. This allows 
the user to have a short calibration run and then a long test without 
changing his fixed data. The tape can be rewound if the user desires. 
The GDAS then initializes itself and starts the whole sequence again. 
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4. 3. 3. 2.0 Subroutine FETCH 

4. 3. 3. 2.1 The purpose of this subroutine is to output a 
message requesting numeric input and convert the 
subsequent input to a binary integer. 

4. 3. 3. 2. 2 The calling sequence is LDX Address of message 
JST FETCH 

Returns here if input in error 
Returns here if input is okay 
Binary integer is in the A register. 


4. 3. 3. 2. 3 Software/Software Interfaces. 

FETCH Interfaces with the following subroutines - 

CRLF Subroutine which issues a carriage return - line 

feed to the teletype 

0TL Subroutine which writes a line to the teletype 

TYPIN Subroutine which fills the teletype input buffer. 

CONVR Subroutine v/hich converts ASLII integers to a 

binary number. 


4.3. 3.2.4 Input Data - The address of the message to be 

output is transfered to FETCH through the x register. 
A binary number is returned from C0NVR in the A 
register, and if an error was detected, the over- 
flow bit will be set. The address of the teletype 
buffer must be available in location A4 16 of scratch- 
pad. 


4. 3. 3. 2. 5 Output Data - Subroutine FETCH passes the terminating 
character of the message requesting input through 
the A register. The address of the first word of the 
message is stored in the first location after the 
call to the message writer, 0TL. 

FETCH transfers a blank character to the teletype 
input routine, TYPIN, through the A register. 

The address of the teletype buffer is transfered to 
the ASLII to binary conversion routine, C0NVR, through 
the x register. 

The binary number that FETCH was requested to input 
is transfered to the cal ling program through the A 
regi ster . 

4. 3. 3. 2. 6 Storage Required - Subroutine FETCH requires D^g or 
13 iQ locations. 
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4. 3. 3. 2. 7 


Description 


When FETCH is called, the address of the message 
requesting a*i input number is in the x register. FETCH stores that 
address in ;he first location after the call to the message writer, 
0TL. FETCH then issues a carriage return-line feed by calling the 
routine, CRLF. After returning from CRLF, FETCH loads the terminating 
character into the A register before calling, 0TL. The terminating 
character is always a period (.) for any message from FETCH. After 
the message is printed by 0TL, FETCH calls the routine TYPIN to fill 
the teletype input buffer. Upon return from TYPIN, FETCH loads the 
address of the teletype buffer into the x register and calls the 
conversion routine, CONVR. When CONVR returns control to FETCH, a 
binary number will be in the A register. The overflow bit will be 
set (i.e., 1) if CONVR was unable to convert the number. If the 
overflow is set, FETCH returns control to the first location after 
the calling program. If the overflow is reset (i.e., 0) the return 
is to the second location after the call. 
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4. 3. 3. 3.0 

Subroutine CL0CK and STOPIT 

4. 3. 3. 3.1 

The purpose of this subroutine is to respond to the 
real time clock interrupt and to start the data 
acquisition. The purpose STOPIT is to reset the 
clock subroutine. 

4. 3. 3. 3. 2 

There is no calling sequence since this is an 
interrupt subroutine. 

4. 3. 3. 3. 3 

Software/Software Interfaces. 


CLOCK calls internal subroutine. 

STRTAD to reset and start the A/D Convertor. 
Internal subroutine ST0PIT is called by external 
programs to reset the CL0CK routine variables to 
their initial state. 

4. 3. 3. 3. 4 

Input Data-There are no data inputs thruugh calling 
sequences. However, the following variables must 
be available in scratchpad (locations 0-256 jq ) 

LOCATION 

VARIABLE 

DESCRIPTION 

81 16 

N0RATE 

Number of Rates (1, 2, or 3) 

82 16 

N0L0 

Number of Low Speed A/D channels 

83 16 

LOCNT 

Number of High speed reads per low speed 
read 

84 16 

NOMED 

Number of medium speed channels. 

85 16 

MEDCNT 

Number of high speed reads per medium 
speed read 

86 16 

N0HS 

Number of high speed channels 

8716 

BUFZER 

Starting address for storing data into 
buffer zero. 

LO 

tH 

CO 

CO 

BUFZND 

Ending address for storing data into 
buffer zero. 

8916 

BUF0NE 

| Address for buffer one. 

8Ai6 

NDBUF1 


8Ci6 

MAXBUF 

Data buffer length 

8016 

TP0UT 

Tape output flag 
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LOCATION 

VARIABLE 

8Ei6 

DONCLK 

*o 

lx. 

00 

BGNCHN 

91 16 

D0NMS 

92!, 

DNSEC 

93 16 

DNMIN 

94 16 

DNHRS 

8B 1 6 

BUFUSE 


DESCRIPTION 

Systron Donner time conversion flag. 

Beginning A/D channel number 

Temporary storage for 
Systron Donner times 
(Hours, minutes, seconds, 
and milliseconds) 

Flag indicating which data buffer is 
currently being filled (buffer zero 
or one) . 


Subroutine STOPIT sets the following variables to their initial state. 
BUFEND End of Buffer sentinel set to 0. 

D0NCLK Systron Donner convert flag set to 0 
CNTMED Medium A/D rate internal counter to -1 
CNTLO Low A/D rate internal counter to -1 

AOCNUM Number of A/D chanels to read to -1 


4. 3. 3. 3. 5 


4. 3. 3. 3. 6 


4. 3. 3. 3. 7 


Output Data - CLOCK outputs the start codes to the 
A/D convertor and sets up the automatic input loca- 
tions for the digital input values. 

Storage Required - The CLOCK program requires 66 i 6 or 
102 10 locations. This includes 14 i 0 locations 
required by ST0P1 f . 

Description - When a clock interrupt occurs, the 
CLOCK subroutine responds by clearing the interrupt. 
Since the time period between high speed A/D reads 
was used as a clock count, each interrupt signals 
the start time for high speed data acquisition. 

The low and medium speed counters are checked (and 
incremented) to determine if the low and medium 
speed channels are to be read on this particular 
interrupt. After the number of channels to be read 
is determined, CLOCK checks the buffer to insure 
that all the data to be converted will fit. If so, 
the acquisition is started. If the data will not 
fit, the buf feres are switched and the A/D convertor 
is started. The systron donner times are read and 
the donner time convert and tape output flags are set. 
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In either case, control is returned to the inter- 
rupted program. 

When the test is complete, STOP IT is called by an 
external program. This subroutine stops the A/D 
convertor and the clock to prevent future conver- 
sions. STOPIT also resets counters interal to CLOCK. 
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4. 3. 3. 4.0 


Subroutine TEST 


4. 3. 3. 4.1 The purpose of this subroutine is to control the 

background tasks during a data acquisition. 

4.3. 3.4.2 The calling sequence is JST * TEST 

4. 3. 3. 4. 3 Software/Software interfaces. 


Subroutine TEST calls the following programs. 

CLOCK illed for the first data acquisition 

All additional acquisitions are done with clock 
interrupts. 

TDUMP Called to start the data transfer to magnetic tape 

anc to write end of files on the tape. The A register 

contains the minus number of words to be transfered and 
the X register contains the buffer address minus one. 
The A register is loaded with -1 for writing end of 
fi les . 

UNPACK Called to convert the Systron Donner times from packed 

BCD to binary numbers. The A register contains the 

number to be converted. 

ST0PIT Called to stop the data acquisition and reset the 

CLOCK routine. 

OTL Called to output messages to the teletype 

ODEC Called to convert a binary number to decimal and 

output it to the teletype 

CRLF Called to output a carriage return-line feed. 


4. 3. 3. 4. 4 Input Data - There are no data passed as calling argu- 

ments. However, the following variables must be 
available in scratchpad (locations 0-255! 6 ) 

LOCATION VARIABLE DESCRIPTION 

8Bi 6 BUFUSE Buffer in use flag. 

8Eie D0NCLK Donner Clock time conversion flag. 

80i6 CLCNT Clock count (time between interrupts) in 

microseconds . 
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4. 3. 3. 4. 5 


Output Data - The following messages are output 
from TEST. 


1 - DONE - ERRORS — PARITY XXXX TIMING. XXXX 

2 - XXXX RECORDS WRITTEN. 

3 - RATE T 00 FAST! 

Message one and two give numeric values for the parity and timing errors and 
the number of records written. 

Message three indicates that one buffer is full before the previous buffer 
has been completely output. This is an abort condition. 


4. 3. 3.4.6 Storage required. TEST Requires AB 16 or 171 10 locations. 


4. 3. 3. 4. 7 Description. When TEST is called, the appropriate 

counters (parity error, timing, tape output, etc.) are set to zero and the 
appropriate I/O devices are set up. The clock is started and subroutine 
CLOCK is called to start LLe data acquisi tion. TEST then enters a loop. 

The first step of the loop is to check for a buffer ready for output and 
tape not busy. If a buffer is ready, it is output with a call to TDUMP. 

If not, then the second step of the loop is executed. The second step 
checks for Systron Donner times to be converted. If there are times to 
be converted, then the routine UNPACK is called. After the times are 
unpacked and stored for output, the test stop times are checked to see 
if the test has exceeded the time. If so, the test is stopped. If not, 
the third step of the loop is executed. This step looks for an 'X' char- 
acter from the teletype. If an X has been input, the test is stopped. If 
an X has not been input, the loop is repeated. 

After the test has been completed, subroutine TEST writes two end of files 
on the tape and backspaces over the second. It then informs the user of 
any tape errors and of the total number of records written on tape. Control 
is then returned to the calling program. 
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80, 6 

TP0UT 

Tape output flag. 

91,6 

DONMS 

| 

92 , 6 

DNSEC I 

[ Temporary storage for Systron 

93,6 

DNM1N , 

Donner times. 

94,6 

DNHRS 

I 

99,6 

STPHR 


9A, 6 

ST PM IN 1 

Test stop times in hours, minutes, 

9B, 6 

STPSEC | 

[ seconds, and milliseconds. 

9C, 6 

ST0PMS 

I 

9D,6 

9E,6 

BUFHDI [ 
BUFHDO ) 

Address of the data buffers. 

9F,s 

MBUFL6 

Minus data buffer length. 

Al,6 

PAR 

Tape parity error counter. 

u> 

CsJ 

< 

TIMG 

Tape timing error counter. 

A6,6 

TAPEND 

Tape end of operation flag. 

In addition. 

the addresses of all the external subroutines are stored in 

scratchpad . 



LOCATION 

ROUTINE 


3B, 6 

CL0CK 


71,6 

CRLF 


72,6 

0TL 


73,6 

UNPACK 


75, 6 

TDUMP 


76,e 

0DEC 


7B, 6 

ST0PIT 
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4. 3. 3. 5.0 

Subroutine REST0R 


4. 3. 3. 5.1 

The purpose of this subroutine is to reset all I/O 
devices and to restore all scratchpad variables to 
their initial state. 

4. 3. 3. 5. 2 

Calling Sequence. 

The calling sequence is JST*REST0R 

4. 3. 3. 5. 3 

Software/Software 

Interfaces - none 

4. 3. 3. 5. 4 

Input Data 

- none 


4. 3. 3. 5. 5 

Output Data - The following list of scratchpad 
variables are restored to the indicated values. 

LOCATION 

VARIABLES 


VALUE 

80 16 

CLCNT 


0 

81,6 

NORATE 


0 

82,6 

N0L0 


0 

83,6 

L0CNT 


0 

84,6 

N0MED 


0 

85,6 

MEDCNT 


0 

86,6 

N0HS 


0 

87,6 

BUFZER 


77C i 6 

88,6 

BUFZND 


B64 i 6 

89 , 6 

BUF0NE 


B7D i 6 

8A, 6 

NDBUF1 


F65 i 6 

8C,6 

MAXBUF 


1025 

8F, 6 

BGNCHN 


0 

90,6 

ADCRD 


0 

95,6 

STHP 


0 

96,6 

STMIN 


0 

97,6 

STSEC 
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LOCATION 

VARIABLE 

VALUE 

98 16 

ST RMS 

0 

99,6 

STPHR 

0 

9A 1S 

STPMIN 

0 

9B 1S 

STPSEC 

0 

9C , 6 

STOPMS 

0 

9D 1S 

BUFHDI 

B64 1 6 

9E , 6 

BUFHDO 

763 1 6 

9F 16 

MBUFLG 

-1025 

AO, s 

MXBDUF 

1000 

A3,6 

HDBFLG 

25 

A6, 6 

TAPEND 

1 


The twenty five header words in each data buffer are also filled with zeros. 


4. 3. 3. 5. 6 Storage Required - This routine requires 4 1 1 6 or 65 i o 

locations . 

4. 3. 3. 5. 7 Description - The I/O devices are reset first. Then 
the appropriate values are restored by using constants defined in REST0R 
or by constants defined by load A register immediate instructions. 
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4. 3. 3. 6.0 


Subroutine EINS, ZWEI, and DREI 


4. 3. 3. 6.1 


The purpose of these subroutines is to input Data 
generic to the number of A/D conversion rates. 


4. 3. 3. 6. 2 Calling Sequence - The calling sequences are - 

I JST * DREI 

I + 1 N0P 

1+2 NOP 

I + 3 RETURNS HERE 

I JST * ZWEI 

I + 1 RETURNS HERE 

I JST * EINS 

I + 1 N0P 

I + 2 RETURNS HERE 

4.3. 3.6.3 Software/5.oftware Interfaces - Subroutine FETCH is 

cal led to i;ype the data request and read the input 

value. 


4. 3. 3.6.4 Input Data - The only data input is through the tele- 

type and is returned in the A register by FETCH. 

The following input values are required by the indicated 
entry point. 

ENTRY POINTS VARIABLES REQUIRED 


EINS, ZWEI, DREI 
ZWEI, DREI 
DREI 

ZWEI, DREI 
DREI 


N0HS number of High speed A/D channels. 

N0MED number of medium speed a/D channels 

N0L0 number of low speed A/D channels 

MEDCNT number of times to read high speed channels 
before reading medium speed channels 

L0CNT number of times to read high speed channels 
before reading low speed channels. 


4. 3. 3.6.5 Output Data - The Data which are teletype inputs are 

stored in scratchpad (locations 0-255). In addition, 
descriptive messages requesting input are passed to 
FETCH to be output to the teletype. 
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4. 3. 3. 6. 6 


S torage Required - The total storage required by 
the three entry points is 71 i e or 1 1 3 i 0 locations. 


4.3. 3.6. 7 Description - Each subroutine only types out that 

information which is generic to its specific rate. The X register is 
loaded with the address of the message to be output. Subroutine FETCH 
is called to type out the message and read the input number. If the 
number is in error ( i . e . , couldn't be converted from ASCII), FETCH returns 
to the next instruction past the call. The message loop is then repeated 
If the value is valid, it is stored and the next call to FETCH is executed. 
All the entry points return control to the calling program. 


H-24 



Subroutine HEADER 


4. 3. 3. 7.0 


4. 3. 3. 7.1 

4. 3. 3. 7. 2 

4. 3. 3. 7. 3 


4. 3. 3. 7. 4 

9Di6 BUFHDI 

9E i 6 BUFHDO 

A4 1 6 BUFTTY 


The purpose of this subroutine is to accept user 
input data for the data buffer header words 6-25. 

The Calling sequence is JST * HEADER 

Software/Software Interfaces - 

Subroutine HEADER calls: CRLF to output carriage 

return-line feeds; OTT to output a single character 

from the A register; 0TL to output a message; and 

FETCH to output a message and to read a number for 

storage. 

Input Data - The following data must be available 
in scratchpad - 

Address of the buffer header for 
buffer one and zero. 

Address of the teletype input buffer. 


The only other input data is from the teletype for storage in the headers. 


4. 3. 3. 7. 5 Output Data - up to 20 words of information may be 

stored in the buffer headers. Output message addresses 
are passed to FETCH and 0TL. A single character to 
be output is transferred to 0TT via the A register. 

4. 3. 3. 7. 6 Storage Required - HEADER requires 57 or 87 locations. 

4. 3. 3. 7. 7 Description - HEADER calls 0TL to output the message: 

HEADER DATA INPUT TO TERMINATE, INPUT 'A/'. The subroutine FETCH is called 

to write the word number of the header in to which the next input value will 

go. If the value input is in error, the FETCH is called to print the mes- 
sage again. If the value is valid, it is stored in the headers. The header 
word pointers are incremented and FETCH is called again with the next word 
number. The process is repeated until either 20 values are input or the 

quantity "A/" is input. HEADER then returns control to the calling program. 
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4.3. 3.8.0 


4. 3. 3. 8.1 


4. 3. 3. 8. 2 


4. 3. 3. 8. 3 


4. 3. 3. 8. 4 

4. 3. 3.8.5 

4. 3. 3. 8. 6 


Subroutine PFAIL 

The purpose of the subroutine is to provide an 
orderly shutdown of the CF16A during a power 
fail ure. 

There is no calling sequence since this is an 
interrupt subroutine. 

Software/Software Interfaces 

This subroutine doesn't call any other subroutines. 
However, it does store a jump to the restart rou- 
tine at the interrupt location, 00i 6 . 

Input Data - None 

Output Data - An instruction is stored at the 
restart interrupt location 

Storage Required - This routine requires 7 16 or 7 10 
locations . 


4. 3. 3. 8. 7 Description - When a power failure is detected, an 

interrupt is generated. This subroutine responds to that interrupt by 
halting an tape I/O. It also stores a jump to the restart routine at 
location 0. The interrupts are enable and the CFIGA is halted. 
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4. 3. 3. 9.0 


Subroutine PWR0N 


4. 3. 3. 9.1 The purpose of PWR0N is to insure an orderly restart 

after a power failure. 


4. 3. 3. 9. 2 There is no calling sequence since this routine 

responds to an interrupt at location 0. 


4. 3. 3. 9. 3 Software/Software Interfaces - This subroutine calls 

REST0R to reset the system variables to their ori- 
ginal state. It also calls ST0PIT to reset the 
CLOCK data acquisition routine to its original state. 


4. 3. 3. 9. 4 Input Data - None 


4. 3. 3. 9. 5 Output Data - This subroutine stores a jump to the 

spurious interrupt handler at location 0. 


4. 3. 3.9.6 Storage Required - This routine requires 9i 6 locations. 


4. 3. 3. 9. 7 Description - When a restart interrupt occurrs at 

location 0, this routine calls REST0R, ST0PIT, and 
then restores the spurious interrupt subroutine call 
at location 0. (When a device requests an interrupt 
and does not respond with an address, the CF16A 
interrupts to location 0.) PWR0N then begins execu- 
tion at location 111 is. the start of GDAS. 
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4.3.3.10.0 

4.3.3.10.1 

4.3.3.10.2 

4.3.3.10.3 

4.3.3.10.4 

4.3.3.10.5 

4.3.3.10.6 

4.3.3.10.7 


Subroutine BADRPT 

The purpose of this subroutine is to alert the 
user that a spurious interrupt has occurred. 

There is no calling sequence since this is an 
interrupt subroutine. 

Sof tware/Sof tware Interfaces - This subroutine 
calls 0TL to output the message - 

"SPURIOUS INTERRUPT! CYCLE TO CONTINUE." 

Input Data - None 

Output Data - The ASCII character ' is loaded 
into the A register before calling 0TL. 

Storage Required - This subroutine requires 1 C i 6 
or 27 10 locations. 

Description - When a device requests an interrupt 
and does not respond with an interrupt address, the 
CF16A interrupts to location 0. This subroutine 
issues a halt 1/0 to the tape controller and then 
outputs a message to alert the user of the interrupt 
and halts. Although the user can restart at the 
point of interruption, this interrupt usually signi- 
fies a hardware fault and is likely to recur until 
the hardware is repaired. 
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4.3.3.11 .0 


4.3.3.11.1 


4.3.3.11.2 


4.3.3.11 .3 


4.3.3.11.4 
4.3.3.11 .5 


4.3.3.11.6 


Subroutine TPCHK 

The checks the tape unit to insure that a tape 
is mounted and not wirte protected. 

Calling sequence - the calling sequence is 
JST * TPCHK. 

Sof tware/Software Interfaces - TPCHK calls CRLF 
to issue a carriage return-line feed to the tele- 
type. Subroutine 0TL is called to output the 
message. 

"TAPE NOT ON-LINE." 

" TAPE IS WRITE PROTECTED." 

Input Data - None 

Output Data - The message termination character is 
transfered to 0TL via the A register. The address of 
the message to be output is transfered as an argume^? 
immediately after the call. 

LDA termination Character 
JST * 0TL 

DATA message address. 

Storage Required - This subroutine requires 30 16 or 
48 i o locations. 


4.3.3.11.7 Description - Using the sense instruction, TPCHK 

checks to insure that the tape is on-line and not 
write protected. If either of these conditions are 
not met, a message is output. TPCHK then enters a 
loop until the error condition is corrected. After 
the error conditions are corrected, TPCHK returns to 
the cal ling program. 
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4.3.3.12.0 

4.3.3.12.1 

4.3.3.12.2 

4.3.3.12.3 

4.3.3.12.4 

4.3.3.12.5 

4.3.3.12.6 

4.3.3.12.7 


Subroutine ADCHk 

Subroutine ADCHK checks the A/D convertor for 
automatic mode. 

Calling Sequence - The calling sequence is 
JST *ADCHK 

Software/Software Interfaces - This subroutine 
calls CRLF to output a carriage return- line feed 
and 0TL to output the message 

"ADC IN MANUAL M0DE." 

Input Data - None 

Output Data - The termination character of the message 
is transferred to 0TL via the A register. The address 
of the message is transferred as a calling argument. 

LDA Termination Character 
JST * 0TL 

DATA Message Address 

Storage Required - ADCHK required 1 6 i 6 or 22 i 0 loca- 
tions . 

Description - Using the sense instruction, ADCHK 
determines whether or not the A/D convertor is in 
manual mode. If so, ADCHK outputs a message to 
the user and waits for the A/D converter to be 
put in automatic mode. ADCHK returns to the call- 
ing program. 
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